Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,822,340 A→C 17.7% Q12H (CAA→CAC) ‡ xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,822,3400AC17.7% 300.8 / 24.5 198Q12H (CAA→CAC) ‡xanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (100/59);  new base C (6/33);  total (106/92)
Fisher's exact test for biased strand distribution p-value = 6.99e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTG  >  CP009273/3822193‑3822487
                                                                                                                                                    |                                                                                                                                                   
  aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAAcc                                                                                                                                                   >  1:901506/1‑148 (MQ=255)
    aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCgg                                                                                                                                                 <  1:3135955/148‑1 (MQ=255)
     aaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CACCCCGt                                                                                                                                                >  2:3026305/1‑148 (MQ=255)
      aaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGG‑CACCCGCtt                                                                                                                                               >  1:3355885/1‑148 (MQ=255)
         tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTgcg                                                                                                                                            >  2:1237966/1‑148 (MQ=255)
         tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAAACGGGTgcg                                                                                                                                            >  1:636623/1‑148 (MQ=255)
          gCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGG‑GCACCGGGTgcgc                                                                                                                                           >  1:2405190/1‑148 (MQ=255)
          gCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTgcgc                                                                                                                                           <  1:840695/148‑1 (MQ=255)
            aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAg                                                                                                                                         <  1:2098233/148‑1 (MQ=255)
              ttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGAc                                                                                                                                       >  2:1151748/1‑148 (MQ=255)
                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCGGTTGCGCCGACTCa                                                                                                                                    >  1:2247620/1‑148 (MQ=255)
                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCa                                                                                                                                    >  1:1138155/1‑148 (MQ=255)
                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCa                                                                                                                                    >  1:188999/1‑148 (MQ=255)
                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTCCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCa                                                                                                                                    >  1:2095304/1‑148 (MQ=255)
                  ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCCGTTGCGCAGACTCaa                                                                                                                                   >  2:1271605/1‑148 (MQ=255)
                   tcATTTAAATATTTAGATTAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGATAAACTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCaaa                                                                                                                                  <  2:2241173/148‑1 (MQ=255)
                   tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCCCTTGCGCAGACTCaaa                                                                                                                                  >  1:1259311/1‑148 (MQ=255)
                     ttttAAATATTTAGATGAGGTTGTGGCAAACGTTTGCTTTTCTATTTCGGTCAGATAAAATGCCCGTTTTGCTACCACGGTTTTTTTTTCGTTTTGTGTTTCCACCCTCGTGTCAGAAATTGCG‑CACCCGGTTGCGCAGACTCAAAAc                                                                                                                                <  1:564102/147‑1 (MQ=255)
                     aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAAc                                                                                                                                >  2:1757478/1‑148 (MQ=255)
                       ttAAATATTTAGATGTGGTTGTGGCAACCGTTTGCTTTTCTCGTTCGGTCAGCTAAATTGCCCGCTGTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCATAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAg                                                                                                                              <  2:3465719/148‑1 (MQ=255)
                         aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCg                                                                                                                            >  2:436700/1‑148 (MQ=255)
                         aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCg                                                                                                                            >  2:2705942/1‑148 (MQ=255)
                            tatTTAGATGTGTATGTGGCAAACTTTTGCTTTTCTCTTTCGTTCAGCTAAAATGCCCGCTTTGCTACCAGGGGATTTTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAATTGCG‑CAACGGGTGGCGCAGACTCAAAACAGCGAAc                                                                                                                         <  2:2/148‑1 (MQ=255)
                             atTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACAACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CACCCGGTTGCCCAGACTCAAAAAAGCCaaat                                                                                                                        >  2:3405605/1‑146 (MQ=255)
                             atTTAGATGAGGATGTGGCAAACGTTTGCTTGTCTCTTTCGGTCCGCTAAAAGGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACt                                                                                                                        <  2:2264277/148‑1 (MQ=255)
                                 aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCCAACAGAtt                                                                                                                    >  2:1219392/1‑148 (MQ=255)
                                   aTGAGGTTGTGGAAAAGGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGATACCAGGGTATTTTTTTCGTTGTCTTTTTCAACCCTCGTGTCAGAAATTGCG‑CAACGGTTTGCGCAGACTCAAAACAGCGAACTGATTTa                                                                                                                  <  2:606244/148‑1 (MQ=255)
                                   aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                                                                                                  <  2:3522458/148‑1 (MQ=255)
                                      aGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCAAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCAACCCG‑AGGGCCGAATCAAAAAACCGAACCCATTTAcca                                                                                                               >  1:3326529/1‑147 (MQ=255)
                                      aGGATGTGGCAAACGTTTGCTTTTCTCTGTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCg                                                                                                               <  2:3253047/148‑1 (MQ=255)
                                         aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGAGCAGACTCAAAACAGCGAAATGATTTAACGACt                                                                                                            >  2:3533854/1‑148 (MQ=255)
                                              gCAAAGGTTTGCTTTTCTCTTTCGGACAGCTAAACGGCCCGCTTTGCTACCTGGGGTTTTTTTTCGTTGTGTGTGCCCCCCCTCGAGTCAGAAATGGCC‑CACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                                                                                                       <  1:419606/148‑1 (MQ=255)
                                               cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                      <  2:3277233/148‑1 (MQ=255)
                                                 aaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                    >  2:546891/1‑148 (MQ=255)
                                                   cttttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCTCTTTGCTACCAGGGGATTGTTTTCGATGTCTGTTTCCACCCTCTAGTCAGAAAATGCG‑CACCCGGTTGCGCAGACTAAAAACAGCGAACTGATTTACCGTCTTGAAGATCGt                                                                                                  <  2:173543/146‑1 (MQ=255)
                                                   cGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGt                                                                                                  >  2:1196009/1‑148 (MQ=255)
                                                    gTTTGCTTTTCTCTTTCGGTCAGAAAAACTGCCCGCTTTGCTACCAGGGGTTTGTTTTCGATGTCTGTTACCACCCTCGAGACACAAAATGCG‑CCCCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                 <  1:2079342/148‑1 (MQ=255)
                                                     tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCCGGTGCGCCGACTCAAAAAATCCAAAAGATTTAACGTCTTGAAGATCGTcc                                                                                                >  1:117886/1‑148 (MQ=255)
                                                      ttGCTTTTCTCTTTCGGTCAGCTAAATTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                               <  2:3225725/148‑1 (MQ=255)
                                                      ttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGGCTTGAAGATCGTccg                                                                                               >  2:644165/1‑148 (MQ=255)
                                                      ttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCCAACAGATTTAACGTCTTGAATATAGTccg                                                                                               >  1:2951308/1‑148 (MQ=255)
                                                      ttGCTTTTCTCGTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                               <  2:2960273/148‑1 (MQ=255)
                                                        gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCCGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                                                                                             >  1:1452791/1‑148 (MQ=255)
                                                            ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCCCCGGTTGCGCAGACTCAAAACACCCAACAGATTTTACGTCTTGAAGATCGTCCGCCGCtt                                                                                         >  1:3040579/1‑148 (MQ=255)
                                                                ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGG‑CCACCCGTTGCGCAGACTCAAAACAGCGAACTGATTTACAGTCTTGAAGAACGTCCGCCGCGTCCTc                                                                                     >  1:681041/1‑148 (MQ=255)
                                                                 tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGGTGCGCAGACTCAAAACCGCGACCTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa                                                                                    >  2:1283954/1‑148 (MQ=255)
                                                                 tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATGGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGGTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa                                                                                    >  2:2678809/1‑148 (MQ=255)
                                                                 tgtCGGTCAGCTAAAATCCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAAACGGTTGCGCAGACTCAAAACAGAGAAATGATTTACCGAAATGTAGAAAGTCCGCCGCTTCCTCa                                                                                    >  2:2679119/3‑148 (MQ=255)
                                                                  ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaa                                                                                   <  2:874391/148‑1 (MQ=255)
                                                                     ggTCAGCTAAAATGCCCGCTGTGCTACCACGGGATTGTTTTCGATGTCTGTTGCCACCCTCGAGCCAGAAAAGGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAcc                                                                                <  2:2295829/148‑1 (MQ=255)
                                                                      gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCCGTTGCGCAGACTCAAAACAGCGAACAGATTTAACGTCTTGTATATCGTGCGCCGCGTCCTCAAAacc                                                                               >  1:1329685/1‑148 (MQ=255)
                                                                        cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCCCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg                                                                             >  2:1591545/1‑148 (MQ=255)
                                                                           cTAAACTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGttt                                                                          <  2:3355229/148‑1 (MQ=255)
                                                                            tAATATGCCCGTTTTGATACCTGGTGTTTGTTTTCGTTGTCTGTTACCACCCTCGATACATAAATGGCG‑CACGCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTg                                                                         <  1:1383524/148‑1 (MQ=255)
                                                                            tAAAATGCCCGTTTTGCTACCAGGGGATTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCATAAAATGCG‑CAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTg                                                                         <  2:1504634/148‑1 (MQ=255)
                                                                            tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTg                                                                         <  2:2657417/148‑1 (MQ=255)
                                                                               aaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGGTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACGCAAAACCGCGAACTGATTTACCGTCTTGTAGATCGTCCGCCGCTTCCTCACACCCTGTTTgccg                                                                      >  2:2140005/1‑148 (MQ=255)
                                                                               aaTGCCCGCTTTGCGACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCGGTTGCGCAGACTCAAAACACCGAACTGATTTAACGTCTTGTAGATCGTCCGCCGCGTCCTCAAACCCAGTTTggcg                                                                      >  1:2161044/1‑148 (MQ=255)
                                                                                 tGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgcc                                                                    <  2:547239/148‑1 (MQ=255)
                                                                                 tGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgcc                                                                    >  2:2461588/1‑148 (MQ=255)
                                                                                   cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTg                                                                  <  2:342437/148‑1 (MQ=255)
                                                                                    ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGt                                                                 >  2:260858/1‑148 (MQ=255)
                                                                                      gCTTTGCTACCCCGGGATTGTTGTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATCCG‑CAACCGGTTGCGCAGACTCAAACCAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                               <  2:1828630/148‑1 (MQ=255)
                                                                                        tttGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGc                                                             <  2:3186407/148‑1 (MQ=255)
                                                                                         ttGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CCACCCCTTGCGCACACTCAAAAAAACGAACAGATTTAACGTGTTGAATATCGTCCGCCGCGTCCTCAAACCCAGTTTGGCGCCGGTCAGCa                                                            >  1:723498/1‑148 (MQ=255)
                                                                                           gCTACCAGGTGTTTGTTTTCGTTGTCTTTTTCCCCCCTCGAGAAAGAAATTGCG‑CAACGGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                          <  1:203675/148‑1 (MQ=255)
                                                                                           gCTACCAGGGTATTGTTTTCGATGTCTTTTCCCACCCTCGAGTCATAAAAGGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                          <  1:1199903/148‑1 (MQ=255)
                                                                                           gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                          <  2:2095000/148‑1 (MQ=255)
                                                                                           gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                          <  2:2095339/148‑1 (MQ=255)
                                                                                             tACCACGGTTTTGTTTTCGTTGTCTGTTCCCCCCCCCGAGCCACAAAAGGCG‑CACCCGGTCGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctg                                                        <  1:352885/148‑1 (MQ=255)
                                                                                             tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctg                                                        <  2:3041940/148‑1 (MQ=255)
                                                                                                  cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCCGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCCTCCGCCGCTTCCTCAAAACCAGGTTGCCGCCTGTCAGCAACTGCTGGc                                                   >  1:750350/1‑148 (MQ=255)
                                                                                                   gggTTTGTTTTCTATGTCTGTTTCCACCCTCGAGTAAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCg                                                  <  2:723498/148‑1 (MQ=255)
                                                                                                       ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGt                                              <  2:719966/148‑1 (MQ=255)
                                                                                                        tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACGCTGTGTGCCGCCGGTCAGCATCTGCTGGCGAGGtt                                             >  2:1441089/1‑148 (MQ=255)
                                                                                                        tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGtt                                             >  2:2355/1‑148 (MQ=255)
                                                                                                          ttttAGTTGTCTTTTTCCACCCTCGGGTAAAAAAATGCG‑CACCCGTGGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCg                                           <  1:537530/148‑1 (MQ=255)
                                                                                                               gATGTCTGTTTCCACCCTCGAGCCAGAAAATGCG‑CCCCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                      <  1:2647065/148‑1 (MQ=255)
                                                                                                                   tgtgtgTCCACCCTCGAGTCATAAATTGCC‑AAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                  <  1:3149003/146‑1 (MQ=255)
                                                                                                                   tCTGTTCCCCCCCTCGAGTCACAAAAGGCG‑CACCGGGTGGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                  <  1:286415/148‑1 (MQ=255)
                                                                                                                     tGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                <  2:1705575/148‑1 (MQ=255)
                                                                                                                     tGTTTCCACCCTCGACTCACAAAATGCC‑CACCCGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                <  1:580156/148‑1 (MQ=255)
                                                                                                                       tttCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCGGTTTTCCGCCTGTCAGCATCAGCTGGCGCTGTTCGTTGCGGTGATAAc                              >  1:3264755/1‑148 (MQ=255)
                                                                                                                        ttCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGAAGCAAAAAAGCGAACTGATTTAACGTCTTGAAGATCGTGCGCCGCGTCCTCAAAACCTGTTTTCCGCCTGTCAGCATCTGGTGGCGATGTTCGTTGCGGTGATAACg                             >  1:2912599/1‑148 (MQ=255)
                                                                                                                            aCCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                         >  2:1654201/1‑148 (MQ=255)
                                                                                                                              ccTCGAGACAGAAAAGGCG‑CAACCGGTGGCGCAGACTCAAACCAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                       <  2:587720/148‑1 (MQ=255)
                                                                                                                                tCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                     >  2:414503/1‑148 (MQ=255)
                                                                                                                                tCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                     <  2:1550598/148‑1 (MQ=255)
                                                                                                                                tCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAACCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                     >  2:1971280/1‑148 (MQ=255)
                                                                                                                                 cGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                    <  2:1793286/148‑1 (MQ=255)
                                                                                                                                    gacataaacTGCG‑CCCCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTa                 <  1:1933202/143‑1 (MQ=255)
                                                                                                                                     tCAGAAAATGCG‑CACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTaa                <  2:87985/148‑1 (MQ=255)
                                                                                                                                       aGAAACTGCG‑CACCCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc              <  1:2430906/148‑1 (MQ=255)
                                                                                                                                        gAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCGTCCTCAAAACCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATAACGCCACCGCTATTAATCt             >  1:1157935/1‑148 (MQ=255)
                                                                                                                                         aaacTGCG‑CACCCGGGTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg            <  1:49011/148‑1 (MQ=255)
                                                                                                                                            aTGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa         >  2:146561/1‑148 (MQ=255)
                                                                                                                                            aTGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa         >  2:2010664/1‑148 (MQ=255)
                                                                                                                                             tGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAg        >  2:2964028/1‑148 (MQ=255)
                                                                                                                                                    aCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg  >  2:1243762/1‑148 (MQ=255)
                                                                                                                                                    aCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg  >  1:2719134/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCG‑CAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTG  >  CP009273/3822193‑3822487

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