Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,857,513 A→C 21.2% T254P (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,5130AC21.2% 99.5 / 4.5 63T254P (ACC→CCC) gpsANAD(P)H‑dependent glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/10);  new base C (2/13);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 1.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGCACGCA  >  NC_007779/3857278‑3857750
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
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:2127498/1‑239 (MQ=255)
                                                                                                                                                                                                                                       aCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGc       >  1:1839687/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATtata                                                       >  1:142431/1‑187 (MQ=255)
                                                                                                                                                                                                                                           gTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAgg             >  1:2169251/1‑229 (MQ=255)
                                                                                                                                                                                                                                            tACCGACAACCAGTCGCGTAACCGCCGTTTTGGCTTGATGCTCGGTCACGTCATGGCTGCACAAGGAGCTTTGGAGATGATTCGTCAGGTGGTGGAAGGCCACCGCAATACGAAAGAAGTCTGTGAACTGGCGACTCGCTTGCGCGATGAAATGCCAATCACCGAGGCAACTTATGATGTATTATATTGCGGAAAAAACGCTAGCGAGGCAGCATTGATTTTCATAGGTCGTGCACGCa  >  2:1414651/1‑239 (MQ=255)
                                                                                                                                                                                                                                            tACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGCACGCa  >  2:1119967/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
AAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTTACTAGGTCGTGCACGCA  >  NC_007779/3857278‑3857750

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: