Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,248,009 G→T 8.2% I219I (ATC→ATA rihB ← ribonucleoside hydrolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,248,0090GT 8.2% 744.1 / 27.9 207I219I (ATC→ATArihBribonucleoside hydrolase 2
Reads supporting (aligned to +/- strand):  ref base G (101/89);  new base T (9/8);  total (110/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 219 total aligned reads displayed.

AGGGCCACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCTTCCGGG  >  CP009273/2247862‑2248156
                                                                                                                                                   |                                                                                                                                                   
aGGGCCACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCatg                                                                                                                                                     <  1:2160314/148‑1 (MQ=255)
aGGGCCACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCatg                                                                                                                                                     >  1:2014769/1‑148 (MQ=255)
  ggCCACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCatgat                                                                                                                                                   >  1:1817646/1‑148 (MQ=255)
    ccACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGt                                                                                                                                                 <  1:1172824/148‑1 (MQ=255)
    ccACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGt                                                                                                                                                 <  2:2403898/148‑1 (MQ=255)
      aCTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCg                                                                                                                                               <  2:2143124/148‑1 (MQ=255)
      aCTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCg                                                                                                                                               <  2:478995/148‑1 (MQ=255)
         gTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTg                                                                                                                                            <  2:3010516/148‑1 (MQ=255)
         gTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTg                                                                                                                                            >  1:2941422/1‑148 (MQ=255)
            gACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAAc                                                                                                                                         <  2:2569135/148‑1 (MQ=255)
             aCGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACa                                                                                                                                        >  2:1033556/1‑148 (MQ=255)
             aCGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACa                                                                                                                                        >  1:2143073/1‑148 (MQ=255)
              cgaccACTTCGACGTACATCTCCTGGGTTTTAATGCCATCACGGTTAATCAGATAACCGAGGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAg                                                                                                                                       <  1:3067623/145‑1 (MQ=255)
                  cACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGctct                                                                                                                                   <  1:3177719/148‑1 (MQ=255)
                   aCTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCtctc                                                                                                                                  <  1:1931831/148‑1 (MQ=255)
                    cTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTcc                                                                                                                                 >  1:730422/1‑148 (MQ=255)
                    cTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTcc                                                                                                                                 <  2:2407686/148‑1 (MQ=255)
                      tCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCcgg                                                                                                                               <  1:831332/148‑1 (MQ=255)
                       cGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCcggc                                                                                                                              <  2:2345598/148‑1 (MQ=255)
                       cGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCcggc                                                                                                                              <  2:1227884/148‑1 (MQ=255)
                           gTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCgggg                                                                                                                          >  1:1206288/1‑148 (MQ=255)
                             aCATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGgcc                                                                                                                        >  1:2119392/1‑148 (MQ=255)
                             aCATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGgcc                                                                                                                        >  2:1768897/1‑148 (MQ=255)
                                tctcCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGgccgcc                                                                                                                     >  1:349847/1‑148 (MQ=255)
                                tctcCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGgccgcc                                                                                                                     >  1:349565/1‑148 (MQ=255)
                                tctcCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGgccgcc                                                                                                                     <  2:1739893/148‑1 (MQ=255)
                                 ctcCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCt                                                                                                                    >  2:33470/1‑148 (MQ=255)
                                      gggTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCt                                                                                                               >  1:2465822/1‑148 (MQ=255)
                                      gggTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCt                                                                                                               >  1:2611663/1‑148 (MQ=255)
                                      gggTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCt                                                                                                               >  1:2831257/1‑148 (MQ=255)
                                        gTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCttt                                                                                                             <  2:1271192/148‑1 (MQ=255)
                                         ttttAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTc                                                                                                            <  1:128358/148‑1 (MQ=255)
                                          tttAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGGTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTcc                                                                                                           >  2:1647623/1‑148 (MQ=255)
                                            tAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCAt                                                                                                         <  2:2034543/148‑1 (MQ=255)
                                                   aTCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCa                                                                                                  <  1:2219842/148‑1 (MQ=255)
                                                   aTCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCa                                                                                                  <  2:2611663/148‑1 (MQ=255)
                                                     cAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAAt                                                                                                >  1:1671669/1‑148 (MQ=255)
                                                     cAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAAt                                                                                                >  1:1010841/1‑148 (MQ=255)
                                                     cAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCATTTCAATCAGAGCAAt                                                                                                >  1:455091/1‑148 (MQ=255)
                                                        ggTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCAc                                                                                             <  2:1708559/148‑1 (MQ=255)
                                                         gTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACg                                                                                            >  1:683581/1‑148 (MQ=255)
                                                             aTCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCg                                                                                        >  2:930201/1‑148 (MQ=255)
                                                             aTCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCg                                                                                        >  2:1618072/1‑148 (MQ=255)
                                                              tCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCgg                                                                                       >  2:1998674/1‑148 (MQ=255)
                                                                aGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCgggg                                                                                     <  2:860922/148‑1 (MQ=255)
                                                                  aTAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACCGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTg                                                                                   >  2:703304/1‑148 (MQ=255)
                                                                    aaCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCa                                                                                 >  2:2620198/1‑148 (MQ=255)
                                                                    aaCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCa                                                                                 >  1:614551/1‑148 (MQ=255)
                                                                    aaCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCa                                                                                 >  2:535843/1‑148 (MQ=255)
                                                                     aCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTcc                                                                                                           <  2:1344026/121‑1 (MQ=255)
                                                                     aCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCaa                                                                                <  2:1444439/148‑1 (MQ=255)
                                                                     aCCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCaa                                                                                >  2:2930972/1‑148 (MQ=255)
                                                                      ccGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCaaa                                                                               >  1:2935573/1‑148 (MQ=255)
                                                                          tGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGt                                                                           <  1:1742657/148‑1 (MQ=255)
                                                                                  gCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCtggttgg                                                                   <  2:2107124/148‑1 (MQ=255)
                                                                                  gCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCtggttgg                                                                   <  2:2107061/148‑1 (MQ=255)
                                                                                    gtcgtGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTg                                                                 <  2:1114359/148‑1 (MQ=255)
                                                                                       gtGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAACCACGTCCGGGGTGCAAACGGTCTGGTTGGTgaga                                                              <  1:3087450/148‑1 (MQ=255)
                                                                                          cACCGGGCCGCCAGCAAGGCCGTAGTTGTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCg                                                           <  2:2656694/148‑1 (MQ=255)
                                                                                           aCCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGa                                                          <  1:950860/148‑1 (MQ=255)
                                                                                           aCCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGa                                                          <  1:1447645/148‑1 (MQ=255)
                                                                                            ccGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAg                                                         >  1:3260050/1‑148 (MQ=255)
                                                                                             cGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAgg                                                        >  1:2221571/1‑148 (MQ=255)
                                                                                             cGGGCCGCCAGCAAGGCCGTAATTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAgg                                                        >  1:1726892/1‑148 (MQ=255)
                                                                                               ggCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGcc                                                      <  2:953362/148‑1 (MQ=255)
                                                                                                  cgccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCcat                                                   <  1:972399/148‑1 (MQ=255)
                                                                                                  cgccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCcat                                                   <  2:1991674/148‑1 (MQ=255)
                                                                                                  cgccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCcat                                                   <  2:1991970/148‑1 (MQ=255)
                                                                                                     cAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCcatcat                                                <  1:851373/148‑1 (MQ=255)
                                                                                                          aGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCAATCATCACTa                                           >  2:1505629/1‑148 (MQ=255)
                                                                                                           ggCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTaa                                          <  2:2790037/148‑1 (MQ=255)
                                                                                                            gCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAAt                                         >  2:1827734/1‑148 (MQ=255)
                                                                                                            gCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAAt                                         >  1:605724/1‑148 (MQ=255)
                                                                                                            gCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTGGGTGAGATCGAGGCCCATCATCACTAAt                                         >  1:1730043/1‑148 (MQ=255)
                                                                                                             ccGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATg                                        >  1:1025025/1‑148 (MQ=255)
                                                                                                             ccGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATg                                        >  1:2678961/1‑148 (MQ=255)
                                                                                                               gTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGa                                      <  1:2049973/148‑1 (MQ=255)
                                                                                                                tAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGaa                                     <  1:1592925/148‑1 (MQ=255)
                                                                                                                      tCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCgg                               <  1:1341977/148‑1 (MQ=255)
                                                                                                                           cTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTg                          >  2:2421562/1‑148 (MQ=255)
                                                                                                                           cTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTg                          >  1:1745723/1‑148 (MQ=255)
                                                                                                                           cTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTg                          >  2:3183369/1‑148 (MQ=255)
                                                                                                                            tGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGa                         >  1:29823/1‑148 (MQ=255)
                                                                                                                            tGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGa                         <  2:310662/148‑1 (MQ=255)
                                                                                                                             gCGTTTTGAGAGTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGaa                        >  2:659058/1‑148 (MQ=255)
                                                                                                                              cGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTACCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAac                       <  1:33470/148‑1 (MQ=255)
                                                                                                                                      gagTGAAGTTCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTAcgcg               <  2:3149617/148‑1 (MQ=255)
                                                                                                                                       agTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGt              >  2:2335297/1‑148 (MQ=255)
                                                                                                                                       agTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGt              >  1:2720802/1‑148 (MQ=255)
                                                                                                                                       agTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGt              >  2:944008/1‑148 (MQ=255)
                                                                                                                                       agTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGt              >  2:120968/1‑148 (MQ=255)
                                                                                                                                          gAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCg           >  1:234261/1‑148 (MQ=255)
                                                                                                                                            aGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGc         <  2:614532/148‑1 (MQ=255)
                                                                                                                                             gTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCt        >  2:2081817/1‑148 (MQ=255)
                                                                                                                                              ttCATTATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCtt       >  1:288262/1‑148 (MQ=255)
                                                                                                                                              ttCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCtt       <  1:755588/148‑1 (MQ=255)
                                                                                                                                              ttCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCtt       >  2:2743631/1‑148 (MQ=255)
                                                                                                                                              ttCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCtt       <  1:2417591/148‑1 (MQ=255)
                                                                                                                                                  tgatgTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCTTCCgg   >  1:197966/1‑148 (MQ=255)
                                                                                                                                                   gatgTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCTTCCggg  >  1:163135/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGGCCACTGTTGACGTCCACTTCGACGTACATCTCCTGGGTTTTAATGCCATCAGGGTTAATCAGATAACCGATGCAGGTGGCGTCGTGCACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGAGGCCCATCATCACTAATGGAACGCCGGAGGTGAACACTACGCGTGCGGCTTCCGGG  >  CP009273/2247862‑2248156

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