Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,159,137 G→A S670S (AGC→AGT ygiQ ← Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,159,1370GA97.1% 389.0 / 0.1 102S670S (AGC→AGTygiQRadical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (50/49);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 103 total aligned reads displayed.

ACGACGCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACC  >  NC_000913/3158992‑3159283
                                                                                                                                                 |                                                                                                                                                  
acgacgCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGc                                                                                                                                                <  1:117650/150‑1 (MQ=255)
 cgacgCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGcc                                                                                                                                               >  2:942392/1‑150 (MQ=255)
 cgacgCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCa                                                                                                                                              <  1:1703388/151‑1 (MQ=255)
   acgCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAAt                                                                                                                                            >  1:1680470/1‑151 (MQ=255)
           tttGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCt                                                                                                                                    >  1:552857/1‑151 (MQ=255)
            ttGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCtt                                                                                                                                   >  1:1247279/1‑151 (MQ=255)
              gCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCtttt                                                                                                                                 >  2:1700260/1‑151 (MQ=255)
              gCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCttt                                                                                                                                  <  2:679224/150‑1 (MQ=255)
                cGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTg                                                                                                                               <  2:938323/151‑1 (MQ=255)
                     cAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCAt                                                                                                                          >  1:568356/1‑151 (MQ=255)
                     cAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCAt                                                                                                                          >  2:10063/1‑151 (MQ=255)
                       ggCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCg                                                                                                                        >  2:449110/1‑151 (MQ=255)
                       ggCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCg                                                                                                                        <  2:600682/151‑1 (MQ=255)
                        gCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGc                                                                                                                       <  1:722420/151‑1 (MQ=255)
                          gTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCtt                                                                                                                     <  2:361884/151‑1 (MQ=255)
                            cTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                   >  1:874471/1‑151 (MQ=255)
                             tGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                   >  2:493180/1‑150 (MQ=255)
                             tGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                   >  1:619596/1‑150 (MQ=255)
                             tGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                   >  2:1665979/1‑150 (MQ=255)
                              gACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                   >  2:1216716/1‑149 (MQ=255)
                              gACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCa                                                                                                                  <  1:350867/150‑1 (MQ=255)
                                cGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAgc                                                                                                                <  1:483475/150‑1 (MQ=255)
                                        tgCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCgg                                                                                                       >  1:1064926/2‑151 (MQ=255)
                                         gCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCgg                                                                                                       >  1:399731/1‑150 (MQ=255)
                                         gCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGa                                                                                                      <  2:1082261/151‑1 (MQ=255)
                                              cGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAAc                                                                                                 >  2:660145/1‑151 (MQ=255)
                                                gCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGTGACTAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAAACGCGACGACTACCAATCAGATGCTTTTTGACCATCGCTTCCAGCGCCTGGCGGATGAACgg                                                                                               <  1:699509/151‑1 (MQ=255)
                                                  gTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGcc                                                                                             >  1:1047563/1‑151 (MQ=255)
                                                        ttCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTg                                                                                       >  1:1523854/1‑151 (MQ=255)
                                                         tCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGc                                                                                      >  2:1171666/1‑151 (MQ=255)
                                                         tCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGc                                                                                      >  1:666440/1‑151 (MQ=255)
                                                          cGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGcc                                                                                     >  2:1659884/1‑151 (MQ=255)
                                                          cGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGc                                                                                      <  1:788814/150‑1 (MQ=255)
                                                          cGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGc                                                                                      <  1:906170/150‑1 (MQ=255)
                                                            tCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCgg                                                                                   >  2:1159057/1‑151 (MQ=255)
                                                             cAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGa                                                                                  <  2:745062/151‑1 (MQ=255)
                                                                         ggTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGc                                                                      <  2:1537290/151‑1 (MQ=255)
                                                                         ggTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGc                                                                      >  1:1634661/1‑151 (MQ=255)
                                                                          gTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGc                                                                      >  2:1224155/1‑150 (MQ=255)
                                                                            aTTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaac                                                                   >  2:1065125/1‑151 (MQ=255)
                                                                            aTTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaa                                                                    <  1:1033767/150‑1 (MQ=255)
                                                                             ttGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGGTCGTGGTAACGcaac                                                                   <  1:1486860/150‑1 (MQ=255)
                                                                              tGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaaca                                                                  >  2:292712/1‑150 (MQ=255)
                                                                              tGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaaca                                                                  >  2:1155423/1‑150 (MQ=255)
                                                                                 ggTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGc                                                              <  1:453846/151‑1 (MQ=255)
                                                                                  gTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCt                                                             <  2:1076528/151‑1 (MQ=255)
                                                                                    tCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCttt                                                           >  2:1653967/1‑151 (MQ=255)
                                                                                     cTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCttt                                                           >  1:182012/1‑150 (MQ=255)
                                                                                       ggCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTAt                                                         <  2:117660/150‑1 (MQ=255)
                                                                                       ggCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATg                                                        >  2:1453476/1‑151 (MQ=255)
                                                                                         cgacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCa                                                      <  1:1489755/151‑1 (MQ=255)
                                                                                          gacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCa                                                      >  1:1075749/1‑150 (MQ=255)
                                                                                          gacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCa                                                      >  1:1345083/1‑150 (MQ=255)
                                                                                           acgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCaaa                                                    <  2:35216/151‑1 (MQ=255)
                                                                                           acgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCaaa                                                    <  1:1700144/151‑1 (MQ=255)
                                                                                           acgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCaa                                                     >  1:883619/1‑150 (MQ=255)
                                                                                            cgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAac                                                   >  1:1107019/1‑151 (MQ=255)
                                                                                             gaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAacg                                                  <  1:702222/151‑1 (MQ=255)
                                                                                                    aCGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTg                                            >  2:537172/1‑150 (MQ=255)
                                                                                                     cGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGtt                                          >  1:1141948/1‑151 (MQ=255)
                                                                                                     cGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGtt                                          <  1:999378/151‑1 (MQ=255)
                                                                                                       cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTg                                        <  2:609865/151‑1 (MQ=255)
                                                                                                       cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTg                                        <  2:607456/151‑1 (MQ=255)
                                                                                                        aTCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGt                                       >  1:828571/1‑151 (MQ=255)
                                                                                                         tctcTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTc                                      >  1:467552/1‑151 (MQ=255)
                                                                                                           tctTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGc                                    >  1:966293/1‑151 (MQ=255)
                                                                                                           tctTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGc                                    >  2:1637043/1‑151 (MQ=255)
                                                                                                            ctTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGt                                       <  1:901779/147‑1 (MQ=255)
                                                                                                            ctTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCg                                     <  1:915924/149‑1 (MQ=255)
                                                                                                              tCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGccc                                  <  2:854965/150‑1 (MQ=255)
                                                                                                              tCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCt                                 >  2:167143/1‑151 (MQ=255)
                                                                                                               cAATGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCtt                                >  1:897659/1‑151 (MQ=255)
                                                                                                                aaTGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCt                                 <  2:1616850/149‑1 (MQ=255)
                                                                                                                aaTGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTc                               <  2:1411042/151‑1 (MQ=255)
                                                                                                                aaTGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTc                               <  1:511763/151‑1 (MQ=255)
                                                                                                                aaTGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGCTTGTCGCCCTTc                               <  2:915755/151‑1 (MQ=255)
                                                                                                                 aTGGTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCg                              <  1:449075/151‑1 (MQ=255)
                                                                                                                   ggTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCgg                             <  1:861396/150‑1 (MQ=255)
                                                                                                                   ggTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGt                            <  2:1168245/151‑1 (MQ=255)
                                                                                                                   ggTTGGCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGt                            >  1:200575/1‑151 (MQ=255)
                                                                                                                       ggCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACg                         <  1:924304/150‑1 (MQ=255)
                                                                                                                       ggCGCAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACg                         <  2:1482859/150‑1 (MQ=255)
                                                                                                                        gcgcAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGaa                       <  1:1453378/151‑1 (MQ=255)
                                                                                                                        gcgcAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGaa                       <  2:619650/151‑1 (MQ=255)
                                                                                                                        gcgcAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGa                        >  1:671966/1‑150 (MQ=255)
                                                                                                                         cgcAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGaa                       <  1:1509218/150‑1 (MQ=255)
                                                                                                                           cAGGCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGAc                    <  2:814553/151‑1 (MQ=255)
                                                                                                                              gCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTc                 <  1:1027093/151‑1 (MQ=255)
                                                                                                                              gCACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTc                 >  2:1218754/1‑151 (MQ=255)
                                                                                                                               cACTAAGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCt                >  2:122773/1‑151 (MQ=255)
                                                                                                                                   aaGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGc            >  2:339897/1‑151 (MQ=255)
                                                                                                                                    aGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGct           <  1:1240590/151‑1 (MQ=255)
                                                                                                                                    aGCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGct           <  1:942320/151‑1 (MQ=255)
                                                                                                                                     gCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGctc          <  1:1665868/151‑1 (MQ=255)
                                                                                                                                     gCAATCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGctc          >  2:187319/1‑151 (MQ=255)
                                                                                                                                       aaTCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGctct         <  2:81660/150‑1 (MQ=255)
                                                                                                                                        aTCGCGACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTta       <  2:1680582/151‑1 (MQ=255)
                                                                                                                                          cgcgACGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTtat      <  1:1298456/150‑1 (MQ=255)
                                                                                                                                             gacgacTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAAcc  >  1:312212/1‑151 (MQ=255)
                                                                                                                                              aCGACTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAAcc  >  2:569700/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACGACGCTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACC  >  NC_000913/3158992‑3159283

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