Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,547,305 T→C 5.6% V737A (GTG→GCG)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,547,3050TC 5.6% 210.2 / 2.8 72V737A (GTG→GCG) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTG  >  NC_000913/4547157‑4547453
                                                                                                                                                    |                                                                                                                                                    
gatgatATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGt                                                                                                                                                      >  4:270297/1‑149 (MQ=255)
    atatTAAGCAGCTATATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGGGCTTGGTAGAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGGGCg                                                                                                                                                   <  2:93547/148‑1 (MQ=255)
     tatTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTa                                                                                                                                                 >  2:349350/1‑149 (MQ=255)
          agcagcTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGa                                                                                                                                             <  4:225986/148‑1 (MQ=255)
             agcTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTGAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTgg                                                                                                                                         >  7:51416/1‑149 (MQ=255)
              gcTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGc                                                                                                                                        <  2:196506/149‑1 (MQ=255)
                tctATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGt                                                                                                                                      >  5:318316/1‑149 (MQ=255)
                   aTTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTAGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGt                                                                                                                                   >  2:55154/1‑149 (MQ=255)
                      aCGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTAt                                                                                                                                >  2:97079/1‑149 (MQ=255)
                         taGTCAGCTGTGGGGTACTGGCTCATGCAAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCCAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGATTATGcc                                                                                                                             <  6:85335/148‑1 (MQ=255)
                            tCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtg                                                                                                                          >  4:116297/1‑149 (MQ=255)
                                   tGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCtt                                                                                                                   <  6:166547/149‑1 (MQ=255)
                                     gggTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTa                                                                                                                  <  1:335128/148‑1 (MQ=255)
                                     gggTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTa                                                                                                                  <  7:293472/148‑1 (MQ=255)
                                     gggTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTAt                                                                                                                 <  1:215626/149‑1 (MQ=255)
                                        tACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGCAAACCAGACGGGGGTGCGTACCGACTGGCGGGGTGATGCCGTGCTGCCTTATGc                                                                                                               >  5:313575/1‑148 (MQ=255)
                                        tACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGc                                                                                                               >  8:259323/1‑148 (MQ=255)
                                        tACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCCACTGGCGTGGTTATGCCGTGCTGCCTTATGc                                                                                                               >  5:38360/1‑148 (MQ=255)
                                             gCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTTGAGTGCTTATGACGTGCTGCATTATCCCACTg                                                                                                          >  5:5884/1‑148 (MQ=255)
                                               tCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGCGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                                                        >  5:204660/1‑148 (MQ=255)
                                               tCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGCGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                                                        >  5:209022/1‑148 (MQ=255)
                                                cATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAata                                                                                                      >  8:343818/1‑149 (MQ=255)
                                                 aTGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                                                                     <  3:270297/149‑1 (MQ=255)
                                                    ccAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                                                                  <  6:318316/149‑1 (MQ=255)
                                                    ccAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCg                                                                                                   <  1:55154/148‑1 (MQ=255)
                                                     cAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                                                                  >  8:74150/1‑148 (MQ=255)
                                                     cAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAAACAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                                                                  >  2:276536/1‑148 (MQ=255)
                                                      aaTGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTAc                                                                                                                                                <  7:147307/101‑1 (MQ=255)
                                                      aaTGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTAc                                                                                                                                                >  8:147307/1‑101 (MQ=255)
                                                        tGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                                                              <  6:64443/149‑1 (MQ=255)
                                                          gCGTAACGCTGTGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTCGCGTGGTTATGCCGTGCTGCCTTATGACACTGAATATCGGGaaaa                                                                                             >  1:50402/1‑148 (MQ=255)
                                                             tAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATaga                                                                                         <  1:349350/149‑1 (MQ=255)
                                                              aaCGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                                                        >  5:261020/1‑149 (MQ=255)
                                                                  cTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                                                                    >  7:52774/1‑149 (MQ=255)
                                                                        cAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTTGCGCTGGa                                                                              >  8:67839/1‑149 (MQ=255)
                                                                         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGa                                                                              >  1:129069/1‑148 (MQ=255)
                                                                           ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                                                           >  2:96626/1‑149 (MQ=255)
                                                                             gTTAAACGATACGGTGGTGCTTGTTGAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                                                                         <  8:51416/149‑1 (MQ=255)
                                                                                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGCGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                                                                    <  6:209022/149‑1 (MQ=255)
                                                                                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGCGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                                                                    <  6:204660/149‑1 (MQ=255)
                                                                                    gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGAc                                                                                                                                            <  8:23278/75‑1 (MQ=255)
                                                                                    gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGAc                                                                                                                                            >  7:23278/1‑75 (MQ=255)
                                                                                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                                                               >  2:169893/1‑147 (MQ=255)
                                                                                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                                                               <  1:169893/147‑1 (MQ=255)
                                                                                                cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGt                                                      >  2:73332/1‑149 (MQ=255)
                                                                                                cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGt                                                      <  1:73332/149‑1 (MQ=255)
                                                                                                 ttGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc                                                     >  7:77540/1‑149 (MQ=255)
                                                                                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATaga                                                                                         >  3:311421/1‑112 (MQ=255)
                                                                                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATaga                                                                                         <  4:311421/112‑1 (MQ=255)
                                                                                                    ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAt                                                  <  4:268574/149‑1 (MQ=255)
                                                                                                       aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt                                                 >  8:324773/1‑147 (MQ=255)
                                                                                                          cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa                                             >  6:280251/1‑148 (MQ=255)
                                                                                                              tGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                                                        >  7:316118/1‑101 (MQ=255)
                                                                                                              tGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                                                        <  8:316118/101‑1 (MQ=255)
                                                                                                                 cgcAAAAGATTCAAAAGTCGAATACCAGACGTGGGTGCGTAGCGACTGTCGTGGTTATGCCTTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATATCGCgg                                     <  6:5884/149‑1 (MQ=255)
                                                                                                                  gcAAAAGATGCAAAAGTCGAAAACCAGCCGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGt                                    <  2:282037/149‑1 (MQ=255)
                                                                                                                  gcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGGGGt                                    >  5:88274/1‑149 (MQ=255)
                                                                                                                   cAAAAGAAGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGCGAAAATAGAGTGGCGCTGGATACCAATTCCCTGGCTGATAACGTCGATTTAGATAACGCGGtt                                   >  1:163961/1‑149 (MQ=255)
                                                                                                                    aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATCCCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTg                                  >  1:357740/1‑149 (MQ=255)
                                                                                                                    aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTg                                  <  7:343818/149‑1 (MQ=255)
                                                                                                                        gATGCAAAAGTCGAAAACCAGACGGGGGGGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTaa                              <  1:97079/149‑1 (MQ=255)
                                                                                                                         aTGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAAc                             >  8:291259/1‑149 (MQ=255)
                                                                                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACgttgt                        >  2:47080/1‑149 (MQ=255)
                                                                                                                               aaGTCGAAACCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACgttgtt                       <  4:241691/149‑1 (MQ=255)
                                                                                                                               aaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACgttgtt                       <  3:116297/149‑1 (MQ=255)
                                                                                                                                aGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTc                      <  8:52774/149‑1 (MQ=255)
                                                                                                                                aGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTc                      <  2:129069/149‑1 (MQ=255)
                                                                                                                                      atagcAGACGGGGGTGCGTACCTACTGGCGTTGTTATGCCGTGCTGCCTTTTGCCACTGAATATCGGGAAAATAGTGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTc                <  6:38360/145‑1 (MQ=255)
                                                                                                                                            gACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTg             <  6:170137/146‑1 (MQ=255)
                                                                                                                                            gACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTg             >  5:170137/1‑146 (MQ=255)
                                                                                                                                               gggggTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGa       <  7:74150/149‑1 (MQ=255)
                                                                                                                                                ggggTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGAt      <  5:280251/149‑1 (MQ=255)
                                                                                                                                                 gtgTGCGTACCGACTGTCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATCCCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATc     <  1:276536/147‑1 (MQ=255)
                                                                                                                                                    tGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGGTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTg  <  1:96626/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTG  >  NC_000913/4547157‑4547453

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑