Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,982 G→T 100% T226T (ACG→ACT topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,9820GT98.6% 644.9 / 0.4 209T226T (ACG→ACTtopADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (106/97);  total (109/100)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01
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 209 total aligned reads displayed.

ATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325836‑1326127
                                                                                                                                                  |                                                                                                                                                 
atGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAc                                                                                                                                                 >  2:759120/1‑149 (MQ=255)
atGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAc                                                                                                                                                 <  2:2484632/149‑1 (MQ=255)
 tGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAcc                                                                                                                                                >  1:638513/1‑149 (MQ=255)
  ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCa                                                                                                                                               <  2:491994/149‑1 (MQ=255)
   gTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCAc                                                                                                                                              >  1:3752180/1‑149 (MQ=255)
       cgccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCa                                                                                                                                          <  2:906845/149‑1 (MQ=255)
         ccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                        <  2:576540/149‑1 (MQ=255)
          cgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCt                                                                                                                                       >  1:3257710/1‑149 (MQ=255)
          cgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCt                                                                                                                                       <  2:229896/149‑1 (MQ=255)
             tgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                    <  1:1936274/149‑1 (MQ=255)
             tgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                    <  1:1936504/149‑1 (MQ=255)
             tgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                    >  2:2685277/1‑149 (MQ=255)
              gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCGTCTGGTg                                                                                                                                   <  2:2573103/149‑1 (MQ=255)
              gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTg                                                                                                                                   >  1:2924801/1‑149 (MQ=255)
               ctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGa                                                                                                                                  <  2:2811264/149‑1 (MQ=255)
               ctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGa                                                                                                                                  >  2:2418230/1‑149 (MQ=255)
                    gAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGTGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGt                                                                                                                             >  1:3549896/1‑149 (MQ=255)
                    gAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGt                                                                                                                             <  1:3852755/149‑1 (MQ=255)
                       aaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACCCCATATGGTGAAGCGTTgg                                                                                                                          >  1:766342/1‑149 (MQ=255)
                         aGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCg                                                                                                                        >  2:828602/1‑149 (MQ=255)
                          gATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGt                                                                                                                       <  2:1776334/149‑1 (MQ=255)
                          gATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGt                                                                                                                       <  2:16445/149‑1 (MQ=255)
                              gCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACa                                                                                                                   >  2:2405581/1‑149 (MQ=255)
                               cTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAg                                                                                                                  <  1:3834188/149‑1 (MQ=255)
                               cTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAg                                                                                                                  <  2:2913507/149‑1 (MQ=255)
                                  gTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTg                                                                                                               <  2:1347085/149‑1 (MQ=255)
                                      ccTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtc                                                                                                           >  1:1266854/1‑149 (MQ=255)
                                      ccTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCGGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTAGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtc                                                                                                           >  1:3131886/1‑149 (MQ=255)
                                       cTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtca                                                                                                          <  2:1871781/149‑1 (MQ=255)
                                       cTGTCTGCCGGTCGTGTGCAGTCGGTGGCGATTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtca                                                                                                          <  2:3070015/149‑1 (MQ=255)
                                       cTGTCTGCCGGTCGTGTGCAGTCGGTGGCGATTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtca                                                                                                          <  2:39233/149‑1 (MQ=255)
                                         gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatc                                                                                                        >  1:4060902/1‑149 (MQ=255)
                                         gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatc                                                                                                        >  1:2516003/1‑149 (MQ=255)
                                             gCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGaa                                                                                                    >  2:35669/1‑149 (MQ=255)
                                             gCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCAGGTGGTCGAGCGTGAGCGTGAAATTAAATCGTTAGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAACACCCCAACTGGTGAAGCGTTGGCGTTACAGGTGACTCAACAGaa                                                                                                    >  1:2806265/1‑149 (MQ=255)
                                              ccGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAAc                                                                                                   <  1:555491/149‑1 (MQ=255)
                                                ggTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGa                                                                                                 <  2:1739521/149‑1 (MQ=255)
                                                     tgtgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAAc                                                                                            >  2:1720478/1‑149 (MQ=255)
                                                      gtgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAAcc                                                                                           >  1:149895/1‑149 (MQ=255)
                                                          aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtc                                                                                       <  2:2592198/149‑1 (MQ=255)
                                                          aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtc                                                                                       <  1:923078/149‑1 (MQ=255)
                                                           gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtcc                                                                                      >  2:2010366/1‑149 (MQ=255)
                                                           gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtcc                                                                                      >  1:2836014/1‑149 (MQ=255)
                                                              ggTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgtc                                                                                   <  2:3655239/149‑1 (MQ=255)
                                                               gTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgtcc                                                                                  <  1:2354288/149‑1 (MQ=255)
                                                                    ggTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTca                                                                             >  1:2932980/1‑149 (MQ=255)
                                                                     gTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            >  1:1030676/1‑149 (MQ=255)
                                                                     gTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            >  2:2502347/1‑149 (MQ=255)
                                                                     gTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            >  1:2214617/1‑149 (MQ=255)
                                                                      ttCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaac                                                                           >  1:1244629/1‑149 (MQ=255)
                                                                       tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaaca                                                                          >  1:3896586/1‑149 (MQ=255)
                                                                           cTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGa                                                                      <  1:2004076/149‑1 (MQ=255)
                                                                           cTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGa                                                                      <  2:1441992/149‑1 (MQ=255)
                                                                               tggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaa                                                                  >  1:1950035/1‑149 (MQ=255)
                                                                                ggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaaa                                                                 <  2:1296642/149‑1 (MQ=255)
                                                                                 gtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAAc                                                                <  2:4066189/149‑1 (MQ=255)
                                                                                    gagcgtgagcgtGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCa                                                             >  1:4009181/1‑149 (MQ=255)
                                                                                       cgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGc                                                          <  2:1075758/149‑1 (MQ=255)
                                                                                       cgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGc                                                          <  1:85850/149‑1 (MQ=255)
                                                                                        gtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCt                                                         <  2:3752180/149‑1 (MQ=255)
                                                                                        gtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCt                                                         >  1:442408/1‑149 (MQ=255)
                                                                                            gcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCgg                                                     >  2:1977383/1‑149 (MQ=255)
                                                                                              gtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTGATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTa                                                   >  1:2102240/1‑149 (MQ=255)
                                                                                                gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                 >  1:1249294/1‑149 (MQ=255)
                                                                                                 aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                >  2:3460097/1‑149 (MQ=255)
                                                                                                    ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  2:1387839/1‑149 (MQ=255)
                                                                                                    ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  2:3510279/1‑149 (MQ=255)
                                                                                                    ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  2:1503224/1‑149 (MQ=255)
                                                                                                     tAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgc                                            >  2:564668/1‑149 (MQ=255)
                                                                                                     tAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgc                                            <  1:3211010/149‑1 (MQ=255)
                                                                                                      aaaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgct                                           <  2:2533197/149‑1 (MQ=255)
                                                                                                      aaaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgct                                           <  1:1555459/149‑1 (MQ=255)
                                                                                                        aGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTgg                                         >  2:203516/1‑149 (MQ=255)
                                                                                                             tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaaa                                    >  1:3191636/1‑149 (MQ=255)
                                                                                                                    ggAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCt                             <  2:276450/149‑1 (MQ=255)
                                                                                                                     gAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                            <  1:740758/149‑1 (MQ=255)
                                                                                                                     gAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                            >  1:2275577/1‑149 (MQ=255)
                                                                                                                     gAAGAGTTCTGGGAAGTCGATGCCAGCACTAACAAGCCATCAGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGATCCGTCCGGTCAAAAAAGAACATACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                            >  1:1027413/1‑149 (MQ=255)
                                                                                                                       agagTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACa                          <  2:1542670/149‑1 (MQ=255)
                                                                                                                        gagTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAg                         >  2:3139270/1‑149 (MQ=255)
                                                                                                                          gTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCg                       <  1:3735464/149‑1 (MQ=255)
                                                                                                                           ttCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGt                      >  1:298949/1‑149 (MQ=255)
                                                                                                                           ttCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGt                      >  1:1690175/1‑149 (MQ=255)
                                                                                                                           ttCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGt                      >  2:3999616/1‑149 (MQ=255)
                                                                                                                           ttCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGt                      >  2:2328617/1‑149 (MQ=255)
                                                                                                                             cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGc                    <  1:3899285/149‑1 (MQ=255)
                                                                                                                              tGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCt                   >  1:3171862/1‑149 (MQ=255)
                                                                                                                                 gAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGa                >  2:3874389/1‑149 (MQ=255)
                                                                                                                                  aaGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa               >  1:2501372/1‑149 (MQ=255)
                                                                                                                                  aaGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa               <  2:715471/149‑1 (MQ=255)
                                                                                                                                    gTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACg             <  2:4058351/149‑1 (MQ=255)
                                                                                                                                     tCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGt            <  1:1073582/149‑1 (MQ=255)
                                                                                                                                       gATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGa          >  1:2726300/1‑149 (MQ=255)
                                                                                                                                        aTGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGaa         <  2:880449/149‑1 (MQ=255)
                                                                                                                                         tGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAg        >  2:495837/1‑149 (MQ=255)
                                                                                                                                           ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGAc      >  2:3218435/1‑149 (MQ=255)
                                                                                                                                           ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGAc      <  1:2418230/149‑1 (MQ=255)
                                                                                                                                           ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCGACAGCGTGCTGGAACGTGAAGAc      <  1:346835/149‑1 (MQ=255)
                                                                                                                                              gCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa   >  1:2809941/1‑149 (MQ=255)
                                                                                                                                              gCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa   >  2:1860977/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325836‑1326127

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: