Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,146,384 A→G V71A (GTA→GCA)  yceF ← m(7)GTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,3840AG95.3% 129.3 / ‑2.3 43V71A (GTA→GCA) yceFm(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (19/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

ACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGAT  >  NC_000913/1146137‑1146628
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCt                                                                                                                                                                                                                                                   <  1:308938/251‑1 (MQ=255)
     aaaaGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCt                                                                                                                                                                                                                                              <  2:181312/250‑1 (MQ=255)
                    aGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTaa                                                                                                                                                                                                                                >  2:133880/1‑250 (MQ=255)
                     gTGCAGGGGATGCCCTTTACGCACGTAAGTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGATCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTGTAGAAGGTGACGATATTGCGGCTGGCTTTGCGTAATTGCGGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCCATAATTAAAt                                                                                                                                                                                                                              <  1:56989/251‑1 (MQ=255)
                         aGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATc                                                                                                                                                                                                                          <  1:310892/251‑1 (MQ=255)
                         aGGGGATGCTATTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAGCAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGCCGGGCATTTTCTTCCGTTAACGGTTTTCCGGTGATTTCGCCATCAAGACCACATACCTGGTCTGAACCAATAATTAAATGATc                                                                                                                                                                                                                          >  2:407719/1‑251 (MQ=255)
                          ggggATGCTCCTTACGCACGTAATTACCAATCTCCGCCTCGCTCAGGTGACGGAAAGTGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACTCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATcc                                                                                                                                                                                                                         <  2:464345/251‑1 (MQ=255)
                           gggATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTAAGGCGACGGAAATGGCCGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAACCAACGACAGTCCAGTATAGAAGGTGACGCTATTGCCGCTGGCTTTGCGTAATTGCAGCCGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGCACACCTGCCTGGTCTGAACTCATCATTAAAttctccg                                                                                                                                                                                                                        >  2:60391/1‑245 (MQ=255)
                           gggATGCTCTTTACGCACGTAATTATCAATCTCCGCCACGCTCAGGCGAAGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACCTGCCTGGTCTGAACCAATAATTAAATGATCCg                                                                                                                                                                                                                        >  1:409501/1‑251 (MQ=255)
                                   ctTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGt                                                                                                                                                                                                                >  1:386994/1‑251 (MQ=255)
                                     ttACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTg                                                                                                                                                                                                               >  2:182060/1‑250 (MQ=255)
                                         gCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGcc                                                                                                                                                                                                          >  2:179935/1‑251 (MQ=255)
                                         gCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGcc                                                                                                                                                                                                          <  1:428983/251‑1 (MQ=255)
                                           aCGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCag                                                                                                                                                                                                        <  1:115439/251‑1 (MQ=255)
                                                aTTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAAGTGACGCCAGAGAtt                                                                                                                                                                                                   >  1:85796/1‑251 (MQ=255)
                                                   aTCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTgc                                                                                                                                                                                                 <  1:182060/250‑1 (MQ=255)
                                                                       ggTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCCGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTATTGTGCCAGTCg                                                                                                                                                                             >  1:196981/1‑250 (MQ=255)
                                                                                   ggACGTCAAAAGGCTCCACTGCTGTGTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCAGCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGt                                                                                                                                                                <  2:114470/251‑1 (MQ=255)
                                                                                           aaaGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTg                                                                                                                                                         >  2:70459/1‑250 (MQ=255)
                                                                                            aaGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGAt                                                                                                                                                       <  1:179935/251‑1 (MQ=255)
                                                                                             aGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGAtt                                                                                                                                                      >  1:178712/1‑251 (MQ=255)
                                                                                                cgccACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATGAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCa                                                                                                                                                                     <  1:174100/231‑1 (MQ=255)
                                                                                                cTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGACGCACCa                                                                                                                                                                     >  2:174100/1‑233 (MQ=255)
                                                                                                            ttGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTggg                                                                                                                                        <  2:126955/250‑1 (MQ=255)
                                                                                                                                     aacaacGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATaa                                                                                                                                                                                                                   <  2:248135/150‑1 (MQ=255)
                                                                                                                                     aacaacGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATaa                                                                                                                                                                                                                   >  1:248135/1‑150 (MQ=255)
                                                                                                                                       caacGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACAtt                                                                                                            <  1:288850/251‑1 (MQ=255)
                                                                                                                                                  cagtATAGAAGGTGCAGATATTGCCGCTGGCTTTGCGTAATTGCAGCCGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCGAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGCTAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCCACTGTAGCGGTGATTCGTCGCTGCGTGGGGATTCGTCGACCTCTGGTGCTGCACCTTCAAAAGAGAtt                                                                                                 >  2:189305/1‑251 (MQ=255)
                                                                                                                                                       taGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTtc                                                                                                                                   <  1:125614/212‑1 (MQ=255)
                                                                                                                                                       taGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCACTGCGTGGGGTTtc                                                                                                                                   >  2:125614/1‑212 (MQ=255)
                                                                                                                                                           aggcgACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGtttt                                                                                        >  1:225469/5‑251 (MQ=255)
                                                                                                                                                                aCGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGCTAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGCTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTcc                                                                                     >  1:34350/1‑249 (MQ=255)
                                                                                                                                                                 cGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTg                                                                                                                                                                 <  2:168897/172‑1 (MQ=255)
                                                                                                                                                                 cGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTg                                                                                                                                                                 >  1:168897/1‑172 (MQ=255)
                                                                                                                                                                                cctttcGTATTTGCAGCCGGGCATTATCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCGGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAt                                                                   <  2:413717/246‑1 (MQ=255)
                                                                                                                                                                                  ttGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAt                                                                   <  1:70459/249‑1 (MQ=255)
                                                                                                                                                                                              aGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGc                                                     >  2:373889/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATAttt                               <  1:264735/251‑1 (MQ=255)
                                                                                                                                                                                                                           accggttgttttgCCATCAAGAACACATGCCTGGTCTTAACCAATAATTAAAGGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTTGGGTTTCGTCGTCCTCTGGTGCTGCACATTCAAAATAGGTTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCa                        <  2:295238/239‑1 (MQ=255)
                                                                                                                                                                                                                           aCCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACTCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCa                        <  1:407719/251‑1 (MQ=255)
                                                                                                                                                                                                                               gTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCGCCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGTAGGCTAAAATAAGTTTAGTCATATTTTTTTCCATCAg                    <  1:133880/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aTCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTAt         <  1:255470/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 acacATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTgatgat  >  1:337049/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGAT  >  NC_000913/1146137‑1146628

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