Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,584,756 A→T 93.3% *305K (TAA→AAA)  yjjM ← yjjMN operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,584,7560AT93.3% 796.5 / 36.1 238*305K (TAA→AAA) yjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base T (118/104);  total (126/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 238 total aligned reads displayed.

CAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGA  >  CP009273/4584609‑4584902
                                                                                                                                                   |                                                                                                                                                  
  ggTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGa                                                                                                                                                >  1:1117592/1‑150 (MQ=255)
     tAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCa                                                                                                                                             >  2:4529903/1‑150 (MQ=255)
      aaaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCAt                                                                                                                                            >  1:1827768/1‑150 (MQ=255)
      aaaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCAt                                                                                                                                            >  1:1788076/1‑150 (MQ=255)
       aaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATa                                                                                                                                           <  2:4342023/150‑1 (MQ=255)
        aGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATaa                                                                                                                                          >  2:1263837/1‑150 (MQ=255)
         gTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAc                                                                                                                                         >  1:4435436/1‑150 (MQ=255)
            ttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAcgcg                                                                                                                                      <  1:338838/150‑1 (MQ=255)
             tAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgt                                                                                                                                     >  2:595314/1‑150 (MQ=255)
             tAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgt                                                                                                                                     <  1:3103445/150‑1 (MQ=255)
                gACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGtt                                                                                                                                  >  2:49447/1‑150 (MQ=255)
                       ctacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTa                                                                                                                           >  1:4500425/1‑150 (MQ=255)
                         acCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAAt                                                                                                                         <  2:1189171/150‑1 (MQ=255)
                         acCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAAt                                                                                                                         <  1:4402995/150‑1 (MQ=255)
                          ccGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATg                                                                                                                        >  1:2534352/1‑150 (MQ=255)
                              tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTTAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                    <  1:4267494/150‑1 (MQ=255)
                              tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                    <  1:3964472/150‑1 (MQ=255)
                              tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                    <  2:3788781/150‑1 (MQ=255)
                              tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                    <  2:1977551/150‑1 (MQ=255)
                                  aaaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                >  1:714549/1‑150 (MQ=255)
                                  aaaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                >  2:2904086/1‑150 (MQ=255)
                                    aaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCa                                                                                                              >  2:4210221/1‑150 (MQ=255)
                                     aaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCat                                                                                                             >  1:833459/1‑150 (MQ=255)
                                     aaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCat                                                                                                             >  1:1025035/1‑150 (MQ=255)
                                         tttttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCattgat                                                                                                         >  1:3508772/1‑150 (MQ=255)
                                                tCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAg                                                                                                  >  1:2592604/1‑150 (MQ=255)
                                                 cGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGc                                                                                                 >  2:2073165/1‑150 (MQ=255)
                                                 cGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGc                                                                                                 >  2:3352894/1‑150 (MQ=255)
                                                 cGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGc                                                                                                 >  1:4053905/1‑150 (MQ=255)
                                                     tGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGaa                                                                                             <  2:4309538/150‑1 (MQ=255)
                                                           tgAGCATAGGTAAGCGGTAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCaaa                                                                                       <  1:3107944/150‑1 (MQ=255)
                                                           tgAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCaaa                                                                                       >  2:2437087/1‑150 (MQ=255)
                                                            gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                      <  2:4226451/150‑1 (MQ=255)
                                                             aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                     <  1:473348/150‑1 (MQ=255)
                                                              gCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATgg                                                                                    <  2:2901370/150‑1 (MQ=255)
                                                               cATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGt                                                                                   >  1:1445472/1‑150 (MQ=255)
                                                                aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                  <  1:2005099/150‑1 (MQ=255)
                                                                aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                  >  2:1338260/1‑150 (MQ=255)
                                                                 tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                 <  1:3804531/150‑1 (MQ=255)
                                                                 tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                 <  1:3804235/150‑1 (MQ=255)
                                                                     tAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGa                                                                             >  1:3599348/1‑150 (MQ=255)
                                                                       aGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTGCGGAgc                                                                           >  2:4490135/1‑150 (MQ=255)
                                                                        gCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAgcg                                                                          <  1:851021/150‑1 (MQ=255)
                                                                            gAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAg                                                                      >  1:2806035/1‑150 (MQ=255)
                                                                             aaTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCTAgg                                                                     >  2:3726805/1‑150 (MQ=255)
                                                                             aaTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAgg                                                                     >  2:4064782/1‑150 (MQ=255)
                                                                              aTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAggg                                                                    <  2:3534006/150‑1 (MQ=255)
                                                                               tttCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGt                                                                   >  1:2736327/1‑150 (MQ=255)
                                                                                        aTTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAg                                                          <  2:814808/150‑1 (MQ=255)
                                                                                           tATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTc                                                       <  2:63556/150‑1 (MQ=255)
                                                                                           tATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTc                                                       <  2:4404800/150‑1 (MQ=255)
                                                                                            aTCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCa                                                      <  2:3287161/150‑1 (MQ=255)
                                                                                               ttttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTg                                                   >  2:3904072/1‑150 (MQ=255)
                                                                                               ttttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTg                                                   >  2:2865136/1‑150 (MQ=255)
                                                                                                tttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                  >  1:2785015/1‑150 (MQ=255)
                                                                                                  tGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCgg                                                >  2:2866110/1‑150 (MQ=255)
                                                                                                   gATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGc                                               >  2:1900182/1‑150 (MQ=255)
                                                                                                   gATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGc                                               >  2:1077863/1‑150 (MQ=255)
                                                                                                    aTAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCa                                              <  1:233793/150‑1 (MQ=255)
                                                                                                     tAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCaa                                             >  2:1446627/1‑150 (MQ=255)
                                                                                                      aaaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCaaa                                            <  1:1349326/150‑1 (MQ=255)
                                                                                                       aaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                           <  2:2237466/150‑1 (MQ=255)
                                                                                                        aaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATc                                          <  2:3499536/150‑1 (MQ=255)
                                                                                                           cGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCaga                                       <  1:3650693/150‑1 (MQ=255)
                                                                                                             tGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgc                                     >  2:1808677/1‑150 (MQ=255)
                                                                                                             tGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgc                                     <  1:864027/150‑1 (MQ=255)
                                                                                                                aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                  <  1:2728144/150‑1 (MQ=255)
                                                                                                                aCCCGCTAATAGTTTGCTGTTATTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                  >  1:1927197/1‑150 (MQ=255)
                                                                                                                 cccGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTg                                 >  2:2646234/1‑150 (MQ=255)
                                                                                                                    gCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgaga                              <  1:944879/150‑1 (MQ=255)
                                                                                                                    gCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgaga                              >  1:4469723/1‑150 (MQ=255)
                                                                                                                       aaTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATg                           <  1:100056/150‑1 (MQ=255)
                                                                                                                       aaTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATg                           >  1:2134885/1‑150 (MQ=255)
                                                                                                                        aTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATgg                          >  1:3157109/1‑150 (MQ=255)
                                                                                                                         tAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGt                         <  2:1735489/150‑1 (MQ=255)
                                                                                                                          aGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTc                        >  2:2195850/1‑150 (MQ=255)
                                                                                                                           gTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                       <  1:3177717/150‑1 (MQ=255)
                                                                                                                           gTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                       >  2:3729661/1‑150 (MQ=255)
                                                                                                                            tttGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCat                      <  2:3756274/150‑1 (MQ=255)
                                                                                                                            tttGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCat                      <  1:2411461/150‑1 (MQ=255)
                                                                                                                             ttGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                     <  1:1548047/150‑1 (MQ=255)
                                                                                                                              tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                    <  2:867544/150‑1 (MQ=255)
                                                                                                                              tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                    <  2:2517788/150‑1 (MQ=255)
                                                                                                                                cTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGc                  >  2:1466876/1‑150 (MQ=255)
                                                                                                                                cTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGc                  >  1:3844135/1‑150 (MQ=255)
                                                                                                                                  gTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                >  2:2524278/1‑150 (MQ=255)
                                                                                                                                  gTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                >  2:1098486/1‑150 (MQ=255)
                                                                                                                                   tttttAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCa               >  2:3294790/1‑150 (MQ=255)
                                                                                                                                   tttttAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCa               <  2:4611192/150‑1 (MQ=255)
                                                                                                                                    ttttAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCAt              <  1:703314/150‑1 (MQ=255)
                                                                                                                                       tAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGa           <  2:2295201/150‑1 (MQ=255)
                                                                                                                                       tAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGa           <  2:1922595/150‑1 (MQ=255)
                                                                                                                                        aaTCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGAc          <  2:2400372/150‑1 (MQ=255)
                                                                                                                                        aaTCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGAc          <  2:1001548/150‑1 (MQ=255)
                                                                                                                                         aTCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACg         >  1:1505901/1‑150 (MQ=255)
                                                                                                                                          tCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACgg        >  2:451973/1‑150 (MQ=255)
                                                                                                                                            cGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCa      <  2:589367/150‑1 (MQ=255)
                                                                                                                                              tACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAAt    <  1:2002069/150‑1 (MQ=255)
                                                                                                                                              tACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAAt    >  2:1025369/1‑150 (MQ=255)
                                                                                                                                                cTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa  <  1:1570713/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGA  >  CP009273/4584609‑4584902

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: