Predicted mutation
evidence seq id position mutation annotation gene description
RA 2 120,670 G→A D146D (GAC→GAT)  unknown ← 16S rRNA (guanine(527)‑N(7))‑methyltransferase (EC 2.1.1.170)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*2120,6700GA100.0% 186.9 / NA 70D146D (GAC→GAT) unknown16S rRNA (guanine(527)‑N(7))‑methyltransferase (EC 2.1.1.170)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/27);  total (43/27)

CAAAAGCACTCAAAAAGTGGTAAATGAACTACAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATGTCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAATAAAACGATCACCTTG  >  2/120434‑120915
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
cAAAAGCACTCAAAAAGTGGTAAATGAACTACAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAg                                                                                                                                                                                                                                          >  1:636483/1‑250 (MQ=255)
        cTCAAAAAGTGGTAAATGAACTACAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAgcg                                                                                                                                                                                                                                  >  1:59752/1‑250 (MQ=255)
                              aCAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAAttt                                                                                                                                                                                                           >  1:333672/1‑251 (MQ=255)
                               cAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAAtttt                                                                                                                                                                                                          >  1:187423/1‑251 (MQ=255)
                                   aaaaGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAgcgc                                                                                                                                                                                                                                 <  1:467110/224‑1 (MQ=255)
                                        aTAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGc                                                                                                                                                                                                                     >  1:322702/1‑231 (MQ=255)
                                              aaCAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGatat                                                                                                                                                                                           <  1:883042/251‑1 (MQ=255)
                                                    aaTTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGc                                                                                                                                                                                                                     <  1:158624/219‑1 (MQ=255)
                                                                 aGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGaa                                                                                                                                                                        >  1:313162/1‑251 (MQ=255)
                                                                   tCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       <  1:751653/250‑1 (MQ=255)
                                                                   tCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       <  1:453481/250‑1 (MQ=255)
                                                                   tCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       <  1:432388/250‑1 (MQ=255)
                                                                   tCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       <  1:255695/250‑1 (MQ=255)
                                                                         gTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTa                                                                                                                                                                >  1:669448/1‑251 (MQ=255)
                                                                           aTTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTaa                                                                                                                                                               >  1:656637/1‑250 (MQ=255)
                                                                            ttAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTcc                                                                                                                                                                                     >  1:690269/1‑227 (MQ=255)
                                                                                  ggCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGtt                                                                                                                                                       >  1:623874/1‑251 (MQ=255)
                                                                                   gCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTa                                                                                                                                                      >  1:174309/1‑251 (MQ=255)
                                                                                       tCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTcc                                                                                                                                                  <  1:838734/251‑1 (MQ=255)
                                                                                       tCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTc                                                                                                                                                   >  1:513287/1‑250 (MQ=255)
                                                                                       tCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTc                                                                                                                                                   >  1:667265/1‑250 (MQ=255)
                                                                                          ccGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       <  1:398383/227‑1 (MQ=255)
                                                                                          ccGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCaa                                                                                                                                               >  1:19948/1‑251 (MQ=255)
                                                                                          ccGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCaa                                                                                                                                               >  1:640941/1‑251 (MQ=255)
                                                                                             aTAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTc                                                                                                                                            <  1:46141/251‑1 (MQ=255)
                                                                                                 gtgtTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAGCATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAgg                                                                                                                                                                                                                                             <  1:477393/150‑1 (MQ=255)
                                                                                                 gtgtTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATAc                                                                                                                                                                                                                       >  1:161629/1‑172 (MQ=255)
                                                                                                 gtgtTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGatat                                                                                                                                                                                           >  1:707957/1‑200 (MQ=255)
                                                                                                 gtgtTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATAttcttc                                                                                                                                                                                      >  1:623568/1‑205 (MQ=255)
                                                                                                 gtgtTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       >  1:573790/1‑220 (MQ=255)
                                                                                                             cAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACg                                                                                                                                                             <  1:568400/218‑1 (MQ=255)
                                                                                                             cAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGa                                                                                                                            >  1:238670/1‑251 (MQ=255)
                                                                                                                  gAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATaaa                                                                                                                       >  1:397240/1‑251 (MQ=255)
                                                                                                                    aaTCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGc                                                                                                                     <  1:69739/251‑1 (MQ=255)
                                                                                                                        acaTTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATAC‑CCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATa                                                                                                                <  1:364036/251‑1 (MQ=255)
                                                                                                                            tttCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGc                                                                                                             <  1:67773/251‑1 (MQ=255)
                                                                                                                                              cACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCaa                                                                                            >  1:257983/1‑250 (MQ=255)
                                                                                                                                                cTTAAATTTTGTGCTTCCTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACg                                                                                                              >  1:175297/1‑230 (MQ=255)
                                                                                                                                                   aaaTTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTAc                                                                                      <  1:204335/251‑1 (MQ=255)
                                                                                                                                                     aTTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAg                                                                                         >  1:787707/1‑246 (MQ=255)
                                                                                                                                                        ttGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACt                                                                                 <  1:162047/251‑1 (MQ=255)
                                                                                                                                                         tgtgCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTg                                                                                <  1:26937/251‑1 (MQ=255)
                                                                                                                                                          gtgCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTg                                                                                >  1:292184/1‑250 (MQ=255)
                                                                                                                                                          gtgCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTg                                                                                >  1:194214/1‑250 (MQ=255)
                                                                                                                                                                  tcCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTATCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAgg                                                                       <  1:311552/251‑1 (MQ=255)
                                                                                                                                                                         gATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGa                                                                                                                                                                         <  1:39112/146‑1 (MQ=255)
                                                                                                                                                                         gATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       >  1:463562/1‑148 (MQ=255)
                                                                                                                                                                         gATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACt                                                                                                                                             <  1:127558/174‑1 (MQ=255)
                                                                                                                                                                         gATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGaa                                                                                                                           <  1:692192/192‑1 (MQ=255)
                                                                                                                                                                         gATAAATACCTTTTAACGCATAAAACATACAATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATg                                                                >  1:95468/1‑251 (MQ=255)
                                                                                                                                                                          aTAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGa                                                               <  1:856189/251‑1 (MQ=255)
                                                                                                                                                                          aTAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGa                                                               <  1:514042/251‑1 (MQ=255)
                                                                                                                                                                                 ccTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAg                                                                        >  1:331904/1‑235 (MQ=255)
                                                                                                                                                                                         cGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAg                                                 >  1:64856/1‑250 (MQ=255)
                                                                                                                                                                                          gCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGc                                                >  1:417249/1‑250 (MQ=255)
                                                                                                                                                                                              aaaaCATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCggg                                            >  1:829858/1‑250 (MQ=255)
                                                                                                                                                                                               aaaCATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAGTAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACaa                                                                                    >  1:443548/1‑209 (MQ=255)
                                                                                                                                                                                                            ggTTTGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAg                                                                                         >  1:139550/1‑191 (MQ=255)
                                                                                                                                                                                                                tGGGTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCt                         <  1:228984/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATa                                                                                                                         >  1:54714/1‑152 (MQ=255)
                                                                                                                                                                                                                   gTAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTAcc                          >  1:170155/1‑247 (MQ=255)
                                                                                                                                                                                                                    tAAATGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACa                      <  1:201050/250‑1 (MQ=255)
                                                                                                                                                                                                                        tGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATg                                                                >  1:662620/1‑204 (MQ=255)
                                                                                                                                                                                                                        tGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAAt                 >  1:813839/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGCCAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAAt                 >  1:562596/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccAGCACCATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAATaa               >  1:900035/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccATTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAAc                                                                                                                                                                       >  1:434932/1‑93 (MQ=255)
                                                                                                                                                                                                                                  aTTCTGTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAATAAAACGATCa       >  1:554870/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTAc                           <  1:621068/226‑1 (MQ=255)
                                                                                                                                                                                                                                       gTCATATCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAATAAAACGATCACCTTg  >  1:552649/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
CAAAAGCACTCAAAAAGTGGTAAATGAACTACAATAAAAGATAGTTAACAGAAATTTTACTTTTTAGTCAAAAGTATTAAATGGCGTTCACCGATAAGTGTTGGAACCTCAATAGAAATCACATTTTCAATCTGAAATTGATCACTTAAATTTTGTGCTTCTTCCTGCTGATAAATACCTTTTAACGCATAAAACATACCATCAGGTTTGGGTAAATGCCAGCACCATTCTGTCATGTCTTTTAAGGAAGCAAAAGCGCGGCTTAATACGCCGTCGAATTTTTCTTCTGGTTGATATTCTTCCACACGGCTTAGAACAGGCTTAACGTTTGTTAGTCCCAACTCACGAACCGCATTACGAATAAAGCTGATACGCTTGCCCAGGCTATCCAACAGTACAAACTGCTTATCAGGGTTAATGATCGCCAACGGCAAGCCGGGCAAGCCAGGACCGGTACCTACATCAATAAAACGATCACCTTG  >  2/120434‑120915

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