Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,557,962 T→C K315K (AAA→AAG)  unknown ← Phosphoglycerate transport system transcriptional regulatory protein PgtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,557,9620TC100.0% 163.5 / NA 57K315K (AAA→AAG) unknownPhosphoglycerate transport system transcriptional regulatory protein PgtA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/28);  total (29/28)

CCGTTTACGCGCGCTACACCCGCCGCGGCGGGTTGGACATCAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATTTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGATGGGTTTATTTAATCT  >  1/1557718‑1558207
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ccGTTTACGCGCGCTACCCCCGCCGCGGCGGGTTGGACATCAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTcccc                                                                                                                                                                                                                                                 <  1:464258/251‑1 (MQ=255)
           cgcTACACCCGCCGCGGCGGGTTGGACGTCAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATaaaaaa                                                                                                                                                                                                                                       >  1:47028/1‑250 (MQ=255)
           cgcTACACCCGCCGCGGCGGGTTGGACATCAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATaaaaaa                                                                                                                                                                                                                                       >  1:629500/1‑250 (MQ=255)
                                        cAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGCAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCat                                                                                                                                                                                                         <  1:265690/251‑1 (MQ=255)
                                         aaCCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAgt                                                                                                                                                                                                                   <  1:99525/240‑1 (MQ=255)
                                         aaCCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCata                                                                                                                                                                                                        >  1:197851/1‑251 (MQ=255)
                                         aaCCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCata                                                                                                                                                                                                        >  1:212675/1‑251 (MQ=255)
                                          aCCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAAt                                                                                                                                                                                                                       >  1:117882/1‑235 (MQ=255)
                                          aCCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCatat                                                                                                                                                                                                       >  1:219872/1‑251 (MQ=255)
                                           cccATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATAtt                                                                                                                                                                                                      <  1:210730/251‑1 (MQ=255)
                                               tACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGAtt                                                                                                                                                                                                  <  1:749168/251‑1 (MQ=255)
                                               tACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGAtt                                                                                                                                                                                                  >  1:330636/1‑251 (MQ=255)
                                                     tCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTg                                                                                                                                                                                            <  1:672412/251‑1 (MQ=255)
                                                         aCTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATa                                                                                                                                                                                                                      <  1:70555/221‑1 (MQ=255)
                                                                        caccGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATc                                                                                                                                                                         <  1:33094/251‑1 (MQ=255)
                                                                        caccGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATc                                                                                                                                                                         <  1:46564/251‑1 (MQ=255)
                                                                              aTTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGaa                                                                                                                                                                   >  1:799605/1‑251 (MQ=255)
                                                                                              gCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGACATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACaggag                                                                                                                                                   >  1:794992/1‑251 (MQ=255)
                                                                                                            cGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGt                                                                                                                                      <  1:218402/250‑1 (MQ=255)
                                                                                                              gTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAATTGTTTTTCat                                                                                                                                                                                                         >  1:430998/1‑181 (MQ=255)
                                                                                                                              gATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGGTTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCg                                                                                                                   >  1:450105/1‑251 (MQ=255)
                                                                                                                              gATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTc                                                                                                                                                                                         >  1:495525/1‑181 (MQ=255)
                                                                                                                                    gACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTg                                                                                                                                                                                            >  1:234608/1‑172 (MQ=255)
                                                                                                                                          aCTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGt                                                                                                       <  1:874333/251‑1 (MQ=255)
                                                                                                                                          aCTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGt                                                                                                       <  1:359274/251‑1 (MQ=255)
                                                                                                                                                  atttataatttataACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAAt                                                                                                                                                                                                                       <  1:689015/131‑1 (MQ=255)
                                                                                                                                                         atttataaCTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAac                                                                                                                                               <  1:361005/196‑1 (MQ=255)
                                                                                                                                                         atttataaCTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCACTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAAt                                                                                                                                                          <  1:607133/185‑1 (MQ=255)
                                                                                                                                                           ttataaCTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCaa                                                                                      >  1:70054/1‑251 (MQ=255)
                                                                                                                                                                   tttCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATaa                                                                                                                                                                                                                     <  1:658689/116‑1 (MQ=255)
                                                                                                                                                                    ttCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACaggag                                                                                                                                                   <  1:476183/181‑1 (MQ=255)
                                                                                                                                                                    ttCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGt                                                                                                       <  1:131133/225‑1 (MQ=255)
                                                                                                                                                                    ttCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTca                                                                              >  1:715235/1‑250 (MQ=255)
                                                                                                                                                                    ttCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTca                                                                              >  1:11615/1‑250 (MQ=255)
                                                                                                                                                                      cTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACAt                                                                           <  1:9069/251‑1 (MQ=255)
                                                                                                                                                                              aCCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATc                                                                                        >  1:325441/1‑230 (MQ=255)
                                                                                                                                                                                   tatTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCAtt                                                               >  1:287133/1‑250 (MQ=255)
                                                                                                                                                                                      tttttCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAgtt                                                                                                                                                                                                                  >  1:490146/1‑100 (MQ=255)
                                                                                                                                                                                             tGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATg                                                                                                                                >  1:43496/1‑175 (MQ=255)
                                                                                                                                                                                                 taagtaaAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCttt                                                           >  1:657869/1‑240 (MQ=255)
                                                                                                                                                                                                         gTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCgcaaccaacaaatt                                        >  1:561246/1‑238 (MQ=255)
                                                                                                                                                                                                                   gggAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGt                                                                                                       <  1:77993/178‑1 (MQ=255)
                                                                                                                                                                                                                    ggAATATTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTg                                                       <  1:553425/225‑1 (MQ=255)
                                                                                                                                                                                                                       atatTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATAc                                                                                                                                                        >  1:831724/1‑125 (MQ=255)
                                                                                                                                                                                                                       atatTGAGGTGTTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCgg                          >  1:251820/1‑251 (MQ=255)
                                                                                                                                                                                                                               gtgtTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGa                  >  1:899136/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gtTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATc                                                                                                                                                                 <  1:523203/106‑1 (MQ=255)
                                                                                                                                                                                                                                 gtTTAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATc                                                                                        >  1:567399/1‑179 (MQ=255)
                                                                                                                                                                                                                                    tAACTACGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAAt                                                                                                                                                          <  1:702060/110‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGATGGGtttattt       <  1:480977/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCatat                                                                                                                                                                                                       <  1:119266/57‑1 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAg                                                                                                                                                                                 <  1:868352/79‑1 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTa                                                                                                                        >  1:179013/1‑136 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAAtt                                               >  1:880468/1‑209 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGa                  <  1:855644/238‑1 (MQ=255)
                                                                                                                                                                                                                                            tCATTAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGATGGGTTTATTTaa     >  1:400613/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ttAATCTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGATGGGTTTATTTAATCt  <  1:560913/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
CCGTTTACGCGCGCTACACCCGCCGCGGCGGGTTGGACATCAACCCATACCGCTCCAACTTTGAACCTTTACCACCGAATTTACGACTGGCACGGCAGTAAAAAAGTGCGGTTAGAAAACGCCATGGATTTTGACCGCACTTTTCGATTTATAATTTATAACTTTTCTTATCCAACCCATATTTTTTCATGCGTAAGTAAAGTTTCTTTCTGGGAATATTGAGGTGTTTAACTACGTCATTAATTTTCCCCTGATAAAAAACCAGAGCATCTTCAATAAGTTGTTTTTCATATTGATTCAATTGTTCATCTAAAGAACTCATCTGCGAATCATTAAATACAGGAGGAACCACCAAAGTAATATGAGAAAGTAAGCCGATTGCATAAAGTTCCGCGACGTTAATCAATTCTTTCACATTTCCATGCCATTCTTTATTGCGTAAATTTTGCAACATTTTTTTATCCGGGATTACGATGGGTTTATTTAATCT  >  1/1557718‑1558207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑