Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,106 A→G N139N (AAT→AAC)  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,630,1060AG86.4% 97.3 / 3.3 46N139N (AAT→AAC) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (16/22);  total (20/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

GTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCT  >  NC_000913/4629863‑4630352
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
gTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGAc                                                                                                                                                                                                                                                 <  1:188476/251‑1 (MQ=255)
  aaCGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACcg                                                                                                                                                                                                                                               >  1:384866/1‑251 (MQ=255)
                accgcccccacgttgacctCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGGTGGTCGGTGCGTCGGTCAGCTGCATGTCTGTTTTTGCCAGCAGCAGGCGGCAAAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCTTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTAAGTTCAGGTTATGACCGTCGTGAGCCTGAAt                                                                                                                                                                                                                                 <  2:45584/233‑1 (MQ=255)
                        caCGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCtt                                                                                                                                                                                                                         <  1:179145/250‑1 (MQ=255)
                        caCGGTGCCCTCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCCCAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCtt                                                                                                                                                                                                                         <  1:168481/250‑1 (MQ=255)
                                cctCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACgg                                                                                                                                                                                                                 <  1:110781/249‑1 (MQ=255)
                                    gAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCtt                                                                                                                                                                                                             >  1:368110/1‑251 (MQ=255)
                                      aGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGCTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGt                                                                                                                                                                                                           >  1:9463/1‑251 (MQ=255)
                                          gTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGAATCCAGGCGGTTGGTCGGTTCGTCGAGCAGCAGAATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGCCGACGCTCACCACCAGAGAGGATAGCGATTTTAGCGTCACAATCCGGCAGAAGTAGCGCATCCGCCGCACGCTCCAGCTGAACGTTCAGCAGATGACCGACGTGAGCCTGACAGACCTCTTCCAAACCGCCTTGTTc                                                                                                                                                                                                         >  2:370394/1‑249 (MQ=255)
                                                   gCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCt                                                                                                                                                                                              <  1:349988/251‑1 (MQ=255)
                                                          agccagGCCACGGATTCGGCATCCCGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTCTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCCGCCTGTCaa                                                                                                                                                                                        >  2:58972/1‑250 (MQ=255)
                                                           gccagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCaaa                                                                                                                                                                                       >  1:229797/1‑250 (MQ=255)
                                                             cagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATc                                                                                                                                                                                    >  1:191384/1‑251 (MQ=255)
                                                                              aTCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGCCGTGAGCCTGCATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAACTCGGCATCCGGCTCGGCGt                                                                                                                                                                   >  2:21168/1‑251 (MQ=255)
                                                                                    gtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAgcg                                                                                                                                                             <  1:275728/251‑1 (MQ=255)
                                                                                          ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATaca                                                                                                                                                       <  1:93909/251‑1 (MQ=255)
                                                                                          ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATaca                                                                                                                                                       <  1:83874/251‑1 (MQ=255)
                                                                                          ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATac                                                                                                                                                        >  1:110205/1‑250 (MQ=255)
                                                                                                     cgAGCAGCCGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACACCGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGCCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAgg                                                                                                                                            <  2:9463/251‑1 (MQ=255)
                                                                                                     cgAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAgg                                                                                                                                            <  1:237776/251‑1 (MQ=255)
                                                                                                      gcgctgcagcaTGTCTGGTTTTTCCAGCAGCAGGCGGCTCAACGTTACGCTCCGACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCATGCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCTTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                                                                                                           <  2:110205/246‑1 (MQ=255)
                                                                                                      gAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                                                                                                           <  1:377042/251‑1 (MQ=255)
                                                                                                       ggaggcagccTGTGTGGTTTGGCCAGCAGCGGGGGGGAGAACGGTACGCGACGAGGCTCGCCACCGGAGAGGGTAGGGATGTTCGCGTCCCAGTCGGGCAGACGTAGGGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAGCTCTGCCAGACGGGCTTGTTCAGCGGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCg                                                                                                                                          <  2:111651/247‑1 (MQ=255)
                                                                                                              gcaatgcTGGTTTTTCCAGCAGCAGGCGGCACAACCCTACGCGACGGCGCTCACCACCGGAAAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCACCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCACGTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATTCACTTCATCCAGGCGTTTCAgg                                                                                                                                   <  1:249593/245‑1 (MQ=255)
                                                                                                                      ggTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCAGCCGCCGCACGCTCCAGCTGTACGTTCAGATTTTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAAc                                                                                                                           <  1:55967/251‑1 (MQ=255)
                                                                                                                              cagcagcagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAgc                                                                                                                                                              <  1:260804/208‑1 (MQ=255)
                                                                                                                                    ccggcggcAAAACGCTACGAGACGACGCTCACCACCGGAGAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCCGGTTATGACCGTCGAGAGACTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGCAGGACAGCGCATACACTTCATCCAGGCGTTTCAGGACGTTAACCACGTCAGAAAc                                                                                                               >  2:170772/3‑249 (MQ=255)
                                                                                                                                         ggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGcttct                                                                                                        <  1:21168/251‑1 (MQ=255)
                                                                                                                                                   tACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGTCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACt                                                                                              >  1:267630/1‑251 (MQ=255)
                                                                                                                                                       caacgcctctccccaacGGAGAGTTTAGCTATTTCTGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTGCAGGGTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTcacg                                                                                          <  2:369613/235‑1 (MQ=255)
                                                                                                                                                        gacgacGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTATGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTcacgc                                                                                         >  1:90104/1‑251 (MQ=255)
                                                                                                                                                        gacgacGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCATTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTcacgc                                                                                         >  1:25917/1‑251 (MQ=255)
                                                                                                                                                                             ggTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGc                                                                    >  1:185306/1‑251 (MQ=255)
                                                                                                                                                                                tAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCAAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCCGCGGCCAGCTTGTCACAATCGGCATCCGGATCGGCGTCCAGCGCATACACTTCATCCCGGAGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTAAGCTgg                                                                 >  2:122815/1‑250 (MQ=255)
                                                                                                                                                                                        ttCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTATGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTctt                                                          <  2:90104/250‑3 (MQ=255)
                                                                                                                                                                                        gtctcttcccatcccGGCAGACGATGCGCATCCGCCGCCCGCTCCAGTTGTACGTTCAGCTTATGACCGTCGTGAGCCTGAATGCTCTCTACCATACGGCCTTTTTCAGCGGCCAGTTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACGCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTTAACCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCTGGTTCAGCCGCGGTTctt                                                          <  2:25917/237‑3 (MQ=255)
                                                                                                                                                                                             gTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCa                                                    <  1:203501/251‑1 (MQ=255)
                                                                                                                                                                                                gctcttccgatctcGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggt                                                 >  1:214976/14‑249 (MQ=255)
                                                                                                                                                                                                 cAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggta                                                <  1:122815/251‑4 (MQ=255)
                                                                                                                                                                                                 cAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCTTGCGGCAGATa                                                <  1:58972/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGa                                     <  1:171118/250‑1 (MQ=255)
                                                                                                                                                                                                            aCGTAGCGCATCCGCCGCACGCTCCAGCTGAACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGACGTACAGCGCATACACTTCATCCACGCGTTTCAGGGCGTTAACCACTTCAGATACCGTTTCTTCCATCGACTCAAGCACGGTGTGTTCCGGGTTCAGCTGCGGATCCTCGGGCAGGTAACCAATCtta                                      >  2:250532/1‑249 (MQ=255)
                                                                                                                                                                                                                           ctcACGCTCCAGTTGTACGTCCAGGTTATGACCTGAGTGTGTCTGAATGATCTCTGCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACATCGCAACCACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGTAACCGCTTCTTCAATGGa                                                                                                <  2:335333/175‑1 (MQ=255)
                                                                                                                                                                                                                           cgcACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGa                                                                                                >  1:335333/1‑177 (MQ=255)
                                                                                                                                                                                                                                       cTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCtt          >  1:51801/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCt  >  1:197682/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCT  >  NC_000913/4629863‑4630352

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