Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,181,995 G→A L279L (CTG→TTG)  potD ← spermidine preferential ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,181,9950GA92.9% 119.1 / ‑6.1 45L279L (CTG→TTG) potDspermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (25/14);  minor base C (1/2);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTC  >  NC_000913/1181750‑1182238
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
aaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAgg                                                                                                                                                                                                                                                  >  1:394152/1‑250 (MQ=255)
    ggCGGTAATACCAACGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTATTCAAAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGGTAGTTTCAGAAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGtt                                                                                                                                                                                                                                             >  1:399969/1‑251 (MQ=255)
       ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGAt                                                                                                                                                                                                                                          >  1:415221/1‑251 (MQ=255)
        gTAATACCAACGCCCGCTTGCTGAAATAACGTCCTGCATTCAGCTTCAGATCATACTCTTCATAAATGCTGCTGGCTACGCCAACGTCATTATGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGGAGAGTGTTTTATAGTTCGACACATCTGGCCTTACCAGCTTACGCGCCGCAAGGTTGGGAGTTGGATAACCGATAGTTTAAGCAACCTGTTTTGCAAAATCCGGGCGAACCAGGAGGATGATc                                                                                                                                                                                                                                         >  2:36037/1‑251 (MQ=255)
            tACCACCGCCCGCTTGCTGAATTAACGTCCTGATTTCAGCTTCTGCTAATACTCTTCATACCTGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTCATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCACGGTTTGGCGTTGGATACCCGTTAATTTCAGCAACCTTTTTTGCAACATCCGGGCGCAACTGGCAGTTGATCAAtt                                                                                                                                                                                                                                     >  2:43279/1‑251 (MQ=255)
              ccaccGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTc                                                                                                                                                                                                                                   >  2:109279/1‑251 (MQ=255)
                  cgcccgcTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAgcg                                                                                                                                                                                                                               <  1:107908/251‑1 (MQ=255)
                     ccgcTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAgcgc‑gc                                                                                                                                                                                                                           >  2:396224/1‑251 (MQ=255)
                          tGCTGAATTAACGTCCTGCTTTCAGCTTCTGAGAATACTCTTCATAAATGCGGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCTTTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCctt                                                                                                                                                                                                                        <  1:486228/249‑1 (MQ=255)
                           gCTGAATTAACGTCCCGCTTTCAGCTTCTGATAATACCCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCcttctt                                                                                                                                                                                                                     <  2:217351/251‑1 (MQ=255)
                                aTTAAAGTCCTGCTTTCAGCTTCTGCTAATACTCTTCATAACTGCTGCTGGCTGCGCCAACGTCATTCTGCAATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGCGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGCTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTGCTTTGttt                                                                                                                                                                                                                >  2:69918/1‑251 (MQ=255)
                                        ccTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGTTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCAt                                                                                                                                                                                                         <  1:23042/250‑1 (MQ=255)
                                                  gCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATc                                                                                                                                                                                              <  1:72387/251‑1 (MQ=255)
                                                         ataataCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGc                                                                                                                                                                                       >  1:170046/1‑251 (MQ=255)
                                                               ctctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTCACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGCTCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGtcc                                                                                                                                                                                  >  2:85071/1‑250 (MQ=255)
                                                               ctctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGtccc                                                                                                                                                                                 >  2:131076/1‑250 (MQ=255)
                                                                               ctggctgCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATgccgcc                                                                                                                                                                 >  1:148032/1‑251 (MQ=255)
                                                                                       cccAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGCCTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCg                                                                                                                                                         <  2:374082/250‑1 (MQ=255)
                                                                                             gTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTCGAGTGTTTTATCGTTAGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGCTAGTTTCAGCAACCTGTTTTGACACATCCGGGCGCAACAGGAAGTTGATCAATTTAAGCAC‑GCCTTATTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGTAAATACCGCATTCTTTCGGCcacc                                                                                                                                                   >  2:11766/1‑250 (MQ=255)
                                                                                                         ttCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGACGCAAGGTTTGGCGTTGGATAACCGATAGATTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTg                                                                                                                                       >  1:65743/1‑251 (MQ=255)
                                                                                                          tCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGACCACACCACGTCAATTgg                                                                                                                                      >  1:288765/1‑251 (MQ=255)
                                                                                                                   ttAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCg                                                                                                                              <  2:170046/250‑1 (MQ=255)
                                                                                                                                               ttttATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                       >  2:190830/1‑245 (MQ=255)
                                                                                                                                               ttttATCGTTCGCCACTTCTGGACTTAACAGCGTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                       <  1:190830/245‑1 (MQ=255)
                                                                                                                                                                  tGGACTTAACAGCTTACGCGCCGCAAGGATTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCACTTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGTACATGCCGAGGTTa                                                                               >  2:48826/1‑250 (MQ=255)
                                                                                                                                                                    gACTTAACAGCTTACGCGCCGCACGGATTGGCGTTGGGTAACCAATCGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTAACCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTCACt                                                                            >  1:174481/1‑251 (MQ=255)
                                                                                                                                                                           cagctaccgcgCCGCAAGGTGTGGCGTTGGATACCCGATCGTTTCAGCAACCTGTTTTGCCCCATCCGGGCGCACCAGGAGGTTCATCAATGTCAGCGC‑CCCTTCGTTGTTTCTGTCAGTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCTGCTTTCTTTCGTCCACACCCCATCAATTGTCGTACCCGCCTGTCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACCCCGCCtt                                                                     <  2:95749/244‑1 (MQ=255)
                                                                                                                                                                                tACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTAGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATg                                                                >  1:249845/1‑251 (MQ=255)
                                                                                                                                                                                     gccgcAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCACTTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCACTTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACgg                                                           >  2:310609/1‑251 (MQ=255)
                                                                                                                                                                                         aaaGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATTCGGGCGCACCAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTc                                                       <  1:276418/250‑1 (MQ=255)
                                                                                                                                                                                          aaGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCg                                                      >  1:412253/1‑251 (MQ=255)
                                                                                                                                                                                                tGGCTTTGGAAAACCGATATTTGCAGCACTCTGTTTTGCCACATCCGGGCGCAACAGGAGGTTGACCAATTTCAGCGC‑GCCTCCTTTGTTTTGGGCATTTGCCGGGATCGCCAGGCTGTCCAGCCAGAAAATGCCGCCTTCTTACGGCCACACCACGTCAATGGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCggg                                                 <  2:418492/250‑1 (MQ=255)
                                                                                                                                                                                                    gtttgtaaaCCGTTAGTTTCAGCAACCTGTTTTGCCCCACCCGGGCGCTACAGGAAGTTGATCAATTTCA‑CGCCGCCTTCTTTGTTTTTGTCATTTGCCGGGATCGCCAGCCTGTCCATCCAGAAAATGCCGCCTGCTTTCGGCCACACCCCTTCAATTGCCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTAt                                             <  2:291416/243‑1 (MQ=255)
                                                                                                                                                                                                       ggATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGa                                         <  1:109279/251‑1 (MQ=255)
                                                                                                                                                                                                          tAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGtt                                      >  1:95588/1‑251 (MQ=255)
                                                                                                                                                                                                             ccGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTaa                                    >  2:409852/1‑250 (MQ=255)
                                                                                                                                                                                                                          cAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAGGGCTGCCACGttt                      >  2:107144/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cacaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCa                                                                                                          >  2:100558/1‑154 (MQ=255)
                                                                                                                                                                                                                                       cacaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCa                                                                                                          <  1:100558/154‑1 (MQ=255)
                                                                                                                                                                                                                                        ccaTCCTGGCGCAACAGTAAGTTGATCAATTTCAGCGC‑GCCTTCTTTTTTTTTGGCATTTGCCGGGAGCGCCAGCCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCCCGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGTCTGCCACGTTTGGCATCAGtttttt        <  1:423391/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ccaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGtttttt        <  1:131076/250‑1 (MQ=255)
                                                                                                                                                                                                                                         caTCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTAGCATTTGCGGGGATCGCCAGGCTGTCCATCGAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATGGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGATTGAAAGGCTGCCACGTTTGGCATCAGTTTTTTc       <  2:249845/251‑1 (MQ=255)
                                                                                                                                                                                                                                         caTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCCGTTTTTTc       >  1:41937/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTCCGGGCGCCACAGGAAGTTGATAAATTTAAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTCCGGGTTCGCCGGGTTATCGGAGTTAACGGCTGCCACGTTTGGCATCAGTTTTTTc       >  2:227247/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gggCGCATCAGGAGGTTGATTCATTTCAGCGC‑GCCTTCTTTGTTTTTGTCATTTGCCGGGATCGCCAGTCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACCCCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTc  <  1:420541/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
AACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGC‑GCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTC  >  NC_000913/1181750‑1182238

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: ‑