Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 831,431 C→T S187F (TCC→TTC)  rhlE → ATP‑dependent RNA helicase RhlE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913831,4310CT98.1% 173.6 / ‑6.3 52S187F (TCC→TTC) rhlEATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/25);  total (26/26)
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.93e-01

CGTTCGCTGGTGGTGTTTGGTGGTGTGAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTCCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCATCTGGCTGAACAGC  >  NC_000913/831184‑831664
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
cgttcgCTGGTGGTGTTTGGTGGTGTGAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGGGGCAACCCCGGGACGCTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGCTATCCGTCGTGCGTTAACAAAACTACCTGCGAAGCGCCAGAACATGTTATTCTTcgc                                                                                                                                                                                                                                        >  2:146605/1‑251 (MQ=255)
           ggtgTTTGGTGGTGTGAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTctct                                                                                                                                                                                                                              <  1:376064/250‑1 (MQ=255)
            gtttttGGTGGTTTGAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGTCGTTTGCTGGATCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAACTCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCTCAAGAACCTGTTCGTCTCCGCGACCTTCTCTGa                                                                                                                                                                                                                            <  2:461800/248‑1 (MQ=255)
                       tgtgAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTaa                                                                                                                                                                                                                   >  2:121245/1‑249 (MQ=255)
                        gtgAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAg                                                                                                                                                                                                                 <  2:151673/250‑1 (MQ=255)
                          gAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGAGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGccc                                                                                                                                                                                                              >  2:375622/1‑251 (MQ=255)
                           aGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCCACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCt                                                                                                                                                                                                             >  2:276711/1‑251 (MQ=255)
                              aTTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCGCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGc                                                                                                                                                                                                          <  2:507381/251‑1 (MQ=255)
                                tAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCg                                                                                                                                                                                                         <  1:216875/250‑1 (MQ=255)
                                        aGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTg                                                                                                                                                                                                >  2:438456/1‑251 (MQ=255)
                                        aGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTg                                                                                                                                                                                                >  1:10010/1‑251 (MQ=255)
                                                aaaCTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGcaca                                                                                                                                                                                         >  2:222636/1‑250 (MQ=255)
                                                 aaCTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATAGCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAAc                                                                                                                                                                                       <  2:542304/251‑1 (MQ=255)
                                                                 tGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGa                                                                                                                                                                         >  1:333196/1‑249 (MQ=255)
                                                                  gATGTGCTGGTGGCAACCCCGTGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAg                                                                                                                                                                       <  2:643672/250‑1 (MQ=255)
                                                                   aTGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTg                                                                                                                                                                     >  1:217577/1‑251 (MQ=255)
                                                                     gtgCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGTCATGGGCTTTATCCACGATATCCGTCGCGTTTTAACAAACCTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGTCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGc                                                                                                                                                                   <  1:40217/251‑1 (MQ=255)
                                                                        cTGGTGGCGACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCCCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGACATCGAAGTGGNAcc                                                                                                                                                                >  2:590872/1‑250 (MQ=255)
                                                                            tggCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCAcgccgc                                                                                                                                                            >  1:490744/1‑251 (MQ=255)
                                                                                                  ggACCTGAAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTTTTAACAAAACTGCCTGCGAAGCGCCATAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGAAAAAACTGTTGCACAACCCGCTGTAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGt                                                                                                                                        <  2:515446/249‑1 (MQ=255)
                                                                                                   gACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGAc                                                                                                                                     >  2:192138/1‑251 (MQ=255)
                                                                                                     catgatacagcAGAATGTAGTGACGGTGGATCAGGTTGACATCCTCGATCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCGCGATATCCGTCGCGTGTTAACAAAACTCCCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCTGAAAAACTGTTGCTTAACCCGCTGGAAATCGAAGTGGCACGCCGCAATTCCGCGTCTGATCAGGTGAccc                                                                                                                                   <  2:104143/241‑1 (MQ=255)
                                                                                                        ggAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAg                                                                                                                                 <  1:252432/250‑1 (MQ=255)
                                                                                                           aCATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACg                                                                                                                             >  1:181618/1‑251 (MQ=255)
                                                                                                             aTCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGtt                                                                                                                           <  1:411640/251‑1 (MQ=255)
                                                                                                               cAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTc                                                                                                                          >  2:75303/1‑250 (MQ=255)
                                                                                                                     gCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTg                                                                                                                    <  1:139100/250‑1 (MQ=255)
                                                                                                                     gCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGt                                                                                                                   <  2:10010/251‑1 (MQ=255)
                                                                                                                      cAGTGAAGCTGGATCAGGTTGAAATCCTCGGCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGCTATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTc                                                                                                                  <  2:508800/251‑1 (MQ=255)
                                                                                                                            aGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCACTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAg                                                                                                            <  1:653004/251‑1 (MQ=255)
                                                                                                                              cTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCCGCACGTTCACTTTGTCGATAAGa                                                                                                           >  2:353113/1‑250 (MQ=255)
                                                                                                                              cTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGa                                                                                                           >  2:469998/1‑250 (MQ=255)
                                                                                                                                ggATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGaa                                                                                                          <  1:276711/249‑1 (MQ=255)
                                                                                                                                          gAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAAcgcg                                                                                               <  1:325155/250‑1 (MQ=255)
                                                                                                                                             aTCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAAtt                                                                                           <  1:587344/251‑1 (MQ=255)
                                                                                                                                                        cgATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAACCGCGAATTGCTGTCGcac                                                                                 >  2:392274/1‑250 (MQ=255)
                                                                                                                                                        cgATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGcac                                                                                 <  1:222636/250‑1 (MQ=255)
                                                                                                                                                                  gACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGa                                                                       >  1:157216/1‑250 (MQ=255)
                                                                                                                                                                          gCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGAGAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAAc                                                               <  2:549287/250‑1 (MQ=255)
                                                                                                                                                                             cGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTg                                                             <  1:375622/249‑1 (MQ=255)
                                                                                                                                                                              gACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGgca                                                          <  2:333196/251‑1 (MQ=255)
                                                                                                                                                                                  tGGGCTTTCTCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGgcagc                                                        >  1:136234/1‑249 (MQ=255)
                                                                                                                                                                                                gATATACGTCACGTGTTAACAAAACTACCTGCGAAGAGCCCGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACCTGATTGGGACAGGGAACTGGCAGCAGGTGATGGTGTTTAc                                        >  2:293089/1‑251 (MQ=255)
                                                                                                                                                                                                           cctgtTAACAAAACTACCTGCGAAGCGCCAGATCCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGATCCAGCACGTTCTCTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTCCCCGTACCAAac                             <  2:637795/249‑1 (MQ=255)
                                                                                                                                                                                                                tAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGgc                        >  2:392166/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaaCTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTaa                   <  2:227065/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCTACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTCCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGAGCTGGTGTTTACCCGTACCAAACACGGCGCTAAcc                 >  2:88946/1‑251 (MQ=255)
                                                                                                                                                                                                                         tACCTGCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCAt               >  1:399590/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCCTGCGAAGGGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCAt               <  1:167875/250‑1 (MQ=255)
                                                                                                                                                                                                                              gCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTCAAGCCCTGGCGGAAAAACTGTTGCACAACACGCTGGAAATCGAAGTGGCACGCCGCAATACCGAGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGAGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCATCTGGc          >  1:615315/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCGAAGCGCCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCATCTGGc          >  1:67094/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gcCAGAACCTGTTATTCTTCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCATCTGGCTGAACAGc  >  1:520362/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
CGTTCGCTGGTGGTGTTTGGTGGTGTGAGTATTAACCCGCAGATGATGAAACTGCGTGGCGGCGTTGATGTGCTGGTGGCAACCCCGGGACGTTTGCTGGACCTGGAACATCAGAATGCAGTGAAGCTGGATCAGGTTGAAATCCTCGTCCTCGATGAAGCTGACCGCATGCTCGACATGGGCTTTATCCACGATATCCGTCGCGTGTTAACAAAACTACCTGCGAAGCGCCAGAACCTGTTATTCTCCGCGACCTTCTCTGACGATATTAAAGCCCTGGCGGAAAAACTGTTGCACAACCCGCTGGAAATCGAAGTGGCACGCCGCAATACCGCGTCTGATCAGGTGACCCAGCACGTTCACTTTGTCGATAAGAAACGCAAACGCGAATTGCTGTCGCACATGATTGGGAAAGGGAACTGGCAGCAGGTGCTGGTGTTTACCCGTACCAAACACGGCGCTAACCATCTGGCTGAACAGC  >  NC_000913/831184‑831664

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