Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,857,507 A→C 19.5% T252P (ACC→CCC)  gpsA → NAD(P)H‑dependent glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,857,5070AC19.5% 90.1 / 2.4 59T252P (ACC→CCC) gpsANAD(P)H‑dependent glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (29/14);  new base C (2/10);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 2.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGC  >  NC_007779/3857270‑3857745
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
tGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTAc                                                                                                                                                                                                                                                <  1:19831/239‑1 (MQ=255)
   ggCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAGCGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCt                                                                                                                                                                                                                                              <  2:1162520/238‑1 (MQ=255)
        aaGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCg                                                                                                                                                                                                                                        <  1:1411100/239‑1 (MQ=255)
            ttCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGGTGGCGGTCTGGCGCGTCCCCTTTCCGACac                                                                                                                                                                                                                                    >  1:1229145/1‑238 (MQ=255)
             tCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGTAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCAGGGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGGGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTGTTGGTGCGGCGCGAGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAAc                                                                                                                                                                                                                                   <  2:27259/239‑1 (MQ=255)
             tCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCGTGGCGGTCTGGTGCTTACCCTTAACGACaa                                                                                                                                                                                                                                    >  1:279379/1‑238 (MQ=255)
                                       aTTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCCTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGGTCTTACCCTCACCGCCAACCAACCGCGAAACCCCCGCTTTTGCat                                                                                                                                                                                                        >  1:2137145/1‑239 (MQ=255)
                                              tGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGTTGGGGGTCTGGTGCTTACCTGTACCGGCAACCAGTCGCGTAACCGGCGGTTTGGACTGAAGCTc                                                                                                                                                                                                 >  1:664682/1‑239 (MQ=255)
                                                             cgGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGTCTGGTGCTTACCTGTTCCGACAACCAGTCGCGGAACCGCCGTTTTGGGATGGTGCTCGGTCCGGGCAttgga                                                                                                                                                                                  >  1:971405/1‑235 (MQ=255)
                                                                         ttatggcggtgggTGGGGGGCTGCCGGAGGGTATGGGTTTTGGTGGGAATGGGCGTGCGGCGTTGATACCCGGTGGGCTGGTGGACAGGTGGGTTCTGGGTGGGGCGCTGGGTGCGGCCCCTGCCCCCTTTATGGGCGGGGGGGGGCTGGGCGTTGTGGTCCTTCCCTGTACCGACACCGAGGGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAgcgcg                                                                                                                                                                      <  2:130708/228‑1 (MQ=255)
                                                                                 aTTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGGCGCGTAACCGCCGTTTTTGCATGATGCCCGGTCAGGGCATGGATGTACAAAGCGCGCAGGggaa                                                                                                                                                              >  1:1079500/1‑239 (MQ=255)
                                                                                 aTTGGTGCGGGGATGTCCGACGGTATCGGGTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGGCTTGGTGCGGCGGTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTGGGGGACGCCGGGCCTTCCTCTACCGACCCCCAGTCTCGGAACCGCCGCTTTGGCATGATGGCCGGGTGCGGACTGGATGCACAAAGCGCGCAGgagga                                                                                                                                                              >  2:701213/1‑237 (MQ=255)
                                                                                          gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAAACAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCa                                                                                                                                                     >  1:1664963/1‑239 (MQ=255)
                                                                                              tgcccGATGGTTTGGTTTTTGGTGAGAGGGGGCTGACGGCGCGGATGCCCGGGGGGCGGGCTGAAAGGGCGCGTTTTGGTGCGGCGGGGGTTGCCGACCCTGCCCCTTTTGGGGGCGTGGCGGGGCTTGGCGTTGTGGTTCTTCCCTGTACCGACACCCGGTCGCTTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAggt                                                                                                                                                  <  2:710893/235‑1 (MQ=255)
                                                                                                tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGGCTGGTGCCTACCTGTACCGACCACCCGACGCCGAACCGCCGTGTTGGGCTGTCGCTCGGGCAGGGGCTGGAGTGACAAAGAGCGCCGGGGAAGATggggcggggggg                                                                                                                                               >  2:357496/1‑230 (MQ=255)
                                                                                                tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCCTTGCGGGCTGGGGCCTTCCCCTCCCGCCAACCCATCGCGGAACCCCCCGTTTTGCCGTGCGCCCCGTGCGGGGAGGGGCGTACCAAACGCGCGGGGGAAGATTGGTCAggtggg                                                                                                                                               >  2:411035/1‑238 (MQ=255)
                                                                                                    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGCTCTGGTGCTTACCTGTACCGACAAACAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACa                                                                                                                                                                             >  1:39796/1‑205 (MQ=255)
                                                                                                          tCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGGTGGGGCTCTGGGGCCTACCTTGCCCCCCAACCCATCGCGGAACCCCCGTTTTGGGGTGAGGCTCAGGCCGGGGATGGGGTTAAAAAGCGGCCAGGGGAAAGTTGGTGCGGGGGGGGGAGGcca                                                                                                                                      >  2:1421156/1‑236 (MQ=255)
                                                                                                                  ggTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGTTGGCGGTCTGGTGCTTACCTTTAACGACACCCAGGCGCGGAACCGCCGGTTTGGGCTTGTGCCCGGTCCGGGCACGGATTGTCAAAGCGGGCCGGAGAAGATTGGGCGGGGGGGGGGAGGCTaaccccacaa                                                                                                                             >  2:1621805/1‑230 (MQ=255)
                                                                                                                  ggTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGGCTGGTGCTTACCCGTACCGGCAACCAGTCGCGTAACCCGCGTTTTTGCATGATTCTCGGTCCGGGCAAGGAGGTACAAAGCCCGCAGGGGAAGATTGGGCAGGGGGGGGAAAGGCACCGCAATAc                                                                                                                             >  2:702415/1‑239 (MQ=255)
                                                                                                                                       cTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGgc                                                                                                        >  1:459263/1‑239 (MQ=255)
                                                                                                                                              cccGTGGGCTGGCTGAAAGGTCGCGTCTGGGGGCGGCGCGGGGGGCCGACCCTGCCACCTTTAGGGGCATGGGGGGGTTTGGCGTTCTGGTGCTTACCTGTCCCGACAACCAGTCGCGTAACCGTCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGc                                                                                                 <  1:1621805/239‑1 (MQ=255)
                                                                                                                                                         gctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCTGACCCTGCCACCTTTATGGGCATGGCGGGGGTTTGGGGGCTGGGGCTTTCCCTGTCCCACCAACAGTCGCGGTACCGCCGCTTTTGGATGATGGTCGGTGCGGGCACGGAAGGACAACACGCGCCGGGGAAAAATGGGCAGGGGGGGGGAGGGTAACGCCATACGAAAGAAGACCCCCAAATCGCGCATCCCGTCgggccgga                                                                                      >  2:446638/1‑233 (MQ=255)
                                                                                                                                                           tgAAGTGTCGCGTCTGGGTGGGGCGCTGGGGGCGGACCCTGCCCCTTTTATGGGCATGGCGGGGTTTGGCGATTTGGTGCTTCCCTGTCCCGACAACCAGTCGGGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGaa                                                                                     <  1:357496/238‑1 (MQ=255)
                                                                                                                                                                tgtgGCGTCGTGGGCGGGCGCTGGGCGCCGCCCCCGCCCGTTTTGGGGGCGGGGGGGGGCTGGGCTTGTTGGTGCTCACGGCTCCCGACAACCAGCGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGcc                                                                                <  2:783496/238‑1 (MQ=255)
                                                                                                                                                                   cgcgTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTTGCGATCTGGGGCTTACCTGTACCGACAACCAGGCGCGGAACCGCCGTTTTGGCATGATGCTCGGGCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATa                                                                            >  1:1021312/1‑239 (MQ=255)
                                                                                                                                                                   cgcgTCTGGGTGGGGCGCGGGGCGCCGCCCCTGCCCCCTTTATGGGCATGGGGGGGTTTGGCGTTCTGGTGCTTCCCTGTACCGACACCCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAAt                                                                             <  1:1640028/238‑1 (MQ=255)
                                                                                                                                                                                   gtggggTCCCGCCCCTGCCCCTTTTATGCGCAGGGGGGGGTTTGCCGTTTTGTTCCTTCCTGGTCCCAAAAACCGGTGGGGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  1:609057/236‑1 (MQ=255)
                                                                                                                                                                                    tggggTCCCGCCCCTCCCACCTTTAGGGGCATGGCGGGGTTTGGGTGTTTGGTGCTTCCCGTTACCAACAACCTGTCGCGAAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTACGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  2:459263/236‑1 (MQ=255)
                                                                                                                                                                                    gtGGGTGCCGACGCCGCCTGCTTTGGGGGGGGGGGGGGGTGTGGCTATTTGGTGCTTGTCTGTCCAGACGACCGGTCGCTAAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  2:360207/237‑1 (MQ=255)
                                                                                                                                                                                    cgggggccggccccTCCCCCTTTTATGGCCGTGGGGGGCTTGGGGGGTTTGGTGTTTACCGGACCCGAGACCCGGTCGGTTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  2:1079500/227‑1 (MQ=255)
                                                                                                                                                                                    cTGGGTGCCGCCCCTGCCCCCTTTAGGGGCATGGCGGTGTTTGGGGTTTTGTTGCTTCCCTGTCCGGACACCCAGTGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  1:1881746/238‑1 (MQ=255)
                                                                                                                                                                                                 cTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGc                                              >  1:4742/1‑239 (MQ=255)
                                                                                                                                                                                                       ccTTTATGGGCAGGGCGGGGCTTGGCGTTCTGGGGTTTACCTGTCCCGACACCCAGTGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                            <  1:1397037/219‑1 (MQ=255)
                                                                                                                                                                                                          ttATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCCTGGATGTACAAAGCGCGCAGGGGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAATACGAAAGAAGGCCGCGAACAGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAc                                     >  1:1561342/1‑239 (MQ=255)
                                                                                                                                                                                                              gggCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGGAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAACGCGCAGGAGAAGATTGGTCAGGTGGGGGAAGGGTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGGTGAAATTCCAATAACCGAGGAAATTTATCAAGGATTATATTTCg                                             >  1:918902/1‑227 (MQ=255)
                                                                                                                                                                                                                gCAGGGCGGGGTTTGGCGATCGGGTGCTTCCCTGTACCGACACCCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAcgcgcg                                <  1:1804312/238‑1 (MQ=255)
                                                                                                                                                                                                                  aTGGCGGGGCTTTGCGATCTGGGGCTTACCTGTACCGACAACCAGTCGCGTAAACGCCGGTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAAAACGGAAGAAGTCCGCGAACAGGCGCATCGCTTCGGCGTTGAAATGCCAAAAACCGAGGAAAATTATCAAGACTTATTTTGTGGAAAAAACGCGCGCGa                             >  2:1487958/1‑239 (MQ=255)
                                                                                                                                                                                                                  aTGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGCCAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAgg                                                                                                                                         >  1:500107/1‑131 (MQ=255)
                                                                                                                                                                                                                  aTGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATaa                                                                           >  1:1366090/1‑193 (MQ=255)
                                                                                                                                                                                                                  aTGGCGGGGCTTGGCGATCTGGGGGTTACCTTTACCGGCAACCAGTCGCGTAACCGCCGTTTTTGGATGATGCTCGGTCAGGGCATGGGTGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGGGGAAGGGTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAAAAACCGAGGAAATTTATCCAGTATTATATTTCGGGAAAAAcgcgcccca                             >  2:1804312/1‑235 (MQ=255)
                                                                                                                                                                                                                    ggCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGGGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAgg                           >  1:145675/1‑239 (MQ=255)
                                                                                                                                                                                                                    ggCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAACGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAg                            >  1:959259/1‑238 (MQ=255)
                                                                                                                                                                                                                           gttgggcattgtGGTGTTTGCCTGTCCCGACAACCAGTGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGgcagca                      <  1:611183/227‑1 (MQ=255)
                                                                                                                                                                                                                           gtgtggggttgtGGTGCTTAGCTGGACCGACAACCGGGAGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGgcagca                      <  2:971405/227‑1 (MQ=255)
                                                                                                                                                                                                                            ttGGCTATGTGGTGCTTACCTGTACCAACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCAtt                    <  1:1313524/238‑1 (MQ=255)
                                                                                                                                                                                                                             tGGCGATCTGGTGCTTACCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa                  <  1:1674859/239‑1 (MQ=255)
                                                                                                                                                                                                                              tggggttttGGTGTTTACCTGACCCAACAACCAGTCGTGTAGCCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTCGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa                  <  1:1711999/231‑1 (MQ=255)
                                                                                                                                                                                                                               ggggttgtGGTGCTTCTCGGTGCGGAGAAGCAGGCGGGTGACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa                  <  2:1884150/231‑1 (MQ=255)
                                                                                                                                                                                                                               ggCGATCTGTTGCTTCCCTGTACGCACCACCAGCGGGGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa                  <  1:1497452/238‑1 (MQ=255)
                                                                                                                                                                                                                                 cGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACCAAACGCGCAGGAGAAGATTGGGCCGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCCTCGCTTCGGCGTTTAAATGCCAATAACCGAGGAAATTTATCAAGAATTATATTGCGGAAAAAACGCGCGCGAGGCAGCAtttcctt               >  2:1674859/1‑234 (MQ=255)
                                                                                                                                                                                                                                   aTCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGcc                                                                                >  1:1808495/1‑172 (MQ=255)
                                                                                                                                                                                                                                    tCTGGTGCTTACCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCa                                                                                                                                                                                       <  1:1333309/68‑1 (MQ=255)
                                                                                                                                                                                                                                    tCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCa                                                                                                                                                                                       >  2:1333309/1‑68 (MQ=255)
                                                                                                                                                                                                                                    tCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACttt              >  1:1775218/1‑237 (MQ=255)
                                                                                                                                                                                                                                        gTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAgg        >  1:2050101/1‑239 (MQ=255)
                                                                                                                                                                                                                                          tgtttACCTGTCCCGACACCCAGGCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAgg        <  1:411035/235‑1 (MQ=255)
                                                                                                                                                                                                                                              tACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGagaa                                                                                                                                                              >  1:891459/1‑84 (MQ=255)
                                                                                                                                                                                                                                              tACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGAAGCATTGACTTTACTAGGTCGt    >  2:199123/1‑238 (MQ=255)
                                                                                                                                                                                                                                               aCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGc  >  1:1839687/1‑239 (MQ=255)
                                                                                                                                                                                                                                               aCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGc  >  1:2127498/1‑239 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
TGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGC  >  NC_007779/3857270‑3857745

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