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

AACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTC  >  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)
                     ccgcTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAgcgcgc                                                                                                                                                                                                                           >  2:396224/1‑251 (MQ=255)
                          tGCTGAATTAACGTCCTGCTTTCAGCTTCTGAGAATACTCTTCATAAATGCGGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCTTTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCctt                                                                                                                                                                                                                        <  1:486228/249‑1 (MQ=255)
                           gCTGAATTAACGTCCCGCTTTCAGCTTCTGATAATACCCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCcttctt                                                                                                                                                                                                                     <  2:217351/251‑1 (MQ=255)
                                aTTAAAGTCCTGCTTTCAGCTTCTGCTAATACTCTTCATAACTGCTGCTGGCTGCGCCAACGTCATTCTGCAATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGCGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGCTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTGCTTTGttt                                                                                                                                                                                                                >  2:69918/1‑251 (MQ=255)
                                        ccTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGTTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCAt                                                                                                                                                                                                         <  1:23042/250‑1 (MQ=255)
                                                  gCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATc                                                                                                                                                                                              <  1:72387/251‑1 (MQ=255)
                                                         ataataCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGc                                                                                                                                                                                       >  1:170046/1‑251 (MQ=255)
                                                               ctctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTCACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGCTCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGtcc                                                                                                                                                                                  >  2:85071/1‑250 (MQ=255)
                                                               ctctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGtccc                                                                                                                                                                                 >  2:131076/1‑250 (MQ=255)
                                                                               ctggctgCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATgccgcc                                                                                                                                                                 >  1:148032/1‑251 (MQ=255)
                                                                                       cccAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGCCTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCg                                                                                                                                                         <  2:374082/250‑1 (MQ=255)
                                                                                             gTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTCGAGTGTTTTATCGTTAGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGCTAGTTTCAGCAACCTGTTTTGACACATCCGGGCGCAACAGGAAGTTGATCAATTTAAGCACGCCTTATTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGTAAATACCGCATTCTTTCGGCcacc                                                                                                                                                   >  2:11766/1‑250 (MQ=255)
                                                                                                         ttCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGACGCAAGGTTTGGCGTTGGATAACCGATAGATTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTg                                                                                                                                       >  1:65743/1‑251 (MQ=255)
                                                                                                          tCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGACCACACCACGTCAATTgg                                                                                                                                      >  1:288765/1‑251 (MQ=255)
                                                                                                                   ttAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCg                                                                                                                              <  2:170046/250‑1 (MQ=255)
                                                                                                                                               ttttATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                       >  2:190830/1‑245 (MQ=255)
                                                                                                                                               ttttATCGTTCGCCACTTCTGGACTTAACAGCGTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                       <  1:190830/245‑1 (MQ=255)
                                                                                                                                                                  tGGACTTAACAGCTTACGCGCCGCAAGGATTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCACTTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGTACATGCCGAGGTTa                                                                               >  2:48826/1‑250 (MQ=255)
                                                                                                                                                                    gACTTAACAGCTTACGCGCCGCACGGATTGGCGTTGGGTAACCAATCGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTAACCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTCACt                                                                            >  1:174481/1‑251 (MQ=255)
                                                                                                                                                                           cagctaccgcgCCGCAAGGTGTGGCGTTGGATACCCGATCGTTTCAGCAACCTGTTTTGCCCCATCCGGGCGCACCAGGAGGTTCATCAATGTCAGCGCCCCTTCGTTGTTTCTGTCAGTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCTGCTTTCTTTCGTCCACACCCCATCAATTGTCGTACCCGCCTGTCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACCCCGCCtt                                                                     <  2:95749/244‑1 (MQ=255)
                                                                                                                                                                                tACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTAGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATg                                                                >  1:249845/1‑251 (MQ=255)
                                                                                                                                                                                     gccgcAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCACTTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCACTTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACgg                                                           >  2:310609/1‑251 (MQ=255)
                                                                                                                                                                                         aaaGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATTCGGGCGCACCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTc                                                       <  1:276418/250‑1 (MQ=255)
                                                                                                                                                                                          aaGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCg                                                      >  1:412253/1‑251 (MQ=255)
                                                                                                                                                                                                tGGCTTTGGAAAACCGATATTTGCAGCACTCTGTTTTGCCACATCCGGGCGCAACAGGAGGTTGACCAATTTCAGCGCGCCTCCTTTGTTTTGGGCATTTGCCGGGATCGCCAGGCTGTCCAGCCAGAAAATGCCGCCTTCTTACGGCCACACCACGTCAATGGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCggg                                                 <  2:418492/250‑1 (MQ=255)
                                                                                                                                                                                                    gtttgtaaaCCGTTAGTTTCAGCAACCTGTTTTGCCCCACCCGGGCGCTACAGGAAGTTGATCAATTTCACGCCGCCTTCTTTGTTTTTGTCATTTGCCGGGATCGCCAGCCTGTCCATCCAGAAAATGCCGCCTGCTTTCGGCCACACCCCTTCAATTGCCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTAt                                             <  2:291416/243‑1 (MQ=255)
                                                                                                                                                                                                       ggATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGa                                         <  1:109279/251‑1 (MQ=255)
                                                                                                                                                                                                          tAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGtt                                      >  1:95588/1‑251 (MQ=255)
                                                                                                                                                                                                             ccGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTaa                                    >  2:409852/1‑250 (MQ=255)
                                                                                                                                                                                                                          cAACCTGTTTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAGGGCTGCCACGttt                      >  2:107144/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cacaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCa                                                                                                          >  2:100558/1‑154 (MQ=255)
                                                                                                                                                                                                                                       cacaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCa                                                                                                          <  1:100558/154‑1 (MQ=255)
                                                                                                                                                                                                                                        ccaTCCTGGCGCAACAGTAAGTTGATCAATTTCAGCGCGCCTTCTTTTTTTTTGGCATTTGCCGGGAGCGCCAGCCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCCCGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGTCTGCCACGTTTGGCATCAGtttttt        <  1:423391/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ccaTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGtttttt        <  1:131076/250‑1 (MQ=255)
                                                                                                                                                                                                                                         caTCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTAGCATTTGCGGGGATCGCCAGGCTGTCCATCGAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATGGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGATTGAAAGGCTGCCACGTTTGGCATCAGTTTTTTc       <  2:249845/251‑1 (MQ=255)
                                                                                                                                                                                                                                         caTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCCGTTTTTTc       >  1:41937/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTCCGGGCGCCACAGGAAGTTGATAAATTTAAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTCCGGGTTCGCCGGGTTATCGGAGTTAACGGCTGCCACGTTTGGCATCAGTTTTTTc       >  2:227247/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gggCGCATCAGGAGGTTGATTCATTTCAGCGCGCCTTCTTTGTTTTTGTCATTTGCCGGGATCGCCAGTCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACCCCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTc  <  1:420541/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
AACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTTTTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTC  >  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: