Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,208,550 T→A 5.3% K95* (AAA→TAA)  mug ← G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,208,5500TA 5.3% 1006.8 / 33.4 303K95* (AAA→TAA) mugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (149/138);  new base A (10/6);  total (159/144)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 311 total aligned reads displayed.

AAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGT  >  CP009273/3208402‑3208694
                                                                                                                                                    |                                                                                                                                                
  aCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAAtttttt                                                                                                                                                >  1:4123400/1‑149 (MQ=255)
      cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAAt                                                                                                                                            <  2:2655929/149‑1 (MQ=255)
       tGGGATTTTGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGACCAGTATAGCCAACGCCTGCTGCTGTTAATCTTCAATTTTTTCAATc                                                                                                                                           <  1:76201/149‑1 (MQ=255)
       tGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATc                                                                                                                                           >  1:1294072/1‑149 (MQ=255)
           aTTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCt                                                                                                                                       <  2:4787570/149‑1 (MQ=255)
           aTTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCt                                                                                                                                       <  1:570523/149‑1 (MQ=255)
           aTTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCt                                                                                                                                       <  1:1673181/149‑1 (MQ=255)
            tttGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCtt                                                                                                                                      >  2:494710/1‑149 (MQ=255)
               ggCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAcg                                                                                                                                   >  2:4472474/1‑149 (MQ=255)
                 cagcaCCCAAATCTGCGTCGAACCACTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAcgcc                                                                                                                                 <  2:1117635/149‑1 (MQ=255)
                     aCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAcgcccgcc                                                                                                                             <  1:2044010/149‑1 (MQ=255)
                     aCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAcgcccgcc                                                                                                                             >  1:4135705/1‑149 (MQ=255)
                       ccAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAACTTACGCCCGCCTg                                                                                                                           >  2:3286980/1‑149 (MQ=255)
                               gCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGc                                                                                                                   >  1:2848603/1‑149 (MQ=255)
                                 gTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTc                                                                                                                 >  1:2631016/1‑149 (MQ=255)
                                 gTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAAGCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTc                                                                                                                 >  1:1731968/1‑149 (MQ=255)
                                  tCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTcc                                                                                                                <  1:855569/149‑1 (MQ=255)
                                   cGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCt                                                                                                               <  1:3270702/149‑1 (MQ=255)
                                   cGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCt                                                                                                               <  2:3231069/149‑1 (MQ=255)
                                    gAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTg                                                                                                              >  2:4235053/1‑149 (MQ=255)
                                     aaCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGc                                                                                                             <  2:2793446/149‑1 (MQ=255)
                                       ccAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                           <  2:75716/149‑1 (MQ=255)
                                        cAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCttt                                                                                                          >  1:4112020/1‑149 (MQ=255)
                                         aaTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTg                                                                                                         >  1:4481929/1‑149 (MQ=255)
                                          aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGa                                                                                                        >  1:1626051/1‑149 (MQ=255)
                                            ggTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaaa                                                                                                      >  2:190766/1‑149 (MQ=255)
                                            ggTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaaa                                                                                                      >  2:1027786/1‑149 (MQ=255)
                                              tGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACt                                                                                                    >  2:2203035/1‑149 (MQ=255)
                                               gAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACtt                                                                                                   <  2:3641985/149‑1 (MQ=255)
                                               gAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTCTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACtt                                                                                                   >  1:3674904/1‑149 (MQ=255)
                                                 gCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCa                                                                                                 <  1:2340632/149‑1 (MQ=255)
                                                 gCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTGTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCGTTGAAACTTCa                                                                                                 >  2:414223/1‑149 (MQ=255)
                                                  cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAt                                                                                                <  2:406303/149‑1 (MQ=255)
                                                        tttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCtt                                                                                          >  2:1179004/1‑149 (MQ=255)
                                                        tttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCtt                                                                                          <  1:142509/149‑1 (MQ=255)
                                                          tCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGc                                                                                        >  2:213336/1‑149 (MQ=255)
                                                           cccccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCa                                                                                       <  1:2678771/149‑1 (MQ=255)
                                                           cccccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCa                                                                                       <  1:1896539/149‑1 (MQ=255)
                                                            ccccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCAc                                                                                      >  2:4269618/1‑149 (MQ=255)
                                                                aCTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                  >  2:3169058/1‑149 (MQ=255)
                                                                   gtgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGa                                                                               >  1:985027/1‑149 (MQ=255)
                                                                   gtgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGa                                                                               >  2:4731731/1‑149 (MQ=255)
                                                                    tgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGAc                                                                              <  2:1743436/149‑1 (MQ=255)
                                                                       aCCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGt                                                                           <  1:4016894/149‑1 (MQ=255)
                                                                        ccGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTc                                                                          <  2:3758175/149‑1 (MQ=255)
                                                                         cgcgCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCt                                                                         <  2:2040578/149‑1 (MQ=255)
                                                                           cgcTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTACTTCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTAc                                                                       <  1:1443011/149‑1 (MQ=255)
                                                                            gctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTAcc                                                                      >  1:4595358/1‑149 (MQ=255)
                                                                              tggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                    >  1:2313768/1‑149 (MQ=255)
                                                                              tggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                    <  2:65499/149‑1 (MQ=255)
                                                                                gctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGtt                                                                  <  1:3405110/149‑1 (MQ=255)
                                                                                 ctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGttt                                                                 <  1:282540/149‑1 (MQ=255)
                                                                                  tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTg                                                                >  1:4730502/1‑149 (MQ=255)
                                                                                   gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                               <  1:3168012/149‑1 (MQ=255)
                                                                                   gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                               >  1:398070/1‑149 (MQ=255)
                                                                                    aaTCCCTGTTCATATGCTTGTTTGCCCAGAATCGTCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                              >  1:2553491/1‑149 (MQ=255)
                                                                                    aaTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                              >  2:3629873/1‑149 (MQ=255)
                                                                                            ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcac                                                      >  1:2007753/1‑149 (MQ=255)
                                                                                            ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcac                                                      >  2:1448673/1‑149 (MQ=255)
                                                                                             tCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                     >  2:1328164/1‑149 (MQ=255)
                                                                                              cATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACaa                                                    <  2:3822166/149‑1 (MQ=255)
                                                                                               atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAAc                                                   <  1:1027103/149‑1 (MQ=255)
                                                                                                tatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                  >  1:1790650/1‑149 (MQ=255)
                                                                                                tatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                  >  1:1172115/1‑149 (MQ=255)
                                                                                                 atGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGa                                                 <  2:1681946/149‑1 (MQ=255)
                                                                                                 atGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGa                                                 >  1:657412/1‑149 (MQ=255)
                                                                                                  tGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGAt                                                <  1:494710/149‑1 (MQ=255)
                                                                                                    cTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATaa                                              <  2:4483641/149‑1 (MQ=255)
                                                                                                     ttGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAAt                                             >  2:3682937/1‑149 (MQ=255)
                                                                                                     ttGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAAt                                             >  1:805345/1‑149 (MQ=255)
                                                                                                      tGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATc                                            <  1:17603/149‑1 (MQ=255)
                                                                                                        tttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCca                                          >  1:4181371/1‑149 (MQ=255)
                                                                                                         ttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCcag                                         <  2:1549690/149‑1 (MQ=255)
                                                                                                          tGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCcagc                                        <  1:2427159/149‑1 (MQ=255)
                                                                                                             ccAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                     <  2:1845896/149‑1 (MQ=255)
                                                                                                              cAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGAt                                    >  1:2105972/1‑149 (MQ=255)
                                                                                                               aGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATg                                   >  2:2732513/1‑149 (MQ=255)
                                                                                                                gAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGc                                  <  1:440715/149‑1 (MQ=255)
                                                                                                                 aaTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCt                                 <  2:3458546/149‑1 (MQ=255)
                                                                                                                 aaTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCt                                 >  2:4414315/1‑149 (MQ=255)
                                                                                                                   tCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtgt                               >  2:2451763/1‑149 (MQ=255)
                                                                                                                        aaCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                          <  1:3029587/149‑1 (MQ=255)
                                                                                                                         aCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                         <  1:4269619/149‑1 (MQ=255)
                                                                                                                          cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcct                        <  2:2848603/149‑1 (MQ=255)
                                                                                                                          cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcct                        <  2:1357140/149‑1 (MQ=255)
                                                                                                                            ccTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGc                      >  1:3135986/1‑149 (MQ=255)
                                                                                                                             cTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                     >  1:1012527/1‑149 (MQ=255)
                                                                                                                             cTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                     >  1:1267953/1‑149 (MQ=255)
                                                                                                                              tGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCgg                    >  2:3878597/1‑149 (MQ=255)
                                                                                                                                cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCt                  <  2:2815795/149‑1 (MQ=255)
                                                                                                                                    tGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCaa              >  2:4307401/1‑149 (MQ=255)
                                                                                                                                      aTAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt            >  2:742195/1‑149 (MQ=255)
                                                                                                                                      aTAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt            >  1:1578421/1‑149 (MQ=255)
                                                                                                                                      aTAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt            <  2:350808/149‑1 (MQ=255)
                                                                                                                                            ttCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggt      >  2:4031140/1‑149 (MQ=255)
                                                                                                                                            ttCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggt      >  2:2321428/1‑149 (MQ=255)
                                                                                                                                            ttCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggt      >  2:1980018/1‑149 (MQ=255)
                                                                                                                                             tCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtc     >  1:3521419/1‑149 (MQ=255)
                                                                                                                                              cAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcg    <  1:4537123/149‑1 (MQ=255)
                                                                                                                                               aaTTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcgg   >  1:629621/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGT  >  CP009273/3208402‑3208694

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: