Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,449 A→G 5.7% intergenic (+57/‑68) yejH → / → rplY putative ATP‑dependent DNA or RNA helicase/50S ribosomal subunit protein L25

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,4490AG 5.7% 315.2 / 6.7 105intergenic (+57/‑68)yejH/rplYputative ATP‑dependent DNA or RNA helicase/50S ribosomal subunit protein L25
Reads supporting (aligned to +/- strand):  ref base A (47/51);  new base G (4/2);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 105 total aligned reads displayed.

CGCCCGCATGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAGTTCCCG  >  NC_000913/2282302‑2282597
                                                                                                                                                   |                                                                                                                                                     
gtcccTCATGAAAGGCCTGTACTGGCAGGTGCGTGATAAAGTGTTCGATTACTAAGGTCGTTTTCTTCTGGAGCACTAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATa                                                                                                                                                     <  1:330248/148‑1 (MQ=255)
   ccGCATGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGc                                                                                                                                                  <  2:805793/150‑1 (MQ=255)
    cGCATGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑a                                                                                                                                                <  1:683337/150‑1 (MQ=255)
    cGCATGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑a                                                                                                                                                <  1:718161/150‑1 (MQ=255)
    cGCATGAAAGGCCAGTACTGGAAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑a                                                                                                                                                >  1:1201283/1‑150 (MQ=255)
        tGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGc                                                                                                                                            <  1:648844/150‑1 (MQ=255)
        tGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGc                                                                                                                                            >  2:818105/1‑150 (MQ=255)
          aaaGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAg                                                                                                                                          >  2:679645/1‑150 (MQ=255)
           aaGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGa                                                                                                                                         <  2:1484148/150‑1 (MQ=255)
            aGGCCAGTACTGGCAGGTGCGTGAAAAACTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCGTATGC‑aaaagcggn                                                                                                                                        <  2:56020/150‑7 (MQ=255)
             ggCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATc                                                                                                                                       <  1:400562/150‑1 (MQ=255)
             ggCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATc                                                                                                                                       >  1:1277715/1‑150 (MQ=255)
             ggCCAGTACTGGCAGGGGCGAGAAAAAGTGTTCGATTACGAAGGGCGGTTTCGTCGGGCGCACGAATTTCGCGGTTAATGGTAGTTTTGATTGAGGTATAAGGCGATTGAGTAAAGAATCTCGCCCGCTTTTGCGTATGC‑aaaagcgggc                                                                                                                                       <  1:1563460/150‑8 (MQ=255)
              gCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCa                                                                                                                                      <  2:1578885/150‑1 (MQ=255)
                cAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACt                                                                                                                                    >  2:558592/1‑150 (MQ=255)
                  gTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTa                                                                                                                                  >  2:886619/1‑150 (MQ=255)
                    aCTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGAATGATGTATAAGGCGATTGAGTTTAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTAcc                                                                                                                                <  2:1703325/150‑1 (MQ=255)
                     cTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCACTTTTGCATACGC‑AAAGCAGATCACTTACCt                                                                                                                               >  1:1141823/1‑150 (MQ=255)
                       ggCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGt                                                                                                                             >  1:1539475/1‑150 (MQ=255)
                        gCAGGTGCGTGATAAAGTGTTCGTTTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTATTGGTAGTTTTGATTTATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGtt                                                                                                                            <  1:1170047/150‑1 (MQ=255)
                         cAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTg                                                                                                                           <  1:1210942/150‑1 (MQ=255)
                         cAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTg                                                                                                                           >  2:773966/1‑150 (MQ=255)
                            gTGCGTGAAAAAGTGTTCGATTACGTAGGTCGTTTTCGTCTGGCGCACGATATACGCGGTTAATGGTAGTTTTGATTGATGTATTAGGCGATTGAGTATAGTATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTg                                                                                                                        <  1:301042/150‑1 (MQ=255)
                               cGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCACGCTTTTGCATACGC‑AAAGCAGATCACTTAACTGTTGCTGGGt                                                                                                                     >  1:1119658/1‑150 (MQ=255)
                                   aaaaaGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGcc                                                                                                                 <  1:1232559/150‑1 (MQ=255)
                                    aaaaGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCt                                                                                                                >  1:772011/1‑150 (MQ=255)
                                       aGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGGTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTa                                                                                                             <  1:1315807/150‑1 (MQ=255)
                                           ttCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAg                                                                                                         <  1:1389543/150‑1 (MQ=255)
                                            tCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAgg                                                                                                        <  2:844569/150‑1 (MQ=255)
                                             ggATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGTTTGATGTATAAGGCGATTGAGTATAGTATCTCGCCCGCTTTTGCATACGC‑AATGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGa                                                                                                       <  1:10771/149‑1 (MQ=255)
                                             cGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGa                                                                                                       >  2:1693751/1‑150 (MQ=255)
                                              gATTACGAAGGTCGTTTTCGTCTGGCGCTCGAATTACTCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTAGTGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGAt                                                                                                      <  1:113634/150‑1 (MQ=255)
                                                ttACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTa                                                                                                    >  1:396460/1‑150 (MQ=255)
                                                 tACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTaa                                                                                                   >  1:10937/1‑150 (MQ=255)
                                                   cGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAAtt                                                                                                 <  1:230604/150‑1 (MQ=255)
                                                    gAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAAttt                                                                                                <  2:493441/150‑1 (MQ=255)
                                                     aaGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTa                                                                                               >  1:453188/1‑150 (MQ=255)
                                                     aaGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTa                                                                                               >  2:1525241/1‑150 (MQ=255)
                                                       ggTCTTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAag                                                                                             <  1:1432730/150‑1 (MQ=255)
                                                       ggTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAGTTTAag                                                                                             <  2:757188/150‑1 (MQ=255)
                                                         tagTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAagag                                                                                           <  2:1165494/148‑1 (MQ=255)
                                                         tCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAagag                                                                                           <  2:1262126/150‑1 (MQ=255)
                                                          cGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAagaga                                                                                          >  2:1468395/1‑150 (MQ=255)
                                                               tCGTCTGGCGCACGAATTACGCTGTTAATGGTAGTTTTGAATGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGagaaag                                                                                     <  2:149808/150‑1 (MQ=255)
                                                                 gTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGagaaagaa                                                                                   <  2:466710/150‑1 (MQ=255)
                                                                  tCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAgaaagaaa                                                                                  <  1:557637/150‑1 (MQ=255)
                                                                      gcgcACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGtt                                                                              >  1:1156335/1‑150 (MQ=255)
                                                                        gcACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTa                                                                            >  2:826605/1‑150 (MQ=255)
                                                                         cACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGAGATTGATTATAGAATCTCGCCCGCTTTTGCATACGC‑AAATCAGATCTCTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAATAGAGAAAGAAATGTTTAc                                                                           <  1:911882/150‑1 (MQ=255)
                                                                            gAATTACGCGGTTTATGGTAGTTTTGATTGATGTATAAGGCGATTGATTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTAt                                                                        <  1:1341431/150‑1 (MQ=255)
                                                                             aaTTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATc                                                                       >  1:908253/1‑150 (MQ=255)
                                                                               ttACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGGCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCac                                                                     >  1:775264/1‑149 (MQ=255)
                                                                                 aCGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACg                                                                   <  2:1670305/150‑1 (MQ=255)
                                                                                  cgcgGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGc                                                                  <  2:1373117/150‑1 (MQ=255)
                                                                                     ggTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGa                                                               <  1:460913/150‑1 (MQ=255)
                                                                                      gTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGaa                                                              >  2:750014/1‑150 (MQ=255)
                                                                                       ttAATGGTAGTCTTGATTGATTTATAATGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTATCTGTTTCTGGGTCGCCTGTATCAGGATTAATTTAAGATAGAAAGAAATGTTTACTATCAACGCAGAAg                                                             <  2:1101935/150‑1 (MQ=255)
                                                                                         aaTGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTa                                                           >  1:660012/1‑150 (MQ=255)
                                                                                          aTGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTAc                                                          <  2:1665398/150‑1 (MQ=255)
                                                                                           tGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACg                                                         >  2:598102/1‑150 (MQ=255)
                                                                                              tAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTaa                                                      >  2:1061503/1‑150 (MQ=255)
                                                                                                 ttttGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAaga                                                   >  1:1102270/1‑150 (MQ=255)
                                                                                                   ttgattgaTGTATAAGGCGATTGAGTATAGAAGCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGc                                                 <  2:1285124/150‑1 (MQ=255)
                                                                                                      tttgatGTATAAGGCGATTGAGTATATAATCTCGCCAGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAgg                                              <  2:291239/149‑1 (MQ=255)
                                                                                                       ttgatGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCTGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGTAAGAAATGTTTACTTTCATCGCTGAAGTACGTTAAGAGCAggg                                             >  1:194285/1‑150 (MQ=255)
                                                                                                        tgatGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGt                                            >  2:1780899/1‑150 (MQ=255)
                                                                                                         gatGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGAAGAAGTACGTAAAGAGCAGGGTa                                           >  1:1163555/1‑150 (MQ=255)
                                                                                                           tGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAg                                         >  1:926940/1‑150 (MQ=255)
                                                                                                            gTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACATGTTGCTGGGTCGCCTGTAGAAGGATTAATTTAAGATAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAgg                                        >  1:1596135/1‑150 (MQ=255)
                                                                                                             tataAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAggg                                       <  2:660012/150‑1 (MQ=255)
                                                                                                             tataAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCAATTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAggg                                       <  2:117062/150‑1 (MQ=255)
                                                                                                               taAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTTCCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTg                                     <  2:1632578/150‑1 (MQ=255)
                                                                                                               taAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTg                                     >  1:381452/1‑150 (MQ=255)
                                                                                                                 aGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCg                                   >  1:773328/1‑150 (MQ=255)
                                                                                                                   gCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAg                                 <  1:444338/150‑1 (MQ=255)
                                                                                                                    cGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCCCTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgc                                >  1:1118051/1‑150 (MQ=255)
                                                                                                                    cGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgc                                <  1:728047/150‑1 (MQ=255)
                                                                                                                      atggAGTTTTGAATCTCGCCCGGTTTTGGATTCGC‑AAAGGAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCTGGATTAATTTAAGAGAGAAAGAATTGTTTACTATCATCGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgccg                              <  2:1165527/147‑1 (MQ=255)
                                                                                                                      aTTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgccg                              <  1:886619/150‑1 (MQ=255)
                                                                                                                       ttGAGTATAGAATCTCGCCCGCTTTTGCCTACGC‑AAAGCAGCTCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGAGTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgccgc                             >  1:300788/1‑150 (MQ=255)
                                                                                                                         gAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATATAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAgccgccg                           >  1:186685/1‑150 (MQ=255)
                                                                                                                             ataGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTg                       <  2:782476/150‑1 (MQ=255)
                                                                                                                             ataGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGAAGGGTAAGGGTGCGAGACGCCGCCTg                       >  2:312139/1‑150 (MQ=255)
                                                                                                                                gAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGTGTAAGGGTGCGAGCCGCCGCCTGCGt                    <  1:773966/150‑1 (MQ=255)
                                                                                                                                 aaTCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTg                   <  1:1693751/150‑1 (MQ=255)
                                                                                                                                     tcGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGc               <  2:772011/150‑1 (MQ=255)
                                                                                                                                       gCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGTTGCGAGCCGCCGCCTGCCTGCCGCTa             >  2:228509/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTa             >  1:127878/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTa             >  2:788108/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGCTCGCCTGTAGCAGCATTAATATAAGAGAGAAAGAAATGTAAACAATCAACGCAGAAGTACGTAAAGAGAAGGGAAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTa             >  2:223099/1‑150 (MQ=255)
                                                                                                                                        cccGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTaa            >  1:1646268/1‑150 (MQ=255)
                                                                                                                                        atctgctttGCGTATGCAAAAGCAGATCAGTTACCTGTTGCTGGGTCGCCTGTAGCTGGATTAATTTAAGAGAGAAAGAAATGGTTACTATCAACGCAGAAGTACGTAAAGCGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTa             >  1:80940/7‑150 (MQ=255)
                                                                                                                                        atctgctttGCGTATGCAAAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTa             >  1:1244587/7‑150 (MQ=255)
                                                                                                                                          cGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACa          <  2:1418809/150‑1 (MQ=255)
                                                                                                                                            gctttGCGTATGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGACAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAg        >  2:106080/3‑150 (MQ=255)
                                                                                                                                              tttGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAGtt      >  1:648574/1‑150 (MQ=255)
                                                                                                                                                tGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTACGTGCCACTAACAAGTTcc    <  2:1539475/150‑1 (MQ=255)
                                                                                                                                                  cATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGCAGTACGTAAAGCGCAGGGTAAGGGTCCGCGCCGCCGCCTGCGGCCCGCTACCAAGTTCCCg  >  2:1444122/1‑150 (MQ=255)
                                                                                                                                                  cATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAGTTCCCg  >  1:933712/1‑150 (MQ=255)
                                                                                                                                                  cATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAGTTCCCg  >  2:1581696/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCCCGCATGAAAGGCCAGTACTGGCAGGTGCGTGAAAAAGTGTTCGATTACGAAGGTCGTTTTCGTCTGGCGCACGAATTACGCGGTTAATGGTAGTTTTGATTGATGTATAAGGCGATTGAGTATAGAATCTCGCCCGCTTTTGCATACGC‑AAAGCAGATCACTTACCTGTTGCTGGGTCGCCTGTAGCAGGATTAATTTAAGAGAGAAAGAAATGTTTACTATCAACGCAGAAGTACGTAAAGAGCAGGGTAAGGGTGCGAGCCGCCGCCTGCGTGCCGCTAACAAGTTCCCG  >  NC_000913/2282302‑2282597

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