Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,845 A→C 19.2% V197G (GTG→GGG)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,8450AC19.2% 120.5 / 3.6 62V197G (GTG→GGG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (28/17);  new base C (3/12);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 6.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967611‑2968080
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
cctgTTTTTCCCCCCCACCCTGCCCCCGGGCTTTAAACACCTCCGAAAAAAATTCCCCGACGGAGCGGGGTTTCCGAAGGCGGGTGGGAGAAAAGACCAAAAAGAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATacac                                                                                                                                                                                                                                              <  2:98722/233‑1 (MQ=255)
aCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTGGCGTTCCAGGGCCGGATCCAGCCTTGATGCCTGTTAcccccg                                                                                                                                                                                                                                           >  1:1293154/1‑237 (MQ=255)
            gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTTGTTTGCCGCTTCTTCCACGGCCCCGGTCAATCGTTCGCGGTCCTAGGCCGGCTCCACGCTTGTGCCCTGATAAAACCGTCCCATTaaaag                                                                                                                                                                                                                              >  2:1281508/1‑237 (MQ=255)
            gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCCTTTTCCGCTTCTTCCACCGCCCCGGGCAGTCGTTCGCGTTCCAGCGCCGGATCCCACGTGGATGCCCGATACAACTGGTCCATTAAAgg                                                                                                                                                                                                                              >  2:1062448/1‑239 (MQ=255)
                            tctttttATCAACCCGCGAAAAAAATTCCGGGCGGCGCCGGGTGTCCGAAAGGAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGTTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGc                                                                                                                                                                                                              <  2:1293154/235‑1 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTTCCGCTTCTTCCAGCGCCCCGGGCGGTCTCTCCCGTTCCCGACCCGGCGCCACGTGTGAAGCCTGGCACCTCTGCTCCCTTAAAAGGAAAGGGGcctcctt                                                                                                                                                                                                               >  2:549934/1‑237 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCCTTAAAGGGAACGTGGCATCCTGc                                                                                                                                                                                                              >  1:294342/1‑239 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTg                                                                                                                                                                                                               >  1:436198/1‑238 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAACGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATAAACCTGTTCCATTAAAGGTAACGTGGCACCCTg                                                                                                                                                                                                               >  1:56631/1‑238 (MQ=255)
                                   aTCAACCTGGGCAAACAATCCCGGACCAAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCtt                                                                                                                                                                                                       <  1:683396/239‑1 (MQ=255)
                                         cTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCa                                                                                                                                                                                                  >  1:464027/1‑238 (MQ=255)
                                           cGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTGCCGCTTCTTCCATCGCCCCGGTCAAGCGGTCACGTTCCAGAGCCGGGTCCAGCCTTTATGCCTGATAAACCTGGTCCATTAAAGACACCGAGGCAACCGGCCACCCTTCGACCATc                                                                                                                                                                                               >  2:406947/1‑239 (MQ=255)
                                            ggCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAACCTCTTTTTCCGCTTCTTCCCGCGCCCCGGTCAGGCGTTCGCGTTCCAGATCCGGATCCACCGTTGATGCCTGATACCCCCTTTCCATTCAAGGTACCGCGCCACCCTCCCAccccctggccacg                                                                                                                                                                                               >  2:85209/1‑228 (MQ=255)
                                            ggCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCCGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  1:53215/1‑238 (MQ=255)
                                            ggCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  1:217390/1‑238 (MQ=255)
                                                                        gtCCGAAAGCGGGGGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCCCCGGCGGCAAAGCGTTTGCTGTAGCTGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                                                                                                                    <  1:510801/238‑1 (MQ=255)
                                                                                  gggggCGAGTAAGAATCGAAAATGGCCGCCGTTTCTTTTTGGGCCCCGGCGGCAAAGCCTTTGCTGGAGCGGCGAAACTGGTTTGCCGGTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          <  2:536349/234‑1 (MQ=255)
                                                                                   ggTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCCGTTCCCTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGCCAATCCCCACACCCGGTGCTCCCGGTAATTGCCGGATTctcg                                                                                                                                                          >  1:407001/1‑237 (MQ=255)
                                                                                    gTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          >  1:177079/1‑236 (MQ=255)
                                                                                           tAAAGTCCGAAAATCGCCCCCGTTTTTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                                                                                                                                 <  2:294342/238‑1 (MQ=255)
                                                                                                 ccGAAAATAGCCGCCGTTTCTTGTTGGGCACCGGCGGCAAAGCGTTGGCTGTAGCGGCGAAACTCGGTTGCCGCTGCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                                                                                                                            <  1:85209/236‑1 (MQ=255)
                                                                                                         aGCCCCCGTTTCTTTTTGGGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                                                                                                                  <  1:1062448/239‑1 (MQ=255)
                                                                                                                            gcgcgggggggaaaGCGTTTGCTGTAGCGGCGAACCCCGTTTCCCTCTTCTCCCAGCGCCCCGTCCAGTCGTTCGCGTTCCAGAGCCGGATCCGGCGTTGAGGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCGGCCAGCCTTCGGCAATCGCCACGCCAGGTGCGGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                              <  2:56631/228‑1 (MQ=255)
                                                                                                                            gccccGGGGGCAAAGCGTTTGTTGTGGGGGCCAAATCCTTTTTCCGTTTCTTCCAGCGCCCCGGTGAGTCGGTCGGGTGCAAGAGCCGGATCCAGGGTTGATGCCTGATACACCAGTTCCATTAAAGGTACCGTGGCATCCTGCCAGCCTTCGGCAATTGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                              <  2:1139700/236‑1 (MQ=255)
                                                                                                                                       aaaGCGTTTGCTGGAGCGGCCAATCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTGCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATCCCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCGTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCt                                                                                                    <  1:654057/238‑1 (MQ=255)
                                                                                                                                            gtttgtgttcgcgacgAACATCTTTTGCCGCTTTTTCCGGCCCCCCGGTCAGCCCTTCGTGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACGCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTGCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATc                                                                                              <  2:562472/229‑1 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCCGCGCCCCGGTCAGTCGTTCGAGGTCCCGAGCCGCATCCAGCGTTGCTGCCTGATACCCCTGTTCCATTAAAGGCAACGTGGCATCCTGCCAGCCTTCCGCCATCGCCACACCAAGAGCTGACGGTCATGCCGGGAATCTCGGCTGCCGATATTGCCCCACCAACGCAGACAACTGCTACGGagaaagaaaaacagcat                                                                                                >  2:5396/1‑227 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCCGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGTGATTGCCCAAACAAGCGAGTCAACTGCGACTGAGAAAGAATAACTGCCa                                                                                                >  1:382798/1‑236 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGGAAAGGGGCATCCTCCCAGCCTTCCGCAATCACCACACCCGGTGTTGCAGGCAATGAGGGGATTCGGGGCTGGCGATATTGCCCAACCCAGGCAGTCACCGCGCACTGagaaagaaaaacggccc                                                                                                >  2:939845/1‑227 (MQ=255)
                                                                                                                                                     aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGGGGCa                                                                                     >  1:1219787/1‑239 (MQ=255)
                                                                                                                                                          cgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCGGCAACCTGGGTGGCAAGGGt                                                                                >  1:134949/1‑239 (MQ=255)
                                                                                                                                                                 cGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGGCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagg                                                                         <  1:380559/239‑1 (MQ=255)
                                                                                                                                                                   ttttCCGTTTCTCCCGGCGCCCCGGTCCGTCTGTCGCGTTCCGGAGCCGGATCCGGCGTTGAGGCCTGATCCCCCTGTTCAATTAAAGGTAACGTGGCAGCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  2:924538/239‑1 (MQ=255)
                                                                                                                                                                   tttGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  1:1209132/239‑1 (MQ=255)
                                                                                                                                                                    ttGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGCGAAAGACAAGCTGCCACCTGGGGGGCAGGCGCAACCaggaag                                                                      >  1:1104803/1‑239 (MQ=255)
                                                                                                                                                                    ttGCCGCTTCTTCCAGCGCCCCGGTCAATCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaag                                                                      >  1:728491/1‑239 (MQ=255)
                                                                                                                                                                      gCCCCTTCTTCCCCCACCCCGGTGACCCCTTCGCGGTCCAGGCGCGGACCCCGAGTTTGTGCCTGAACCCCTGGTTCGTTTAAAGGTAACGGGGCACTCTCCAAGCTTGCGGAAATCGGCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCCACGGGCGATATTGCCCAAACAAGGCATTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGTTGGCAAGCGTCACCaggaagga                                                                    <  2:1161628/239‑1 (MQ=255)
                                                                                                                                                                            tcttcCAGCGCCCCGGTCAGTCGTTTGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGACGTCAACTCCGACTGAGAAAGAGTAGCTTCCATCTGGGTGGCAACGGTGACCaggaaggtgtctc                                                               >  1:533069/1‑232 (MQ=255)
                                                                                                                                                                              tcccAGCCCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCAc                                                            <  1:1112422/237‑1 (MQ=255)
                                                                                                                                                                                ccAGCGCCCCGGGCAGTCGTTCGCGTTCAAGGGCGGGATCCCAGGTTGAGGCCTGATACACCTGTCCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtt                                                          <  2:464027/239‑1 (MQ=255)
                                                                                                                                                                                 cAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACttt                                                         <  1:159466/239‑1 (MQ=255)
                                                                                                                                                                                   gcgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACGCCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  1:130456/239‑1 (MQ=255)
                                                                                                                                                                                    ggcCCCCGGCAGTTCGGTCGGGTCCAAGAGCCGGGTCCGGTGTTGATGCCTGAAACACCCTGTCCATTAAAGGGACCGTGGCATCCTGCCGGCCTGCGGCAAGCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  2:53215/238‑1 (MQ=255)
                                                                                                                                                                                       cccGGTACGTCGTTGGCTTTCCAGGGCGGGATCCAGCGTTGAGTCCTGTGACCCCGGGTGCATTAAGGGTAACGTGGCATCCTGCCGGCCTCCGGCATACGCCACCCCAGGTGCTGGCGGTAATGCGCGGATTCGCGACTGGCGATATTGCCCAATCGAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTGCACTTTCGt                                                      <  2:382798/236‑1 (MQ=255)
                                                                                                                                                                                                acggttgccgTTTCAGAGGCAGGAACCCGCGTTAATGCCTGATCCCCCTGTTCCATTAAAGGGAACGGGGCACCCTGCCAGCCTTCTGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCGCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCTCCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                            <  2:709286/231‑1 (MQ=255)
                                                                                                                                                                                                  cGTTCGGGTTCCAGAGCCGGATCCAGCGTGGATGCCTGATCCCCCTGTTCCATTAAGGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCa                                          <  2:882410/238‑1 (MQ=255)
                                                                                                                                                                                                      cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGTCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCACCTGCGACTGAGAAAGAATAGCTGCCATCTTGGTGGCAAGCGTCACCAGGAAGGCTTCTTCCCTTTCGTCATACTGGCCCACctc                                      >  1:1129020/1‑238 (MQ=255)
                                                                                                                                                                                                      cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGATTAGCTGCCATCTGGGTGGCAAGCGCCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                     >  1:43456/1‑239 (MQ=255)
                                                                                                                                                                                                         gTTCCAGACCCGGCTCCGGCTTTGATGCTTGATACCCCTGTTCCATTAAAGGTAACGGGGCATCCGGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                  <  2:436198/239‑1 (MQ=255)
                                                                                                                                                                                                            ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTAACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTg                                >  1:882410/1‑238 (MQ=255)
                                                                                                                                                                                                             cggagCCGGCTCCCGGTTTGATGCCTGATACCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCt                              <  1:509362/237‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGATCCAGCGTTGATGCCTGATACCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  1:455317/238‑1 (MQ=255)
                                                                                                                                                                                                                    cGGATCCAGCGGTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  1:5263/237‑1 (MQ=255)
                                                                                                                                                                                                                        tCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGGATAACTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGACATACTGGCGCAACCCCCGTTGCGGTacaaacac                    >  2:646037/1‑233 (MQ=255)
                                                                                                                                                                                                                        tCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                     <  1:742631/238‑1 (MQ=255)
                                                                                                                                                                                                                         ccAGCGTTGATTCCTGATACCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAg                    <  1:406947/238‑1 (MQ=255)
                                                                                                                                                                                                                          cAGCGTTGATGCCTTATACGCCGGTTCCATTAAAGTTAACGTGGCATCCTGCCAGCCTTCGGAAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                  <  1:1308562/239‑1 (MQ=255)
                                                                                                                                                                                                                            gCGTTGGTGCCTGATACCCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                <  1:435556/239‑1 (MQ=255)
                                                                                                                                                                                                                             cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaccc               >  1:85319/1‑237 (MQ=255)
                                                                                                                                                                                                                             cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAGTTTTCCCTTTCGTCATACTGGCGCCACTCCCGTTGCGGGACCCCCACAacac               >  2:998392/1‑239 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATCTTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAAAGCTGCCATCTGGGAGGCAGGCGTCACCAGGAAGGATTCTTCACTTTCGTCATAAAGGGGCAACTCTCGTTGCCGTACCACCAGTAACCaaacc         >  2:1172694/1‑236 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAct                                       >  1:946905/1‑210 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTCCTGTACCACCAGACCACCAAGc          >  1:546545/1‑239 (MQ=255)
                                                                                                                                                                                                                                        tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGGCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGACCCACCAGAAAACCaacaaacgg     >  2:348963/1‑232 (MQ=255)
                                                                                                                                                                                                                                            ccacCTGTTCCATTAAGGGTAACGTGGCCTCTGGCCAGGCTTCGGGAAGCGCCACAACAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCATACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCg  <  1:796196/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967611‑2968080

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