Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,897 C→T A542A (GCG→GCA ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,8970CT98.6% 218.9 / NA 73A542A (GCG→GCAettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (34/37);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCA  >  NC_000913/4628655‑4629140
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
aCGGTTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGAAGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATTCGCTTCGGCTCCAGTGCGTCTgc                                                                                                                                                                                                                                             <  2:399047/251‑1 (MQ=255)
   cgatccccGCTACAAAATTGCGGAGACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTTCGTCTGCGcc                                                                                                                                                                                                                                          <  1:187064/246‑1 (MQ=255)
         ccGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGcc                                                                                                                                                                                                                                          >  2:187064/1‑245 (MQ=255)
          cgctaCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtg                                                                                                                                                                                                                                   >  1:580402/5‑251 (MQ=255)
           gctaCGAACTTGCGGATACTAAGGCTGTTGTAGGCCGGCAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGCTAAGGCGTTCACGCCGCCTCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTAGGCTCCAGTTCTTCTGCGCCCAgcgtgc                                                                                                                                                                                                                                  >  2:232791/4‑251 (MQ=255)
             taCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAATTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtgcgt                                                                                                                                                                                                                                >  2:48748/2‑251 (MQ=255)
              aCGAAATGGCGGATTCTAAGGCTGTTGGAGGCCGAAAAGGGCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCTACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAAGACGCTTGTACTTGATACGCTTCGGCGCCAGTGCGTCTGCGCCCAgtgtgcg                                                                                                                                                                                                                                 <  2:307218/249‑1 (MQ=255)
                      gcggATACTAAGGCGGTGTTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGGTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGAGAAGGCTTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGTTACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTAc                                                                                                                                                                                                                       <  1:70430/249‑1 (MQ=255)
                          atacTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGAATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTCCGCATTACTTCGCAATCCGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCtt                                                                                                                                                                                                                   >  2:444608/4‑251 (MQ=255)
                              tAAGGCTGTTGTAGGCCGGAAACGGCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGCGATTCACAACGTATTGAAATTTCTTGCTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGTTACGCCTCCGCTCACGATCGTCTGAGCCCAGCGTGCAGTTCTTGTACTCTTCGt                                                                                                                                                                                                                >  2:547331/1‑250 (MQ=255)
                                    tgttgtAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTcag                                                                                                                                                                                                          <  1:431006/250‑3 (MQ=255)
                                     tttgagggccGGAAAAGGCGTTCGCGCCGTATCCGGCATTCGGTGCACAATGCCTGATGGGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTcag                                                                                                                                                                                                          <  1:586396/241‑3 (MQ=255)
                                        gtAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGttcgcagtaaa                                                                                                                                                                                                      <  2:52530/250‑11 (MQ=255)
                                        gtAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTaa                                                                                                                                                                                                       >  2:210169/1‑249 (MQ=255)
                                                   aaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCgaa                                                                                                                                                                                          <  2:163842/246‑1 (MQ=255)
                                                   aaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAATTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCgaa                                                                                                                                                                                          <  1:48748/246‑1 (MQ=255)
                                                        gtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACt                                                                                                                                                                                     <  2:8248/251‑1 (MQ=255)
                                                           cACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGCCATTTCTTGGTTTGTAGGCCGGATAAGACGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCa                                                                                                                                                                                   >  1:54991/1‑250 (MQ=255)
                                                            ccgccgtATCCGGCATGCGTTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATGGAAATTTCTTGGTTTGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACGTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAAc                                                                                                                                                                                 <  1:555744/250‑1 (MQ=255)
                                                            aCGCCGTATCCGGCATTCGCTGCTCAATGCATGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAAc                                                                                                                                                                                 >  1:488235/1‑251 (MQ=255)
                                                             cgccgtATCCGGCATTCGGAGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCTCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACt                                                                                                                                                                                >  2:371586/1‑251 (MQ=255)
                                                              gccgtATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACtt                                                                                                                                                                               >  1:105234/1‑251 (MQ=255)
                                                                   atccggctgctggTGCTCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCGTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATTCGCGTGTACTTGATACGCTTCGGCTCCAGTGCGGCTGCTCCCAGCGTGCGTTGCTTTTACTCTTCGTACTCAGTAAATTTACCGTCGAATAACGCATCTTTAcc                                                                                                                                                                           <  2:72499/239‑1 (MQ=255)
                                                                                 gCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTagtc                                                                                                                                                            <  1:163207/251‑4 (MQ=255)
                                                                                      aTGCCTGATGCGACGCTTACAGCGTCTTCTCAGGCCTGGAGATTCACACCGTATTGAAATTTCTTGGCTTCTAGGCCGGATAAGGCGTTCACGCCGCATACGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTACAGTTACCTTCGAAGCACTCAACTTTACCTTCCTCCTGGTAGTCCAgg                                                                                                                                                        >  2:248854/1‑249 (MQ=255)
                                                                                      aTGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTGGCGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGCAGAACTCAACTTTACCTTCATCCTGGTAGTCCAgg                                                                                                                                                        >  2:68144/1‑249 (MQ=255)
                                                                                          cTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGAtgtg                                                                                                                                                   >  2:395388/1‑251 (MQ=255)
                                                                                            gATGCGACGCTTACCGCGTCTGAGCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGc                                                                                                                                                  <  1:59369/250‑1 (MQ=255)
                                                                                               gCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCGTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTgg                                                                                                                                              <  1:264838/251‑1 (MQ=255)
                                                                                                 gACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTagt                                                                                                                                                             >  1:401956/1‑232 (MQ=255)
                                                                                                 gACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCCGGTagt                                                                                                                                                             <  2:401956/234‑3 (MQ=255)
                                                                                                   cGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTGCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt                                                                                                                                          <  1:357324/251‑1 (MQ=255)
                                                                                                   cGCTTACCGCGTCTTATCAGGCCGGGAGATTCCCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCGGCGCCCAGCGTGCGGTTCTTGTACTCTTCGTACTCAGTAAAGTGACCGTCGAAGAACTCAACTTTACCTTCATCCGGGTAGTCCAGGATGTGCGTGGCGAt                                                                                                                                          <  2:580402/251‑1 (MQ=255)
                                                                                                            cgTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGCACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCg                                                                                                                                   >  1:57901/1‑249 (MQ=255)
                                                                                                             gTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAgg                                                                                                                                >  2:17750/1‑251 (MQ=255)
                                                                                                              tCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAgg                                                                                                                                >  2:513705/1‑250 (MQ=255)
                                                                                                              tCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAgg                                                                                                                                >  1:484227/1‑250 (MQ=255)
                                                                                                                     aggcctggagattcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACg                                                                                                                        >  2:298647/15‑251 (MQ=255)
                                                                                                                     aggcctggagattaaCAACGTATTGAAATTTCTTGATTTGAAAGCCGGCTAAAGAGTTCACGACGCATACGGAATTTTACGCATTACTTCGCAATACGATTGTACTTGATACGCTTCGACTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGACCTCAGTAAAGATACCTTCGAAGAACTCCACTTTACCTTAATCCTGGGCGTCCAGGATTGGCGTGGCGATACGATCGAgcgaccaccg                                                                                                                        >  2:350501/15‑242 (MQ=255)
                                                                                                                               attcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCgag                                                                                                              <  1:442066/247‑1 (MQ=255)
                                                                                                                                              gaaatttcttggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGcaca                                                                                                <  1:169732/239‑1 (MQ=255)
                                                                                                                                                aatttcttggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGc                                                                                              >  1:38127/11‑251 (MQ=255)
                                                                                                                                                 atttcttggtttGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGcc                                                                                             <  2:105234/242‑1 (MQ=255)
                                                                                                                                                    tcttggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTAGGCTCCAGTGCGTATGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCgg                                                                                          >  1:438546/7‑251 (MQ=255)
                                                                                                                                                       tggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGATCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGa                                                                                        >  1:559128/4‑250 (MQ=255)
                                                                                                                                                           ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTc                                                                                    <  2:219658/250‑1 (MQ=255)
                                                                                                                                                                   cGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAggcc                                                                           <  2:477268/251‑3 (MQ=255)
                                                                                                                                                                    ggATAAGGCGTTCCCGCCGCCTCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGACCCGCGTGCGTTTCTTGTAATCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAAATTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGCACCAACGGTCGTGCGAGATAACCATAGCACAGCCCGGGAACTCCAGCAGGGCg                                                                          >  2:371932/1‑251 (MQ=255)
                                                                                                                                                                      aTAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGtt                                                                        >  2:501408/1‑251 (MQ=255)
                                                                                                                                                                          ggCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTcc                                                                    >  2:393528/1‑251 (MQ=255)
                                                                                                                                                                          ggCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTcc                                                                    <  1:210169/251‑1 (MQ=255)
                                                                                                                                                                                cacgccgcATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAgcgcg                                                              <  2:238178/249‑1 (MQ=255)
                                                                                                                                                                                  cgccgcATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAgcgcgcg                                                            >  2:62907/1‑251 (MQ=255)
                                                                                                                                                                                   gccgcATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCTCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAgcgcgcgc                                                           <  1:371586/251‑1 (MQ=255)
                                                                                                                                                                                         tCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGCTACGCTTCGGCTCCCGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGtt                                                     >  2:520315/1‑251 (MQ=255)
                                                                                                                                                                                                 tttACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGCTTCGATATc                                             >  2:353722/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGATTCGATATcc                                            >  2:336607/1‑249 (MQ=255)
                                                                                                                                                                                                      gCATTGCTTCGCTATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTGCTTGGTCTCTTCGTACTCAGGAAAGGTACCTTCGAATAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGGGGCGATACGGTCGAGGATCCAACGGTCGTGCGAGATAACCATCGCATAGCTCGGAAACTGCAGCAGGGCGTTTTCCAGCGCGCGCATGGTGTCGATAGCCAGGt                                        <  2:581536/251‑1 (MQ=255)
                                                                                                                                                                                                      gCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGa                                                                                                                                  <  1:263181/161‑1 (MQ=255)
                                                                                                                                                                                                      gCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGa                                                                                                                                  >  2:263181/1‑161 (MQ=255)
                                                                                                                                                                                                           aCTTCTCAATACGCTTGTACTTGATACGCTGCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCgttg                                   <  2:12297/251‑1 (MQ=255)
                                                                                                                                                                                                           aCTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCgttg                                   <  2:12283/251‑1 (MQ=255)
                                                                                                                                                                                                                 cAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTCCCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTtggttggt                             <  1:185724/251‑1 (MQ=255)
                                                                                                                                                                                                                     tcgcgtGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTtcg                          <  2:359738/245‑1 (MQ=255)
                                                                                                                                                                                                                     aCGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGAGCAGGATTTCGATATCCAGGGCGTTGGTTGGTtcgt                         >  2:401007/1‑251 (MQ=255)
                                                                                                                                                                                                                      cGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCGTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTtcgtc                        <  1:371932/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTtcgtcg                       >  2:31945/1‑249 (MQ=255)
                                                                                                                                                                                                                                 gATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGagcagcagc              <  2:228279/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aTACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGGGCAGCAGCAt            >  1:230386/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGATTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGtt         <  1:401007/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTGTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGGTTTCCAGCGCGCGCAGGGTTTCGATATCCAGTTCGTTTGTTTGTTCGTCGAGCAGCAACATGtt         <  2:392602/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ttCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTgcc      >  2:557084/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCa  <  1:31945/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCA  >  NC_000913/4628655‑4629140

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