Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 837,698 G→A 15.8% intergenic (‑186/‑90) mcbA ← / → rlmF colanic acid mucoidy stimulation protein/23S rRNA m(6)A1618 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273837,6980GA15.8% 634.1 / 89.8 227intergenic (‑186/‑90)mcbA/rlmFcolanic acid mucoidy stimulation protein/23S rRNA m(6)A1618 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (118/73);  new base A (24/12);  total (142/85)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 227 total aligned reads displayed.

GTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATC  >  CP009273/837552‑837845
                                                                                                                                                  |                                                                                                                                                   
gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTa                                                                                                                                                  <  2:2761557/150‑1 (MQ=255)
gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTa                                                                                                                                                  <  2:2291145/150‑1 (MQ=255)
gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTa                                                                                                                                                  <  1:2003576/150‑1 (MQ=255)
  gtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAg                                                                                                                                                <  2:896201/150‑1 (MQ=255)
    aGGTCGCCTGCGCGATGCAAGTGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGt                                                                                                                                              <  2:4291517/150‑1 (MQ=255)
    aGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGt                                                                                                                                              >  1:153889/1‑150 (MQ=255)
       tCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGtttt                                                                                                                                           >  1:391740/1‑150 (MQ=255)
       tCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGtttt                                                                                                                                           >  1:2166160/1‑150 (MQ=255)
        cGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTg                                                                                                                                          <  2:145548/150‑1 (MQ=255)
              cgcgATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTATTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAAt                                                                                                                                    <  1:4436391/150‑1 (MQ=255)
                     cAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAAc                                                                                                                             >  2:1803574/1‑150 (MQ=255)
                      aaaTGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTTGGGAATGATAAACg                                                                                                                            >  1:2848758/1‑150 (MQ=255)
                      aaaTGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACg                                                                                                                            >  2:3892354/1‑150 (MQ=255)
                      aaaTGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACg                                                                                                                            <  1:2816387/150‑1 (MQ=255)
                         tGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAg                                                                                                                         >  2:1471833/1‑150 (MQ=255)
                         tGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAg                                                                                                                         >  2:1378226/1‑150 (MQ=255)
                           aCGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGcc                                                                                                                       <  2:4121536/150‑1 (MQ=255)
                              cAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGa                                                                                                                    <  1:2334469/150‑1 (MQ=255)
                                 aaaaTTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGatat                                                                                                                 >  1:2567496/1‑150 (MQ=255)
                                   aaTTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATAttt                                                                                                               >  2:1626300/1‑150 (MQ=255)
                                        aCAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCCTGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCt                                                                                                          >  2:3782491/1‑150 (MQ=255)
                                        aCAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCt                                                                                                          >  2:3783767/1‑150 (MQ=255)
                                          agagATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGt                                                                                                        >  2:2465250/1‑150 (MQ=255)
                                          agagATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGt                                                                                                        >  1:1230135/1‑150 (MQ=255)
                                          agagATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGt                                                                                                        >  1:4207486/1‑150 (MQ=255)
                                          agagATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGt                                                                                                        >  2:868227/1‑150 (MQ=255)
                                           gagaTGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAATCGCAGCCAGATATTTACCCTGTc                                                                                                       <  2:2648894/150‑1 (MQ=255)
                                           gagaTGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTc                                                                                                       >  2:3810744/1‑150 (MQ=255)
                                              aTGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTcgcg                                                                                                    >  2:4390692/1‑150 (MQ=255)
                                               tgtgCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGt                                                                                                   >  1:2360886/1‑150 (MQ=255)
                                                  gCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTAt                                                                                                <  2:1829140/150‑1 (MQ=255)
                                                     tATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCg                                                                                             >  1:2545704/1‑150 (MQ=255)
                                                     tATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCg                                                                                             >  1:3643093/1‑150 (MQ=255)
                                                      aTTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCgg                                                                                            >  2:4084305/1‑150 (MQ=255)
                                                      aTTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCgg                                                                                            <  1:635550/150‑1 (MQ=255)
                                                      aTTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCgg                                                                                            >  1:676014/1‑150 (MQ=255)
                                                      aTTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCgg                                                                                            >  2:3474727/1‑150 (MQ=255)
                                                       ttttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCggg                                                                                           <  2:3684849/150‑1 (MQ=255)
                                                       ttttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCggg                                                                                           <  1:3251807/150‑1 (MQ=255)
                                                           ttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATc                                                                                       >  2:4075890/1‑150 (MQ=255)
                                                            tGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCt                                                                                      >  1:615762/1‑150 (MQ=255)
                                                             gAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCtt                                                                                     <  2:1262987/150‑1 (MQ=255)
                                                              aaCGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTa                                                                                    <  1:3603785/150‑1 (MQ=255)
                                                               aCGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAc                                                                                   >  2:1767146/1‑150 (MQ=255)
                                                                     tatTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCgct                                                                             >  2:71588/1‑150 (MQ=255)
                                                                          ttACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATTATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTa                                                                        <  2:1878481/150‑1 (MQ=255)
                                                                               aaaTAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGga                                                                   >  1:1368134/1‑150 (MQ=255)
                                                                                 aTAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGgaga                                                                 >  1:4139505/1‑150 (MQ=255)
                                                                                  tAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGaa                                                                >  1:875817/1‑150 (MQ=255)
                                                                                  tAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGaa                                                                >  1:2377459/1‑150 (MQ=255)
                                                                                   aaGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAAt                                                               >  1:712059/1‑150 (MQ=255)
                                                                                   aaGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAAt                                                               >  1:2151000/1‑150 (MQ=255)
                                                                                   aaGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAAt                                                               >  1:2411584/1‑150 (MQ=255)
                                                                                      aCAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCt                                                            >  1:1541613/1‑150 (MQ=255)
                                                                                       cAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTa                                                           >  1:451686/1‑150 (MQ=255)
                                                                                        aactCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTAt                                                          <  1:2520239/147‑1 (MQ=255)
                                                                                        aaaTCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTAt                                                          >  1:1688529/1‑150 (MQ=255)
                                                                                             aaTGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCg                                                     <  1:3444902/150‑1 (MQ=255)
                                                                                              aTGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGc                                                    >  1:4637790/1‑150 (MQ=255)
                                                                                                gACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGccc                                                  <  1:3041876/150‑1 (MQ=255)
                                                                                                 aCATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCa                                                 >  1:3358097/1‑150 (MQ=255)
                                                                                                 aCATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCa                                                 >  2:3244111/1‑150 (MQ=255)
                                                                                                   aTGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGa                                               >  1:2639188/1‑150 (MQ=255)
                                                                                                      cATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAc                                            >  1:3630893/1‑150 (MQ=255)
                                                                                                       aTAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAcc                                           >  1:1477240/1‑150 (MQ=255)
                                                                                                       aTAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAcc                                           >  1:4170474/1‑150 (MQ=255)
                                                                                                         aaTGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCgg                                         >  2:3163651/1‑150 (MQ=255)
                                                                                                           tGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCgggg                                       >  2:1128522/1‑150 (MQ=255)
                                                                                                              tGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTg                                    <  1:4390617/150‑1 (MQ=255)
                                                                                                               gATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGc                                   <  2:4467079/150‑1 (MQ=255)
                                                                                                               gATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGc                                   >  2:4594278/1‑150 (MQ=255)
                                                                                                                aTTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCa                                  >  2:227182/1‑150 (MQ=255)
                                                                                                                aTTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTACa                                  >  1:681755/1‑150 (MQ=255)
                                                                                                                  ttttATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATc                                <  1:1892491/150‑1 (MQ=255)
                                                                                                                     tATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCcgc                             >  2:363068/1‑150 (MQ=255)
                                                                                                                      atgatgAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCcgcg                            <  2:2714228/150‑1 (MQ=255)
                                                                                                                       tgatgaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCgcgc                           >  1:1543562/1‑150 (MQ=255)
                                                                                                                       tgatgaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCgcgc                           >  2:817419/1‑150 (MQ=255)
                                                                                                                        gatgaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCa                          >  2:3098424/1‑150 (MQ=255)
                                                                                                                         atgaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCaa                         >  2:738343/1‑150 (MQ=255)
                                                                                                                          tgaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAAc                        >  1:3080042/1‑150 (MQ=255)
                                                                                                                           gaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAAcc                       <  2:279859/150‑1 (MQ=255)
                                                                                                                           gaATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAAcc                       <  1:3749445/150‑1 (MQ=255)
                                                                                                                            aaTTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCg                      >  2:4419288/1‑150 (MQ=255)
                                                                                                                            aaTTGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCg                      >  2:1422871/1‑150 (MQ=255)
                                                                                                                              ttGCCATAAACCATAGCCATATTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGtc                    >  1:4467517/1‑150 (MQ=255)
                                                                                                                               tGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGtca                   >  1:2658632/1‑150 (MQ=255)
                                                                                                                                 ccATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGtcatc                 <  2:1004447/150‑1 (MQ=255)
                                                                                                                                      aaCCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGc            >  2:412719/1‑150 (MQ=255)
                                                                                                                                      aaCCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGc            <  2:2132108/150‑1 (MQ=255)
                                                                                                                                       aCCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGcc           <  2:4182332/150‑1 (MQ=255)
                                                                                                                                       aCCATAGCCATATTACGGTTTTGGGGAATGATAAAAGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGAGTTGCATCCGCGCAACCGTCATCACAGcc           >  2:286985/1‑150 (MQ=255)
                                                                                                                                        ccATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCg          >  1:1023374/1‑150 (MQ=255)
                                                                                                                                          aTAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCt        >  2:4030995/1‑150 (MQ=255)
                                                                                                                                          aTAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCt        >  1:1739905/1‑150 (MQ=255)
                                                                                                                                          aTAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCt        <  1:2971666/150‑1 (MQ=255)
                                                                                                                                           tAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTa       <  1:1337271/150‑1 (MQ=255)
                                                                                                                                            aGCCATGTTAAGGTTTTGGGGAATTATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTAc      <  1:4096628/150‑1 (MQ=255)
                                                                                                                                               cATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGAt   <  2:660770/150‑1 (MQ=255)
                                                                                                                                               cATGGTAAGGTTTTGGGGAATGATAAACGCAGCCCGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCGGCTTAAAGGAGAATGCTATGTCCGCCCAGACACCGGGGTTGCATCCGCGCAACCGTCATCGCAGCCGCTACGAt   >  2:2943649/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATC  >  CP009273/837552‑837845

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: