Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,877,558 G→C 5.6% intergenic (‑93/+147) yidB ← / ← gyrB DUF937 family protein/DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,5580GC 5.6% 1091.0 / 45.6 302intergenic (‑93/+147)yidB/gyrBDUF937 family protein/DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (149/136);  new base C (3/14);  total (152/150)
Fisher's exact test for biased strand distribution p-value = 5.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 314 total aligned reads displayed.

TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTA  >  NC_000913/3877410‑3877707
                                                                                                                                                    |                                                                                                                                                     
 atTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAgatgat                                                                                                                                                     <  2:3033527/150‑1 (MQ=255)
     ccccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTcc                                                                                                                                                 <  1:2147123/150‑1 (MQ=255)
      cccGCATCGCCTTTCAGACAGGCACCGACACCTTCATCACACAGCCCCATGCTACTCCCCTGTAAAGCGTTACGTGTTGACCCAACCAGTATAGATTTGCGGGATGATAATTGCGGATGGCAGAAATAAAAAGGGCGGAGATGATCTCCg                                                                                                                                                <  2:2252776/150‑1 (MQ=255)
      cccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCg                                                                                                                                                >  1:3412665/1‑150 (MQ=255)
      cccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAAACACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCg                                                                                                                                                <  1:2993230/150‑1 (MQ=255)
       ccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGc                                                                                                                                               <  1:2441383/150‑1 (MQ=255)
        cGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAACAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGcc                                                                                                                                              <  2:3938126/150‑1 (MQ=255)
        cGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGcc                                                                                                                                              <  1:1862907/150‑1 (MQ=255)
         gCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGccc                                                                                                                                             <  2:1808838/150‑1 (MQ=255)
         gCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGccc                                                                                                                                             >  2:2151096/1‑150 (MQ=255)
          cATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCt                                                                                                                                            >  2:3984957/1‑150 (MQ=255)
            tCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCttt                                                                                                                                          >  2:1666608/1‑150 (MQ=255)
              gCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTc                                                                                                                                        >  1:3650314/1‑150 (MQ=255)
               ccTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCt                                                                                                                                       >  2:2445975/1‑150 (MQ=255)
               ccTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCt                                                                                                                                       >  1:1134216/1‑150 (MQ=255)
                cTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAACGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                      <  1:1283331/150‑1 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                     >  1:4105164/1‑150 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                     >  1:3276715/1‑150 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                     <  1:1479591/150‑1 (MQ=255)
                  tacAGAAAGGCACCGACAACTACATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTtat                                                                                                                                    <  1:183621/148‑1 (MQ=255)
                   tCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTtata                                                                                                                                   <  1:3995770/150‑1 (MQ=255)
                     aGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                      <  1:2191480/145‑1 (MQ=255)
                     aGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                      >  2:2191480/1‑145 (MQ=255)
                            cACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGc                                                                                                                          >  2:572240/1‑150 (MQ=255)
                            cACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTGTCTTATAGCTTCTTGc                                                                                                                          >  1:2602827/1‑150 (MQ=255)
                               cGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCgg                                                                                                                       >  1:2175014/1‑150 (MQ=255)
                               cGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCgg                                                                                                                       >  1:3714063/1‑150 (MQ=255)
                               cGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCgg                                                                                                                       >  1:1736581/1‑150 (MQ=255)
                                    aCTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcg                                                                                                                  <  2:1984922/150‑1 (MQ=255)
                                    aCTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcg                                                                                                                  <  2:1151468/150‑1 (MQ=255)
                                     cTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcgg                                                                                                                 <  1:856497/150‑1 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcggc                                                                                                                >  1:472255/1‑150 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcggc                                                                                                                >  1:1539448/1‑150 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGcc                                                                                                                                              >  2:3174048/1‑120 (MQ=255)
                                          tcaAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGa                                                                                                            <  2:2811001/150‑1 (MQ=255)
                                              aaaGCCCCATGAGAATCCCCTGTAAAGCGTTACGTGTTGACCCAAAAAGGATCGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGc                                                                                                        <  1:287012/150‑1 (MQ=255)
                                                aGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCt                                                                                                      <  1:3815579/150‑1 (MQ=255)
                                                  ccccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                    >  1:1678425/1‑150 (MQ=255)
                                                    ccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTaac                                                                                                  >  1:1771592/1‑148 (MQ=255)
                                                    ccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                  >  1:1176722/1‑150 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                 >  1:2289128/1‑150 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                 >  2:3569721/1‑150 (MQ=255)
                                                     cATGAGCATCACCTGGAAAGCGTTACGTGCTGACCCAAAACGTATAGAGTCGCGGGATGATAATTGCGGATTGCAGAAATAAACAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                 <  2:3948385/150‑1 (MQ=255)
                                                      aTGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                <  1:3295030/150‑1 (MQ=255)
                                                       tGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATATGGCGTTCACGCCGCATCCgg                                                                                               <  1:3209215/150‑1 (MQ=255)
                                                        gATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGTGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                              <  2:1477196/150‑1 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                     >  1:2735530/1‑150 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                     >  2:4031673/1‑150 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                     >  1:2034604/1‑150 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGAGGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                     >  1:2529727/1‑150 (MQ=255)
                                                                   gTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                   <  2:1248859/150‑1 (MQ=255)
                                                                   gTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCATGAACGCCTTATCCGGCCTACAAAATCa                                                                                   <  1:618135/150‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                  <  2:701739/150‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                  >  2:2085063/1‑150 (MQ=255)
                                                                          gTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAt                                                                            <  2:258056/150‑1 (MQ=255)
                                                                             aCGTGTTGACCCACAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCa                                                                         <  1:3474763/150‑1 (MQ=255)
                                                                              cGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                                                        <  2:2692555/150‑1 (MQ=255)
                                                                                  ttGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACa                                                                    >  1:2266924/1‑150 (MQ=255)
                                                                                  ttGACCCAAAAACTATAGATTTGCGTGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGCATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACa                                                                    <  2:630599/150‑1 (MQ=255)
                                                                                    gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                      <  2:1611014/130‑1 (MQ=255)
                                                                                    gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                      >  1:1611014/1‑130 (MQ=255)
                                                                                     aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                 <  1:1916398/150‑1 (MQ=255)
                                                                                       ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCa                                                               <  2:931064/150‑1 (MQ=255)
                                                                                       ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCa                                                               >  1:974880/1‑150 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                                                             >  2:4009728/1‑150 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                                                             >  2:1174567/1‑150 (MQ=255)
                                                                                          caaaGTAGAGAGTTGCGGGATGATTAGTGCGGATCGCAGACATAAAACCGGCGGAGCTGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACAAAAGCATGAAAATTCAATCCATTGCAAGa                                                            <  2:3852676/149‑1 (MQ=25)
                                                                                          aaacGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGa                                                            <  2:2551750/150‑1 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                       <  2:3123403/150‑1 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                       <  1:1353072/150‑1 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                       <  2:3650314/150‑1 (MQ=255)
                                                                                                ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAGTGCAAGATTTTCg                                                      >  2:35162/1‑150 (MQ=255)
                                                                                                ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCGTTTCTTATAGCGTCTAGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAGTGCAAGATTTTCg                                                      >  2:4057022/1‑150 (MQ=38)
                                                                                                     tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                 <  1:1449515/150‑1 (MQ=255)
                                                                                                     ttagCGGGATGATAATAGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTAAAGCTTCTTCCAGCATGCGGCCTGAACGCCTTAACCGGCCTACAAAAACATGAAAAATCAATAAAATGCAAGAATTTCGAAgtc                                                 >  1:3915407/4‑148 (MQ=255)
                                                                                                      ttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcc                                                <  2:334630/150‑1 (MQ=255)
                                                                                                       tGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt                                               <  2:3239919/150‑1 (MQ=255)
                                                                                                        gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                              >  1:674414/1‑150 (MQ=255)
                                                                                                          gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAt                                            >  1:4041471/1‑150 (MQ=255)
                                                                                                          gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAt                                            >  2:979336/1‑150 (MQ=255)
                                                                                                              tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc                                        >  1:2781413/1‑150 (MQ=255)
                                                                                                              tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc                                        >  1:2680908/1‑150 (MQ=255)
                                                                                                                atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                      >  2:27396/1‑150 (MQ=255)
                                                                                                                   aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                   <  1:3895341/150‑1 (MQ=255)
                                                                                                                   aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                   <  1:430917/150‑1 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgcg                                  >  2:431114/1‑150 (MQ=255)
                                                                                                                       cGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAt                               >  1:3160245/1‑150 (MQ=255)
                                                                                                                       cGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTAGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAt                               >  2:848089/1‑150 (MQ=255)
                                                                                                                          aTTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAg                            <  2:296091/150‑1 (MQ=255)
                                                                                                                              cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                        >  2:3454495/1‑150 (MQ=255)
                                                                                                                              cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                        <  1:3346017/150‑1 (MQ=255)
                                                                                                                               aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACg                       >  1:3174160/1‑150 (MQ=255)
                                                                                                                                           ggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTa           >  2:203073/1‑150 (MQ=255)
                                                                                                                                           ggCGGAGATGATCTCCGCCATTTTCTTATAGCTTCTGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTa           <  2:1929091/150‑1 (MQ=255)
                                                                                                                                             cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgc         <  2:2530618/150‑1 (MQ=255)
                                                                                                                                               gagaTGATCTCCGCCCTTTTCTTCGCGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACCACCGCAGGAAAATTCAATACCTTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgcgc       <  2:1327409/150‑1 (MQ=21)
                                                                                                                                               gagaTCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgcgc       >  1:2355537/1‑150 (MQ=255)
                                                                                                                                                  atgatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAAGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCAt    <  2:1030588/150‑1 (MQ=255)
                                                                                                                                                    gatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTa  >  2:1420379/1‑150 (MQ=255)
                                                                                                                                                    gatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCCGCGTGAAAGCCATATACGGCCAACCAAAAAATGAAAAATAAATAAAAATCAACATTTTAGTAGGCCTTATAAGCGGAGCCCAACAGGCACGCGCGCATGGTTAGCGcccaaa  >  1:1096129/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTA  >  NC_000913/3877410‑3877707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: