Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,857,502 T→G 16.9% V250G (GTG→GGG)  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,5020TG16.9% 144.4 / 5.4 76V250G (GTG→GGG) gpsANAD(P)H‑dependent glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base G (13/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCACTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTT  >  NC_007779/3857265‑3857733
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
tgcACTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTg                                                                                                                                                                                                                                            <  1:1086026/239‑1 (MQ=255)
tgcACTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTg                                                                                                                                                                                                                                            <  1:29153/239‑1 (MQ=255)
tgcACTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTg                                                                                                                                                                                                                                            <  1:286881/239‑1 (MQ=255)
tgcACTGCGGCAAAAGTTGCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCGTTGGTGCGGGGATGTCCGACGGTGTCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGGCCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTg                                                                                                                                                                                                                                            <  2:1287069/239‑1 (MQ=255)
gtcACTGCGGCAAAAGTGTCCGGGGTTACAGCAAGCCGGATTTCATTGGCGTGCAGCGTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTGGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCGCCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTg                                                                                                                                                                                                                                            <  2:1729359/237‑1 (MQ=255)
   aCTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCtt                                                                                                                                                                                                                                         <  1:881222/239‑1 (MQ=255)
   aCTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCt                                                                                                                                                                                                                                          <  1:1909236/238‑1 (MQ=255)
     tGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGTCTGGTGCTTa                                                                                                                                                                                                                                        >  1:1822324/1‑238 (MQ=255)
       cGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGTTGGCGGTCCGGGGCTTAccc                                                                                                                                                                                                                                     >  2:881222/1‑238 (MQ=255)
       cGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGTCTGGTGCTTa                                                                                                                                                                                                                                        >  1:1849395/1‑236 (MQ=255)
       cGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGTCTGGGGCttc                                                                                                                                                                                                                                        >  1:267963/1‑235 (MQ=255)
       cGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGGGCTTa                                                                                                                                                                                                                                        >  1:1640851/1‑236 (MQ=255)
        ggCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTg                                                                                                                                                                                                                                    <  1:1042785/239‑1 (MQ=255)
        ggCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTg                                                                                                                                                                                                                                    <  1:919406/239‑1 (MQ=255)
        ggCAAAAGTGTCCGAGTGTGCAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTAAGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCt                                                                                                                                                                                                                                     <  2:377529/238‑1 (MQ=255)
          cAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGGGCGGCGCTGGGTGCCGACCCTGCCCCCTTTATGGGCCTGGCAGGTGTTGGGGGTCTGGGTGCTTCCCtct                                                                                                                                                                                                                                   >  2:1732567/1‑237 (MQ=255)
                 ttCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGAc                                                                                                                                                                                                                            <  1:823610/237‑1 (MQ=255)
                 ttCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACaa                                                                                                                                                                                                                          <  1:193780/239‑1 (MQ=255)
                 ttCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACa                                                                                                                                                                                                                           <  1:1199663/238‑1 (MQ=255)
                  tCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAAc                                                                                                                                                                                                                         <  1:933876/239‑1 (MQ=255)
                                 aTCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGGGGTCTGGGGCTTACCCTTACCGACCACCAAGCGCGGAAccg                                                                                                                                                                                                           >  1:1005237/1‑238 (MQ=255)
                                  tCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCCTGGGGATCTGGTGCTTACCTGTTCCGACAACCAGGCCCGTAAccgcc                                                                                                                                                                                                         >  1:647208/1‑239 (MQ=255)
                                            aTTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCCTGGCGGTCTGGGGCTTACCTGGACCGACAAACAGACGCGTAAACGCCCGTTTGGCat                                                                                                                                                                                             >  1:553296/1‑239 (MQ=255)
                                                          tGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGGGCTTACCTGTACCGACAAACAGTCGCGTAACCGCCGGTTTGGCATGATGCTCGGTCAgg                                                                                                                                                                               >  1:135697/1‑239 (MQ=255)
                                                             cggcgCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCCGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGGCGTGGGCTTTCCCTTACCGGCCAACCGCTGGGGGACCCGCCGTTTTGGCGGGATGCCGGGCAgggga                                                                                                                                                                            >  2:558120/1‑237 (MQ=255)
                                                                                 ttGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAg                                                                                                                                                        >  1:165352/1‑239 (MQ=255)
                                                                                        tGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGGCAACCAGTCGCGTAACCGCCGGTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGagaaga                                                                                                                                                 >  1:1521833/1‑239 (MQ=255)
                                                                                                  aTGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCCCCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCCCCTTTATGGGCATGGCGGGGTTTGGCGTTCTGGTGCTTCCCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAgg                                                                                                                                       <  1:425537/238‑1 (MQ=255)
                                                                                                     tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCCTGGCGGTCTGGTGCCTTCCTCTACCGGCAACCAAGCCCGTAACCGCCGCTTTGGCCTGATTCCCGGGCAGGGGCAGGGTGGACAAAGCGCGCAAGGAAAGATTGGGTAgggggt                                                                                                                                   >  2:1898271/1‑239 (MQ=255)
                                                                                                         aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCGGGTGGTGGag                                                                                                                               >  1:506384/1‑238 (MQ=255)
                                                                                                               tCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGGGGTCTGGTGGCTACCTCTACCGACAAACCGTCGCCGAACCCCCGCTTTTGGCTGGTTGTCGGGACGGGCCTGGAGTTACAAAACGCGCCGGGGAAGAGTTGGCCGGTGGTGGaaaggca                                                                                                                          >  2:1805328/1‑233 (MQ=255)
                                                                                                                cGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGGGCTTACCTGTACCGACAACCAGTCGCGGAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAACGCGCAGGAGAAGATTGGTCAGGTGGGGGGAGGCTAcc                                                                                                                        >  1:486106/1‑239 (MQ=255)
                                                                                                                         tGCGAATGCGCGTACGGCGCTGATCCCCCGTGGGCTGGCGAAAAGGTCGCGTCTTGGTGCGGCGCGGGGTGCCGCCCCTGCCACCTTTATGGGCAGGGCGGGGCTTGGCGATCTGGTGTTTACCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAggg                                                                                                                                                                              <  1:270901/177‑1 (MQ=255)
                                                                                                                                 cgcgTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGCTGGCGGTCTGGGGCTTACCTGTACCCGCAACCAGCCGCGGAACCGCCGGTTTTGCCTGTTGCTCGGTCAGGGCATGGCTGCACAAAGCGCGCCGGAGAAGAATGGGCAGGCGGGGGAAGGCGACCGCAAAACGGaggaaag                                                                                                        >  2:1060661/1‑236 (MQ=255)
                                                                                                                                      ggcgggcTGAGCACCGGGGGGCGGGCGAACAGGTCGCGTGTGGGGGGGGCGGTGGGGCCCACCCCTGCCACCTTAGGGGGACGGGGGGGGCTGGGCGGTTGGTTGCGTCCCGGGCCCGACACCGATGGGGGAAGCCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCcgcg                                                                                                   <  2:1012683/233‑1 (MQ=255)
                                                                                                                                                 cccccGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCCCCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAgcg                                                                                                                                                             <  1:499183/168‑1 (MQ=255)
                                                                                                                                                         ggctggctGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCCGGTGGTGGAAGGCTACCg                                                                                                                       >  1:1752615/1‑199 (MQ=255)
                                                                                                                                                           ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCGTGGCGGTCTGGTGCTTTCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCGTGATGCTCGGGCAGGGCATGGATGTGCAAAACCCGCCGGggaa                                                                                                                                                   >  2:499183/1‑170 (MQ=255)
                                                                                                                                                              gctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGGTTGGCGATCTGGGGCTTACCTCTACCGACAAACCGTCCCCCAACCGCCGTTTTGGGGTGGTGGTCGGTCAGGGCATGGGAGGAACAAGCGCGCAGGGGAAGAATGGTCAGGTGGGGGGAGGCTAACGCAATACGGAAAAAGGCCGCGAACCGGCGCATCCCCTCGGCGTTgg                                                                          >  2:1721723/1‑238 (MQ=255)
                                                                                                                                                              gctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGGGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGGCAGGGGGGGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGa                                                                          >  1:323408/1‑239 (MQ=255)
                                                                                                                                                              gctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCCTGGCGAGCTGGTGCTTACCTGTTCCCGCAACCAGCCGCGTAACCGCCCTTTTTGGCTGAGTCCCGGGCCGGGCATGGAGGTACCAAGCGCGCAGGAGAAGAGTGGTCAGGGTGGGGGAGGCGCCCGCCATAAGAAAAAAATCAGCGAACTGGCGCATCGCGTCTGCgcttg                                                                          >  2:31494/1‑235 (MQ=255)
                                                                                                                                                                tgaatgggcgcgTCTGGGTGCGGCGCGGGGGCCCGCCCCTGCCACTTTTATGGGCAGGGGGGGGTTTGGCGTTCTGGTGCTTCCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGaa                                                                         <  1:1810213/230‑1 (MQ=255)
                                                                                                                                                                tgAAATGTCGCTTCTGGGGGCGGCGCGGGGGGCCGCCCCTGCCACCTTTATGGGCAGGGCGGGGCTTGGCGATTTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGaa                                                                         <  1:1898271/238‑1 (MQ=255)
                                                                                                                                                                        cgcgTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCGTGGCGATTTGGTGGTTGACTGTAACGGCAACCAAACCCGTAAACACCGGTTTTGGCTTATGCTTCGTGAGGGCATGGGTGGACAAAGAGCGCCGGGGAAGATTTGGTAGGGGGGCGGAGGACAACACCAAACGAAAAAAATCCCGGAAATGGCGCCTAGGCTTGGGGTTGAAAttgtaaaa                                                                >  2:1756094/1‑232 (MQ=255)
                                                                                                                                                                                ggggcggcgCGGGTCGCCCCCCCTCCCCCTTTTGGGGGCAGGGGGGGGCTGGGGTTTTTGGTGCTTCCTGGTACCCACACCCGGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGa                                                        <  2:135697/236‑1 (MQ=255)
                                                                                                                                                                                        gtggtgtcccGCCCCTGCCCCCTTTAGGGGCGGGGGGGGGGTGGGCGTTTTGTGCCTGACCTGACCGGACAACCAGTGGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  2:165352/231‑1 (MQ=255)
                                                                                                                                                                                        gcggggCGCCGCCCCTCCCACTTTTAGGGGAAGGGGGGGGTTTGGGTTTTTGGTGCTTCCCTGTCCCGACACCCGGGCGCGAAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGCGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  1:1539416/236‑1 (MQ=255)
                                                                                                                                                                                        gcTGGGTGCGGACCCTCCCACCTTTATGGGCATGGGGGGGCTTGGCGTTCTGTTCCTTACCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATACCCGAGGAAATTTATc                                                <  1:1307493/239‑1 (MQ=255)
                                                                                                                                                                                         gtGGGTGCCGCCCCTGCCACCTTTATGGGCATGGGGGGGCTTGGCGTTTTGGTGCTTACCTGTCCCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  1:303185/237‑1 (MQ=255)
                                                                                                                                                                                         gggggTGCCGCCCCTGCCTCCTTTATGGGCAGGGGGGGGTTTGGGGATTTGGTGCTTACCTGTCCGGACACCCAGTCGCGAAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  1:1115358/236‑1 (MQ=255)
                                                                                                                                                                                         cggggTGCCGCCCCTGCCCCCTTGAGGGGGAGGGGGGGGTTTGGCGTTCTGGTGCTTCCCGGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  1:1432156/236‑1 (MQ=255)
                                                                                                                                                                                         aggggTGCGGCCCTTGCCCCCTTTGTGGGCATGGGGGGGTTTGGCGTTTTGGTGTTTCCCTGTACCAACAACCAGTCGGGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  1:1805328/236‑1 (MQ=255)
                                                                                                                                                                                          ggggTGCCCGCCCCGCCAACTTTTTGGGGCATGGGGGTGGTGGGGGTTTTGGGGATGTCGTGTACGGAAAGCCAGTGGCGTATCCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATc                                                <  2:267963/237‑1 (MQ=255)
                                                                                                                                                                                           ggggGCCGCCCCTCCCCCTTTTGGGGGCGTGGCGGGGTTTGGCGATTTGGTGCTTCCCTGTACCGACACCCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCaa                                              <  1:1800002/238‑1 (MQ=255)
                                                                                                                                                                                           gggTGCGGCCCCTGCCACCTTTATGGGCAGGGGGGGGTTTGGCGCTTTGGTGCTTCCCTTTACCGACACCCATTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCaa                                              <  1:558120/238‑1 (MQ=255)
                                                                                                                                                                                           gggTGCCCCCCCTGCCACCTTTATGGGCAGGGGGGGGTTTGGCGTTCTGGTGCTTCCCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCaa                                              <  1:948894/238‑1 (MQ=255)
                                                                                                                                                                                                      cTGCCACCTTTATGGGCATGGCGGGGCTTGGGGGTCTGGTGCCTTCCTGCACCGACAACCAGTCGCGTAAACGCCGTTTTGGCCTGGTGCTCGGTCAGGGCATGGGTGGACAAAGCGCGCAGGAGAAGATTGGGCAGGGGGTGGAAGGGTAACGCAATACGAAAGAAGCTCGCGAACGGGGGCATCGCTTCGGCGGTGAAATGCCAATAACCGAGGAAATTTATCAAGATATATATTGc                                  >  2:643996/1‑239 (MQ=255)
                                                                                                                                                                                                      cTGCCACCTTTATGGGCATGGCGGGGCTTGGCGGTCTGGTGCTTACCTGTAACGTCCACCAGTCGCGGAACCCCCGGTTTGGCATGATGCTCGGGCAGGGCATGGATGTACAAAACGCGCAGGGGAAGATTGGTCAGGGGGGGGGAGGCTACCGCAATACGAAAGAAGTCCGCGAACCGGCGCATCGCTTCGGCGGTGAAAAGCCAATAAACGAGGGAAATTATCAAGAATTTTAttta                                  >  2:874546/1‑237 (MQ=255)
                                                                                                                                                                                                        gCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGAATTATATTGCgg                                >  1:685175/1‑239 (MQ=255)
                                                                                                                                                                                                        gCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGGCAGGTGGGGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCgg                                >  1:247936/1‑239 (MQ=255)
                                                                                                                                                                                                               ttATGGGCATGGCGGGGCTTGGCGATCTGGGGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGGGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAc                         >  1:1571671/1‑239 (MQ=255)
                                                                                                                                                                                                                  tGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAcgc                       >  1:1141653/1‑238 (MQ=255)
                                                                                                                                                                                                                    gggCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAcgcgc                     >  1:1701597/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCAGGGGGGGGCTTGGTGATTTGGTGTTTCCCTGTCCCGCAAACCAGTGGCATAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAAcgcgcg                    <  1:839277/238‑1 (MQ=255)
                                                                                                                                                                                                                       cATGGCGGGGCTTGGCGGTCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTGAGGGCATGGATGTACAAAACGCGCAGGGGAAGATTGGTCAGGGGGGGGAAGGCTACCGCACAACGAAAGAAGACCGCGAACTGGGGCATCGCTTTGGCGGTGAAATGCCCATAACCGAGGAAATTTATCAAGTCTTATATTTTGGAAAAAAagggcgcg                  >  2:408803/1‑239 (MQ=255)
                                                                                                                                                                                                                         tGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATaa                                                               >  1:24763/1‑192 (MQ=255)
                                                                                                                                                                                                                         tGGCGGGGCTTGGCGATCTGGGGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGGGGTGGAAGGCTACCGCAATACGAAAGAAGGCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAATATTATATTGCGGAAAAAACGCGCGCGAg                >  1:1636579/1‑239 (MQ=255)
                                                                                                                                                                                                                            cGGGGCCTGGCGATCTGGTGCGTACCTGTACCGACAACCAGGCGCGGAACCGCCGTTTTGGCATGATGCCCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGGCAGGTGGGGGAAGGCGAACGCAATAAGAAAGAAAGCCGCGAAATGGCGCATCGCCTCGGCCTTGAAAAGCCAATAAACGAGGAAAATTAACAAATATTATATTTCGGAAAAAAACCGCGAGAggga             >  2:1398191/1‑237 (MQ=255)
                                                                                                                                                                                                                               ggtttGGCGTTGTTGTGCTTACCTGTACCGACAACCAGGCTCTTAACCGCCGTTTTGACATGATGCTCGGTCAGTGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCACGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGgcagca          <  2:50126/236‑1 (MQ=255)
                                                                                                                                                                                                                                   tGGCGATCTGGTGCTTACCTGTACCGACAACCAGGCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa      <  1:281073/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gggggATGTGTTGCTTCCCTGTACCAACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa      <  1:1756094/235‑1 (MQ=255)
                                                                                                                                                                                                                                   gggCGGTCTGGTGCTTCCCTGTCCCGCCAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGa      <  1:874546/238‑1 (MQ=255)
                                                                                                                                                                                                                                    tggcaattctgTGCTTACCTCTACGTCCAACCGTTCTCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTAGTCAGGTGGTGGAATGCTACCGCAATACGAAAGAAGTCCGCTAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCTAGGCAGCATTGa      <  1:1097074/229‑1 (MQ=255)
                                                                                                                                                                                                                                     gCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACt    >  1:592350/1‑239 (MQ=255)
                                                                                                                                                                                                                                         tCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACttt  >  1:473930/1‑237 (MQ=255)
                                                                                                                                                                                                                                          cTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGACCTTTTGGCATTATGCTCGGTCAGGCCATGGATGTACAAAGCGCGCAGGAGAAGTTGGGTCAGGGGGGGGAAGGCTACCGCAATACGAAGGAATCCCGCGTACTGGCGCATCGCATCGACGTTGAAATGCCGATCACAGAGGAAATTTATCCAGTATtata                                      >  2:862302/1‑200 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
TGCACTGCGGCAAAAGTTTCCGCGTTTACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGGGCTTGGCGATCTGGTGCTTACCTGTACCGACAACCAGTCGCGTAACCGCCGTTTTGGCATGATGCTCGGTCAGGGCATGGATGTACAAAGCGCGCAGGAGAAGATTGGTCAGGTGGTGGAAGGCTACCGCAATACGAAAGAAGTCCGCGAACTGGCGCATCGCTTCGGCGTTGAAATGCCAATAACCGAGGAAATTTATCAAGTATTATATTGCGGAAAAAACGCGCGCGAGGCAGCATTGACTTT  >  NC_007779/3857265‑3857733

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: