Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,584,756 A→T 74.2% *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,7560AT74.2% 492.5 / 201.0 256*305K (TAA→AAA) yjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (41/25);  new base T (88/102);  total (129/127)
Fisher's exact test for biased strand distribution p-value = 3.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 257 total aligned reads displayed.

AGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATAT  >  CP009273/4584610‑4584905
                                                                                                                                                  |                                                                                                                                                     
aGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATg                                                                                                                                                    <  1:505556/150‑1 (MQ=255)
  gTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAg                                                                                                                                                  <  2:2421578/150‑1 (MQ=255)
  gTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAg                                                                                                                                                  >  1:2062349/1‑150 (MQ=255)
     aaaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCAt                                                                                                                                               >  2:3087183/1‑150 (MQ=255)
      aaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATa                                                                                                                                              <  1:969322/150‑1 (MQ=255)
       aGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATaa                                                                                                                                             >  1:233692/1‑150 (MQ=255)
        gTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAAc                                                                                                                                            >  1:263726/1‑150 (MQ=255)
             aGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCgtg                                                                                                                                       <  2:3794530/150‑1 (MQ=255)
              ggACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTCTTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtgt                                                                                                                                      >  2:3985972/1‑150 (MQ=255)
              ggACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtgt                                                                                                                                      <  2:2343903/150‑1 (MQ=255)
               gACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGtt                                                                                                                                     <  2:4346393/150‑1 (MQ=255)
                aCTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTCCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTtc                                                                                                                                    >  2:2487084/1‑150 (MQ=255)
                 cTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTtct                                                                                                                                   <  1:1742016/150‑1 (MQ=255)
                 cTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTtct                                                                                                                                   <  1:1742291/150‑1 (MQ=255)
                     actacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCAtt                                                                                                                               <  2:2567334/150‑1 (MQ=255)
                            cTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGAt                                                                                                                        <  1:3298769/150‑1 (MQ=255)
                             tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                       <  1:1475443/150‑1 (MQ=255)
                                caAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGaa                                                                                                                    <  1:2808247/150‑1 (MQ=255)
                                caAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGaa                                                                                                                    <  1:2808479/150‑1 (MQ=255)
                                caAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGaa                                                                                                                    >  1:149503/1‑150 (MQ=255)
                                caAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGaa                                                                                                                    <  1:573645/150‑1 (MQ=255)
                                 aaaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                   <  1:137866/150‑1 (MQ=255)
                                  aaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATc                                                                                                                  <  2:841116/150‑1 (MQ=255)
                                  aaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATc                                                                                                                  >  1:4100801/1‑150 (MQ=255)
                                  aaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATc                                                                                                                  >  1:4100779/1‑150 (MQ=255)
                                     aCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCatt                                                                                                               <  1:3781977/150‑1 (MQ=255)
                                        tttttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCattgat                                                                                                            >  2:4282586/1‑150 (MQ=255)
                                            tAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGtt                                                                                                        <  2:3436225/150‑1 (MQ=255)
                                             aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGttt                                                                                                       >  2:3974056/1‑150 (MQ=255)
                                             aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGttt                                                                                                       >  2:717452/1‑150 (MQ=255)
                                             aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGttt                                                                                                       <  1:2436474/150‑1 (MQ=255)
                                             aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGttt                                                                                                       <  1:2436311/150‑1 (MQ=255)
                                               tCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAg                                                                                                     >  1:1101721/1‑150 (MQ=255)
                                                cGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGc                                                                                                    >  1:3507074/1‑150 (MQ=255)
                                                 gACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCt                                                                                                   <  1:4137091/150‑1 (MQ=255)
                                                  aCTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTg                                                                                                  <  1:1007587/150‑1 (MQ=255)
                                                  aCTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTg                                                                                                  >  1:3036783/1‑150 (MQ=255)
                                                     gATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAAt                                                                                               <  2:358494/150‑1 (MQ=255)
                                                     gATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAAt                                                                                               <  2:359592/150‑1 (MQ=255)
                                                        tgtgAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCa                                                                                            <  1:2237532/150‑1 (MQ=255)
                                                        tgtgAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCa                                                                                            >  2:3061326/1‑150 (MQ=255)
                                                         gtgAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCaa                                                                                           <  1:1619180/150‑1 (MQ=255)
                                                           gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                         >  1:1227295/1‑150 (MQ=255)
                                                           gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                         >  1:1423526/1‑150 (MQ=255)
                                                           gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                         >  2:3493614/1‑150 (MQ=255)
                                                           gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                         >  2:3492650/1‑150 (MQ=255)
                                                           gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                         <  2:135958/150‑1 (MQ=255)
                                                             gCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATgg                                                                                       >  2:638274/1‑150 (MQ=255)
                                                               aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                     >  1:3341795/1‑150 (MQ=255)
                                                               aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTGTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                     >  2:2672936/1‑150 (MQ=255)
                                                                tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                    <  1:3174289/150‑1 (MQ=255)
                                                                tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                    <  2:2570710/150‑1 (MQ=255)
                                                                  ggTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACg                                                                                  <  1:3600391/150‑1 (MQ=255)
                                                                    tAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGa                                                                                >  2:487063/1‑150 (MQ=255)
                                                                             aTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAggg                                                                       >  1:370124/1‑150 (MQ=255)
                                                                               ttCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTg                                                                     >  2:2189427/1‑150 (MQ=255)
                                                                                  ccGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCa                                                                  >  1:2780085/1‑150 (MQ=255)
                                                                                    ggTATTTATCTTTTGATAAAACCGTGCACCCGCTAATATTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGTTCGAATCATTGATTGTTTAGCTGAAGTCATATGGTTACGGATCGCAAGGGTGGCATc                                                                <  1:3241406/150‑1 (MQ=255)
                                                                                    ggTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATc                                                                >  2:3271059/1‑150 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGTGCAAGGGTGGCATCaa                                                              >  1:4381369/1‑150 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                              >  1:2366639/1‑150 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                              >  1:2366687/1‑150 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                              >  1:2173634/1‑150 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                              <  2:233801/150‑1 (MQ=255)
                                                                                      tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                              >  2:167434/1‑150 (MQ=255)
                                                                                       aTTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAg                                                             <  1:694889/150‑1 (MQ=255)
                                                                                        tttATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAgg                                                            <  2:2965995/150‑1 (MQ=255)
                                                                                          tATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTc                                                          <  1:4532091/150‑1 (MQ=255)
                                                                                             cTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACt                                                       >  1:1452040/1‑150 (MQ=255)
                                                                                               tttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                     >  1:1233933/1‑150 (MQ=255)
                                                                                                ttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCg                                                    <  2:1129130/150‑1 (MQ=255)
                                                                                                   aTAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCa                                                 >  2:1827205/1‑150 (MQ=255)
                                                                                                     aaaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCaaa                                               >  2:2072171/1‑150 (MQ=255)
                                                                                                     aaaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCaaa                                               >  2:4316899/1‑150 (MQ=255)
                                                                                                            tGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGTCATCAAGGTCACTGCGGCAAATCAGAgc                                        >  2:4503534/1‑150 (MQ=255)
                                                                                                            tGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgc                                        >  1:4287150/1‑150 (MQ=255)
                                                                                                               aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                     <  1:1760952/150‑1 (MQ=255)
                                                                                                               aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                     <  1:2712099/150‑1 (MQ=255)
                                                                                                                  cGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgag                                  >  2:1893948/1‑150 (MQ=255)
                                                                                                                  cGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgag                                  <  1:4282993/150‑1 (MQ=255)
                                                                                                                   gCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgaga                                 >  1:4180478/1‑150 (MQ=255)
                                                                                                                   gCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgaga                                 >  2:2764763/1‑150 (MQ=255)
                                                                                                                   gCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGAGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgaga                                 >  2:2764499/1‑150 (MQ=255)
                                                                                                                      aaTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATg                              >  2:2871446/1‑150 (MQ=255)
                                                                                                                          gTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATTGTCa                          >  2:927219/1‑150 (MQ=255)
                                                                                                                          gTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                          >  2:4384740/1‑150 (MQ=255)
                                                                                                                             tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                       <  2:4189594/150‑1 (MQ=255)
                                                                                                                             tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                       <  2:4191043/150‑1 (MQ=255)
                                                                                                                                    tttAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATc                <  1:1605914/150‑1 (MQ=255)
                                                                                                                                     ttAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCg               >  2:2308929/1‑150 (MQ=255)
                                                                                                                                      tAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGa              >  1:587190/1‑150 (MQ=255)
                                                                                                                                        aTCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACg            <  2:1303013/150‑1 (MQ=255)
                                                                                                                                           cGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCa         <  2:1325416/150‑1 (MQ=255)
                                                                                                                                           cGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCa         <  2:3172806/150‑1 (MQ=255)
                                                                                                                                            gTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCaa        >  2:196590/1‑150 (MQ=255)
                                                                                                                                            gTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCaa        >  2:196554/1‑150 (MQ=255)
                                                                                                                                            gTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCaa        >  2:196604/1‑150 (MQ=255)
                                                                                                                                               cTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa     <  1:675003/150‑1 (MQ=255)
                                                                                                                                                 tAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGata   >  2:1649120/1‑150 (MQ=255)
                                                                                                                                                  aaTGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGatat  >  1:3545912/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATAT  >  CP009273/4584610‑4584905

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: