Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,509,102 G→T 100% S57R (AGC→AGA dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,509,1020GT99.2% 466.9 / ‑4.1 147S57R (AGC→AGAdamDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (84/61);  total (85/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 148 total aligned reads displayed.

GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAC  >  CP009273/3508957‑3509246
                                                                                                                                                 |                                                                                                                                                 
ggaactgaTCCTGGCTTTTTTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTctct                                                                                                                                                <  2:632789/148‑1 (MQ=255)
 gaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTa                                                                                                                                               >  2:2497465/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCTGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              >  1:2593292/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              >  1:172570/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              >  2:3576884/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              >  1:2144417/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              >  1:2056149/1‑148 (MQ=255)
  aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                              <  1:1500237/148‑1 (MQ=255)
      gaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGat                                                                                                                                          <  1:2698519/148‑1 (MQ=255)
       aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCCGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                         >  1:2087094/1‑148 (MQ=255)
        tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                        <  2:2274245/148‑1 (MQ=255)
        tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                        >  2:3257306/1‑148 (MQ=255)
        tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                        <  2:2274041/148‑1 (MQ=255)
              cTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCaa                                                                                                                                  >  2:3182790/1‑148 (MQ=255)
                tttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAgg                                                                                                                                <  2:894545/148‑1 (MQ=255)
                 ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGa                                                                                                                               >  2:2188008/1‑148 (MQ=255)
                  tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGat                                                                                                                              <  2:899599/148‑1 (MQ=255)
                  tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGat                                                                                                                              <  2:2630271/148‑1 (MQ=255)
                      gAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATaa                                                                                                                          <  1:2134521/148‑1 (MQ=255)
                        aCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACg                                                                                                                        >  1:650725/1‑148 (MQ=255)
                          tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgag                                                                                                                      <  2:307237/148‑1 (MQ=255)
                          tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgag                                                                                                                      >  1:453955/1‑148 (MQ=255)
                             tCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGaaa                                                                                                                   <  2:1602841/148‑1 (MQ=255)
                               gcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                                                 <  1:3747795/148‑1 (MQ=255)
                                cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                                                <  1:3059899/148‑1 (MQ=255)
                                cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                                                <  1:2905834/148‑1 (MQ=255)
                                  gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                                                              >  2:3307655/1‑148 (MQ=255)
                                  gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                                                              >  1:3568240/1‑148 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             >  2:1962429/1‑148 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             >  2:1223257/1‑148 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             <  2:2680608/148‑1 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             <  2:759113/148‑1 (MQ=255)
                                   aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGGTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                             >  1:1490211/1‑148 (MQ=255)
                                       gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGtt                                                                                                         >  1:2762372/1‑148 (MQ=255)
                                        atatTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCTGACTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTg                                                                                                        <  1:269136/148‑1 (MQ=255)
                                              aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAac                                                                                                  <  2:848524/148‑1 (MQ=255)
                                              aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAac                                                                                                  <  2:1356555/148‑1 (MQ=255)
                                                      cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                          >  1:527132/1‑148 (MQ=255)
                                                      cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                          >  1:1129912/1‑148 (MQ=255)
                                                          aTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                      >  1:1307528/1‑148 (MQ=255)
                                                                ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGCCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACa                                                                                >  1:356711/1‑148 (MQ=255)
                                                                    ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                            >  1:1893650/1‑148 (MQ=255)
                                                                    ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                            >  2:1501101/1‑148 (MQ=255)
                                                                     gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                           >  2:2321523/1‑148 (MQ=255)
                                                                        caaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                        >  1:111516/1‑148 (MQ=255)
                                                                        caaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                        >  2:454434/1‑148 (MQ=255)
                                                                        caaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                        >  2:1649289/1‑148 (MQ=255)
                                                                          aacaGCTCGCGTGCGGCCTGTACGTACACATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCaa                                                                      <  2:3739619/148‑1 (MQ=255)
                                                                           acaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAc                                                                     >  2:1117441/1‑148 (MQ=255)
                                                                                 cgcgTGCGGCATTTACGTACTCATCAGTACGCATATACACAATGTTATAGAGACTGATCAGGTATCTTTTTATATCGCCAAGGAAATAACGAGAAAAGTCGGTTTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACa                                                               <  1:1595580/148‑1 (MQ=255)
                                                                                  gcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                              <  2:3447051/148‑1 (MQ=255)
                                                                                       gggCCTTTACTTACTAATCATTACTAATGTTAAAATTTTTATAGACACTGATCTGGTATATATTTATATCGGAACGGATATAACTACAAAAGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGcc                                                         <  2:3520168/147‑1 (MQ=255)
                                                                                       cGGCCTGTACGTACTCATCAGTACGCATCTCCACAATGTTATAGAGACTGATCAGGGCTCTATTGATATCGGCAAGGATATAACGAGAAATGTCGGTGTTGAGAAACACCGACCCGGCACCTACAACAGGCTCAACCAGACATTCGcc                                                         <  2:2656404/148‑1 (MQ=255)
                                                                                        ggCCTTTACGTACTCATCAGTACGCATCTACACAATGTTATAGAGACTTATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGccc                                                        <  2:1810154/148‑1 (MQ=255)
                                                                                        ggCCTGTACGTACTCAACAGTACGCATCTTCACAATGTTATAGAGACTGATCCGGTCTCTATTGATATCGGCACGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGccc                                                        <  2:3700511/148‑1 (MQ=255)
                                                                                         gCCTGTACGTAATATTCAGTACGCATCTTCACAATGTTATAGAGACTTATCTGGTCTCTTTTGATATCGGCATGTATATAACGATAAAAGTCGGTTTTGAGAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCt                                                       <  2:650725/148‑1 (MQ=255)
                                                                                           cTGTACGTACTCATCTGTACGAATCTTCAAAATGTTTTAGATCCTGATCAGGTCTATTTTTATATCGGCAAGGATATACCGAGAAAAGTGGGTGTTGAGAAAAACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                     <  2:974611/148‑1 (MQ=255)
                                                                                           cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                     <  2:1864043/148‑1 (MQ=255)
                                                                                           cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                     <  2:1864083/148‑1 (MQ=255)
                                                                                            tGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTgg                                                    >  1:1604376/1‑148 (MQ=255)
                                                                                             gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                   >  1:3519722/1‑148 (MQ=255)
                                                                                                cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaa                                                <  2:1543003/148‑1 (MQ=255)
                                                                                                  tacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAAt                                              >  1:1016182/1‑148 (MQ=255)
                                                                                                    cTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                            >  2:224583/1‑148 (MQ=255)
                                                                                                     tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTTTTGATATCGGCAAGGATATAACGACAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGcc                                           <  2:1085742/148‑1 (MQ=255)
                                                                                                       atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGt                                         >  2:1327892/1‑148 (MQ=255)
                                                                                                         caGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                       >  1:1356122/1‑148 (MQ=255)
                                                                                                          aGTACGAATCTACACATTGTTATAGAGACTGATCAGGTCTCTTTTTATATCGCCAAGTATATAACGACAAAAGTGGGTGTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTa                                      <  1:2769645/148‑1 (MQ=255)
                                                                                                           gaccGCATCTTCACAATGTTATAGACACTGATCCGGTCTCTTTTGATATCGCCCAGTATATAACGATAAATGTCGGTGTTCACAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa                                     <  1:3290567/145‑1 (MQ=255)
                                                                                                            tACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                    >  1:936889/1‑148 (MQ=255)
                                                                                                            tACGAATCTTCACATTTTTATAGAGACTGTTCCGGTATCTTTTTTTTTCGGAATGTATATAACTATAATAGTGGGTTTTGAGAACCCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                    <  1:2287375/148‑1 (MQ=255)
                                                                                                                 aTCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                               <  2:2450763/148‑1 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACTAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             >  2:1820937/1‑148 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             >  2:3711429/1‑148 (MQ=255)
                                                                                                                   cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                             >  1:2989286/1‑148 (MQ=255)
                                                                                                                      cacaATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                          <  2:2056149/148‑1 (MQ=255)
                                                                                                                       acaATTTTATAGATACTGATCAGGTCTTTATTGATATCGGCATGTATATAACGATAAAAGTCGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCa                         <  1:2850109/148‑1 (MQ=255)
                                                                                                                       acaATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCa                         >  2:1712847/1‑148 (MQ=255)
                                                                                                                        caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                        >  2:742650/1‑148 (MQ=255)
                                                                                                                        caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                        >  1:2369764/1‑148 (MQ=255)
                                                                                                                         ccggTTCTAGAGACTGAACCGTCGTCGATTGATATGGGAAAGGATTTAAGGAGAAAACACCGTGTTGAGAAAGACCGACCCGGCACCTACCAATCGCTCAACGATAGATTCGCCATTGGGCAAATGCCTTTTAATATCATCAAGCAgg                       <  2:2030564/145‑1 (MQ=255)
                                                                                                                         attGTTATTGTGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                       <  2:1515707/146‑1 (MQ=255)
                                                                                                                          aTGTTATAGCGACTGATCAGGTCTATATTGATATCGGCATGTATATAACGACAAACGTCGGTTTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                      <  1:2235622/148‑1 (MQ=255)
                                                                                                                          aTGTTATAGAGACTGATCTGGTCTCTTTTTATATCGGCAAGCATATAACGAGAAAAGTCGGTGTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                      <  1:2235863/148‑1 (MQ=255)
                                                                                                                          aTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                      >  1:396379/1‑148 (MQ=255)
                                                                                                                            gTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGa                    <  2:399947/148‑1 (MQ=255)
                                                                                                                              tataGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                  >  1:2562887/1‑148 (MQ=255)
                                                                                                                                  gagaCGTATCTGGTCTCTATGGATAGGGCCAAGGAAAAAACTACAATAGTGGGTGTTGAGAAAACCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTg              <  1:2752224/148‑1 (MQ=255)
                                                                                                                                   agaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc             >  2:1964730/1‑148 (MQ=255)
                                                                                                                                   agaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTGTAATATCATCAAGCAGGGGATCCTTGc             >  2:2918660/1‑148 (MQ=255)
                                                                                                                                     aCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc           >  1:2442799/1‑148 (MQ=255)
                                                                                                                                     aCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc           >  2:1571461/1‑148 (MQ=255)
                                                                                                                                      cTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcccc          <  2:3686003/148‑1 (MQ=255)
                                                                                                                                      cTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTGGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcccc          >  2:3102568/1‑148 (MQ=255)
                                                                                                                                         aTCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg       >  2:252080/1‑148 (MQ=255)
                                                                                                                                           cAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc     >  1:3786391/1‑148 (MQ=255)
                                                                                                                                           cAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc     >  2:343149/1‑148 (MQ=255)
                                                                                                                                            aGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGGCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc    >  2:1231477/1‑148 (MQ=255)
                                                                                                                                            aGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc    >  1:2500753/1‑148 (MQ=255)
                                                                                                                                            aGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc    >  2:1229187/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAC  >  CP009273/3508957‑3509246

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