Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,239 A→T 5.3% intergenic (+332/+40) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2390AT 5.3% 1113.8 / 37.1 321intergenic (+332/+40)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (148/156);  new base T (10/7);  total (158/163)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 330 total aligned reads displayed.

CTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACG  >  NC_000913/1224091‑1224387
                                                                                                                                                    |                                                                                                                                                    
cTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgaga                                                                                                                                                      <  1:1825177/149‑1 (MQ=255)
 tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                     >  3:797086/1‑149 (MQ=255)
 tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                     >  2:2364873/1‑149 (MQ=255)
 tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                     >  4:2413654/1‑149 (MQ=255)
 tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGGGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                     >  1:2295103/1‑149 (MQ=255)
       aGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                                                                                               <  2:767706/149‑1 (MQ=255)
         gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             <  1:1668499/149‑1 (MQ=255)
         gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             <  1:2195153/149‑1 (MQ=255)
            tAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                                                                          <  1:1509530/149‑1 (MQ=255)
             aaaCATGTTGTGTTGTGCGTGGGAACAGTCTTAAGGTGGAAGTGGAGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTTGGCGCGTTACCACAACGTGGATAAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTgtgttg                                                                                                                                         >  4:453505/1‑144 (MQ=255)
              aaCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                                                                        >  3:1387547/1‑149 (MQ=255)
               aCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                                                                       <  1:46129/149‑1 (MQ=255)
               aCATTTTGTTTTGTGCGTGGAAACAGCCATAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGTATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                                                                       <  2:387553/149‑1 (MQ=255)
                  ttttgttttgTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                                                                    >  3:1655838/1‑149 (MQ=255)
                  ttttgttttgTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                                                                    >  2:1121196/1‑149 (MQ=255)
                      gttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                                <  2:2082063/149‑1 (MQ=255)
                            tgCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                                                                          <  2:1386217/149‑1 (MQ=255)
                                 gggAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGcc                                                                                                                                            >  2:1251811/1‑126 (MQ=255)
                                  ggAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                                                                           >  4:2263553/1‑126 (MQ=255)
                                  ggAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                                                                    >  2:2101151/1‑149 (MQ=255)
                                       aGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagttt                                                                                                                                    >  2:2232996/1‑124 (MQ=255)
                                        gCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                                                                              <  1:1056208/149‑1 (MQ=255)
                                         ccTTAAGGGGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGAGGATGAGAGGGATGAAAAACTCACGGCAGAGATAACTCTGCCTGGAAGATAAATGCGCTTTTACAGCGGGCtt                                                                                                             >  4:119677/1‑149 (MQ=255)
                                          cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTa                                                                                                            <  4:1757071/149‑1 (MQ=255)
                                           ttAGGGTGTAAAGGGGGAGGTGGAAATAGCCGTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGAGGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                                                                           >  4:2196672/1‑149 (MQ=255)
                                           ttAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                                                                           >  3:1999973/1‑149 (MQ=255)
                                           ttAAGGTGTAAAGGGGGAGGCGGAAATAGCAATGAGGAGTATCAGCAAGACTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                                                                           >  2:1875119/1‑149 (MQ=255)
                                             aaGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                                                                                                         <  1:2101151/149‑1 (MQ=255)
                                              aGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                                                                        <  1:777834/149‑1 (MQ=255)
                                               ggTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                       <  3:1220618/149‑1 (MQ=255)
                                               ggTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCAGGCTTATTTCa                                                                                                       <  4:1708723/149‑1 (MQ=255)
                                               ggTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGACGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                       <  2:925314/149‑1 (MQ=255)
                                                gtgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                      <  1:717983/149‑1 (MQ=255)
                                                gtgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                      <  2:1566676/149‑1 (MQ=255)
                                                 tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                                                                           >  3:1058585/1‑111 (MQ=255)
                                                 tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                                                                           <  4:1058585/111‑1 (MQ=255)
                                                  gtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGct                                                                                                    >  4:1886949/1‑149 (MQ=255)
                                                      aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                                <  1:290389/149‑1 (MQ=255)
                                                       gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                                                                                               >  3:992414/1‑149 (MQ=255)
                                                       gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                                                                                               <  1:627885/149‑1 (MQ=255)
                                                          ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCt                                                                                            >  4:535005/1‑149 (MQ=255)
                                                             ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                                                                         <  1:2206767/149‑1 (MQ=255)
                                                                ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                                                                     >  4:841346/1‑134 (MQ=255)
                                                                ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                                                                      <  1:1897335/149‑1 (MQ=255)
                                                                ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                                                                      >  2:2488801/1‑149 (MQ=255)
                                                                 gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTa                                                                                     <  2:1321526/149‑1 (MQ=255)
                                                                  aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                                                                    <  4:802169/149‑1 (MQ=255)
                                                                      aGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                                                                                <  3:2459930/149‑1 (MQ=255)
                                                                        cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc                                                                              <  4:1999973/149‑1 (MQ=255)
                                                                        cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc                                                                              <  2:45913/149‑1 (MQ=255)
                                                                        cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc                                                                              >  1:885749/1‑149 (MQ=255)
                                                                          aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                           >  1:1253482/1‑86 (MQ=255)
                                                                            gaggatTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                                                                          <  1:1743718/149‑1 (MQ=255)
                                                                             aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCa                                                                         <  2:1946013/149‑1 (MQ=255)
                                                                               gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGc                                                                       >  1:1407614/1‑149 (MQ=255)
                                                                                      aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                                >  1:1132824/1‑85 (MQ=255)
                                                                                      aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa                                                                >  1:1759278/1‑149 (MQ=255)
                                                                                       gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                                                                         >  4:2118951/1‑91 (MQ=255)
                                                                                       gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                                                                         <  3:2118951/91‑1 (MQ=255)
                                                                                       gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                                                                   <  3:395865/129‑1 (MQ=255)
                                                                                           gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa                                                           <  1:1611326/149‑1 (MQ=255)
                                                                                                 tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                                                                       <  2:1767652/67‑1 (MQ=255)
                                                                                                   gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc                                                   >  4:2091700/1‑149 (MQ=255)
                                                                                                        tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                                                                                                          >  1:2455587/1‑89 (MQ=255)
                                                                                                           aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttttt                                           <  2:548/149‑1 (MQ=255)
                                                                                                           aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttttt                                           <  2:1681995/149‑1 (MQ=255)
                                                                                                           aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttttt                                           <  3:1106875/149‑1 (MQ=255)
                                                                                                            ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                                                                                               <  2:644923/96‑1 (MQ=255)
                                                                                                            ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTg                                          <  1:2029982/149‑1 (MQ=255)
                                                                                                             cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGc                                         >  2:486285/1‑149 (MQ=255)
                                                                                                              acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCt                                        >  3:1514426/1‑149 (MQ=255)
                                                                                                              acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCt                                        <  3:1318478/149‑1 (MQ=255)
                                                                                                                  cGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAg                                    <  4:1425492/149‑1 (MQ=255)
                                                                                                                    tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCt                                  <  3:174219/149‑1 (MQ=255)
                                                                                                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGt                                <  3:1877764/149‑1 (MQ=255)
                                                                                                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGt                                >  3:1428763/1‑149 (MQ=255)
                                                                                                                        tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTAc                              <  1:955260/149‑1 (MQ=255)
                                                                                                                        tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTAc                              <  1:1714854/149‑1 (MQ=255)
                                                                                                                        tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTAc                              <  4:1634762/149‑1 (MQ=255)
                                                                                                                               gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                                                                             >  4:957898/1‑63 (MQ=255)
                                                                                                                               gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc                                                    <  2:2426111/120‑1 (MQ=255)
                                                                                                                               gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCGCAAGCTGTACGGTTAcc                       <  3:2259357/149‑1 (MQ=255)
                                                                                                                                  gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtc                    >  4:688403/1‑149 (MQ=255)
                                                                                                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtct                   >  4:1321440/1‑149 (MQ=255)
                                                                                                                                       aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAgg               <  3:424394/149‑1 (MQ=255)
                                                                                                                                        ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGa              >  3:957568/2‑149 (MQ=255)
                                                                                                                                        cTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTCATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTCTTTCTCAAGCTGTACGGTTACCATCTCAGGa              <  4:1522214/149‑1 (MQ=255)
                                                                                                                                         tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGAt             >  4:1168629/1‑149 (MQ=255)
                                                                                                                                          cAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGc                                         <  3:2321822/120‑1 (MQ=255)
                                                                                                                                            aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCaa          >  3:1084983/1‑149 (MQ=255)
                                                                                                                                            aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCaa          <  2:2186752/149‑1 (MQ=255)
                                                                                                                                             ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAt         >  3:2421930/1‑149 (MQ=255)
                                                                                                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCTAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAtt        >  1:1798437/1‑149 (MQ=255)
                                                                                                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc                                                      <  2:1489529/103‑1 (MQ=255)
                                                                                                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc                                                      >  1:1489529/1‑103 (MQ=255)
                                                                                                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTg      <  2:1844622/149‑1 (MQ=255)
                                                                                                                                                 gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGt     <  3:852567/149‑1 (MQ=255)
                                                                                                                                                  agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTa    <  3:535005/149‑1 (MQ=255)
                                                                                                                                                   gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTAc   <  1:1213632/149‑1 (MQ=255)
                                                                                                                                                   gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCACCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTAc   >  1:512798/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACG  >  NC_000913/1224091‑1224387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: