Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,017,759 A→T 23.5% D17V (GAT→GTT)  xanQ → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,017,7590AT23.5% 597.0 / 194.2 263D17V (GAT→GTT) xanQxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (113/88);  new base T (29/33);  total (142/121)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

AATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCGATGATCG  >  CP009273/3017611‑3017908
                                                                                                                                                    |                                                                                                                                                     
aaTGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGAt                                                                                                                                                      <  2:611480/150‑1 (MQ=255)
   gagaTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGc                                                                                                                                                   >  1:4791275/1‑150 (MQ=255)
   gagaTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGc                                                                                                                                                   >  2:3640016/1‑150 (MQ=255)
   gagaTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGc                                                                                                                                                   >  1:4532841/1‑150 (MQ=255)
      aTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGccct                                                                                                                                                >  1:2056855/1‑150 (MQ=255)
      aTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccct                                                                                                                                                <  2:3283490/150‑1 (MQ=255)
        aaGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccctcc                                                                                                                                              >  2:311408/1‑150 (MQ=255)
        aaGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccctcc                                                                                                                                              >  1:4116433/1‑150 (MQ=255)
         aGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccctccc                                                                                                                                             >  2:2488237/1‑150 (MQ=255)
              ttCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTtca                                                                                                                                        >  1:3790837/1‑150 (MQ=255)
              ttCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCACCCTCCCTTtca                                                                                                                                        >  1:1226786/1‑150 (MQ=255)
                ctcAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTtcatc                                                                                                                                      <  2:453277/150‑1 (MQ=255)
                  cAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAg                                                                                                                                    >  2:970979/1‑150 (MQ=255)
                     ggCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGct                                                                                                                                 <  1:3028675/150‑1 (MQ=255)
                     ggCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGct                                                                                                                                 <  1:1971411/150‑1 (MQ=255)
                       cGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGctct                                                                                                                               >  1:901726/1‑150 (MQ=255)
                           cTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGtt                                                                                                                           <  1:3738320/150‑1 (MQ=255)
                               gcCCCAAAGTGGTGGGAATTTATTTTACTACTTCAGATGATCAAATTTACTTTAAAAGGAGTGACGGGAGGTCTGTTATAAACCCGGCAGGGTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTg                                                                                                                       <  1:2679748/150‑1 (MQ=255)
                                 cccAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGcc                                                                                                                     <  2:107486/150‑1 (MQ=255)
                                    aaaTTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCAtt                                                                                                                  >  2:1115944/1‑150 (MQ=255)
                                      aTTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTAc                                                                                                                >  1:915781/1‑150 (MQ=255)
                                      aTTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTAc                                                                                                                <  1:3718327/150‑1 (MQ=255)
                                         cGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCa                                                                                                             <  2:594343/150‑1 (MQ=255)
                                             ggAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTg                                                                                                         <  2:1300855/150‑1 (MQ=255)
                                               aaTTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGtt                                                                                                       >  1:150067/1‑150 (MQ=255)
                                                atttatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTg                                                                                                      <  2:3768505/150‑1 (MQ=255)
                                                 tttatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTgg                                                                                                     >  1:4050109/1‑150 (MQ=255)
                                                 tttatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTgg                                                                                                     <  2:741549/150‑1 (MQ=255)
                                                  ttatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGc                                                                                                    <  1:4737232/150‑1 (MQ=255)
                                                  ttatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGc                                                                                                    >  2:1890687/1‑150 (MQ=255)
                                                   tatttTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCa                                                                                                   >  2:2844701/1‑150 (MQ=255)
                                                       ttACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAAttt                                                                                               >  2:1924184/1‑150 (MQ=255)
                                                        tACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAAtttt                                                                                              >  2:767203/1‑150 (MQ=255)
                                                         aCTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTc                                                                                             >  2:2581957/1‑150 (MQ=255)
                                                          cTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCg                                                                                            <  2:353537/150‑1 (MQ=255)
                                                          cTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCg                                                                                            >  1:3044706/1‑150 (MQ=255)
                                                            aaTTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGtt                                                                                          >  1:2038883/1‑150 (MQ=255)
                                                               tCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCg                                                                                       <  2:816212/150‑1 (MQ=255)
                                                                cAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGa                                                                                      <  1:2508425/150‑1 (MQ=255)
                                                                 aGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGAt                                                                                     <  2:1631190/150‑1 (MQ=255)
                                                                 aGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACATTATATTTGAACTGGAGGTTCGACCTCCCATTCATCAGGCTCTCGTTGGTGCCATTACCCATATGTTGGCAATTTAAGTTCCTAt                                                                                     >  1:3611579/1‑150 (MQ=255)
                                                                     gatCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTa                                                                                 <  1:1311773/150‑1 (MQ=255)
                                                                     gatCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTa                                                                                 <  1:3235134/150‑1 (MQ=255)
                                                                       tCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTCAc                                                                               >  2:1156078/1‑150 (MQ=255)
                                                                         aaaTTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAAccc                                                                             >  2:1440095/1‑150 (MQ=255)
                                                                              tACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCcgcg                                                                        >  2:1709072/1‑150 (MQ=255)
                                                                              tACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCcgcg                                                                        >  1:869203/1‑150 (MQ=255)
                                                                              tACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCcgcg                                                                        >  2:1076939/1‑150 (MQ=255)
                                                                               aCTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGt                                                                       >  2:2807941/1‑150 (MQ=255)
                                                                                 tttAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTa                                                                     <  1:2367445/150‑1 (MQ=255)
                                                                                 tttAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTa                                                                     <  1:4818026/150‑1 (MQ=255)
                                                                                  ttAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTaa                                                                    >  2:2744969/1‑150 (MQ=255)
                                                                                   tAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAAt                                                                   >  1:2888073/1‑150 (MQ=255)
                                                                                    aaaaGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATc                                                                  >  1:4304418/1‑150 (MQ=255)
                                                                                       aGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTg                                                               >  1:3902104/1‑150 (MQ=255)
                                                                                         gAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTggg                                                             <  1:2594760/150‑1 (MQ=255)
                                                                                            tGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGc                                                          <  2:4471100/150‑1 (MQ=255)
                                                                                             gAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCg                                                         >  2:4195047/1‑150 (MQ=255)
                                                                                              aaGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCgg                                                        >  2:1596652/1‑150 (MQ=255)
                                                                                                gggATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGcc                                                      >  1:3274933/1‑150 (MQ=255)
                                                                                                gggATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCCTCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGcc                                                      <  2:1899170/150‑1 (MQ=255)
                                                                                                gggATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGcc                                                      <  1:79319/150‑1 (MQ=255)
                                                                                                gggATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGcc                                                      <  1:3348409/150‑1 (MQ=255)
                                                                                                  gATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCtt                                                    <  2:134507/150‑1 (MQ=255)
                                                                                                  gATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCtt                                                    >  2:4054042/1‑150 (MQ=255)
                                                                                                   aTGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTa                                                   <  2:541601/150‑1 (MQ=255)
                                                                                                    tGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTAc                                                  <  2:2294928/150‑1 (MQ=255)
                                                                                                     gTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACa                                                 >  1:1648230/1‑150 (MQ=255)
                                                                                                         gATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCtt                                             <  2:912297/150‑1 (MQ=255)
                                                                                                         gATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCtt                                             >  2:4265671/1‑150 (MQ=255)
                                                                                                           tataAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTc                                           >  1:2927093/1‑150 (MQ=255)
                                                                                                              aaaCCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGc                                        >  2:3707157/1‑150 (MQ=255)
                                                                                                              aaaCCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGAAATTTTCGTTCCGATGGTAACACCAGCGATAATCGTGGGTGCGGCCTTACAGCTTTCAGc                                        >  1:2864261/1‑150 (MQ=255)
                                                                                                                aCCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCACGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTg                                      <  2:2924275/150‑1 (MQ=255)
                                                                                                                 ccATGCAGGTTCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGa                                     <  2:4153639/150‑1 (MQ=255)
                                                                                                                  cATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGaa                                    >  2:2961273/1‑150 (MQ=255)
                                                                                                                  cATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGaa                                    >  2:2443504/1‑150 (MQ=255)
                                                                                                                    tGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAaac                                  >  2:2403788/1‑150 (MQ=255)
                                                                                                                    tGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAaac                                  >  2:4360904/1‑150 (MQ=255)
                                                                                                                        ggTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACt                              >  1:421100/1‑150 (MQ=255)
                                                                                                                          ttCTGACCTTATATTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGc                            >  1:658126/1‑150 (MQ=255)
                                                                                                                                  ttATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTg                    >  1:1094782/1‑150 (MQ=255)
                                                                                                                                   tataTTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGt                   >  2:2174817/1‑150 (MQ=255)
                                                                                                                                   tataTTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTGGGTGCCATTACCCATCTGTTGGGAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTGACGGCTTTCCGCTGTACCGACTGCCtctggtgt                   >  2:3276191/1‑143 (MQ=255)
                                                                                                                                    atatTTGAACTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGtt                  >  2:4633875/1‑150 (MQ=255)
                                                                                                                                     tatTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGttt                 <  1:2262394/150‑1 (MQ=255)
                                                                                                                                         tGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTAt             >  2:840977/1‑150 (MQ=255)
                                                                                                                                         tGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTAt             >  1:2785223/1‑150 (MQ=255)
                                                                                                                                            aCTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGc          >  1:1601181/1‑150 (MQ=255)
                                                                                                                                            aCTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGc          <  1:3145759/150‑1 (MQ=255)
                                                                                                                                             cTGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCg         <  1:3102652/150‑1 (MQ=255)
                                                                                                                                              tGGAGGTTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCga        <  1:2700496/150‑1 (MQ=255)
                                                                                                                                              tGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCga        >  1:124048/1‑150 (MQ=255)
                                                                                                                                               ggaggaTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCgat       >  1:1739975/1‑150 (MQ=255)
                                                                                                                                                gaggaTCACCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCgatg      <  2:3822797/150‑1 (MQ=255)
                                                                                                                                                gaggaTCACCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCgatg      <  2:1226786/150‑1 (MQ=255)
                                                                                                                                                 aggaTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCgatga     <  1:1890730/150‑1 (MQ=255)
                                                                                                                                                   gaTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCGATGATc   >  1:2491512/1‑150 (MQ=255)
                                                                                                                                                   gaTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCGATGATc   >  1:1122265/1‑150 (MQ=255)
                                                                                                                                                    aTCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCGATGATCg  <  1:2570743/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAACTGCCTATCTTGTTTCTATGGCGATGATCG  >  CP009273/3017611‑3017908

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: