Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,144,158 C→A 39.2% intergenic (+253/+19) yghA → / ← exbD putative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,144,1580CA39.2% 207.6 / 230.4 194intergenic (+253/+19)yghA/exbDputative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base C (69/49);  new base A (45/31);  total (114/80)
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.99e-01

Only 100 of 199 total aligned reads displayed.

ATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGT  >  CP009273/3144009‑3144305
                                                                                                                                                     |                                                                                                                                                   
aTGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAgc                                                                                                                                                     <  1:3064793/150‑1 (MQ=255)
 tGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAgcg                                                                                                                                                    <  1:4863139/150‑1 (MQ=255)
    tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCAt                                                                                                                                                 >  1:511380/1‑150 (MQ=255)
    tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCAt                                                                                                                                                 <  1:2732769/150‑1 (MQ=255)
      cactAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCa                                                                                                                                               <  2:4811439/147‑1 (MQ=255)
       aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAg                                                                                                                                              >  2:2016899/1‑150 (MQ=255)
       aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAg                                                                                                                                              >  2:3869946/1‑150 (MQ=255)
       aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAg                                                                                                                                              >  1:3830798/1‑150 (MQ=255)
       aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAg                                                                                                                                              >  1:3829701/1‑150 (MQ=255)
         tAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATCGATTTTGTAGGCCTGTTGAGCGGAGCGCATCAGGc                                                                                                                                            >  1:4952342/1‑150 (MQ=255)
         tAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGc                                                                                                                                            >  1:4628337/1‑150 (MQ=255)
         tAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGc                                                                                                                                            <  2:441786/150‑1 (MQ=255)
           gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCaa                                                                                                                                          >  2:2770956/1‑150 (MQ=255)
                  tAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtact                                                                                                                                   >  2:1073752/1‑150 (MQ=255)
                  tAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtact                                                                                                                                   >  2:713620/1‑150 (MQ=255)
                   aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactt                                                                                                                                  <  1:2448523/150‑1 (MQ=255)
                       cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCtactttact                                                                                                                              >  2:2206761/1‑150 (MQ=255)
                       cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCtactttact                                                                                                                              >  2:4887504/1‑150 (MQ=255)
                       cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCACTTCtactttact                                                                                                                              >  1:149800/1‑150 (MQ=255)
                       cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttact                                                                                                                              >  1:1316406/1‑150 (MQ=255)
                            aaGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCt                                                                                                                         >  1:947906/1‑150 (MQ=255)
                            aaGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCt                                                                                                                         >  2:3721197/1‑150 (MQ=255)
                             aGCGTAGCGCATCAGTCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGACTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                        <  2:3269707/150‑1 (MQ=255)
                             aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                        >  2:3684508/1‑150 (MQ=255)
                             aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                        >  2:1404756/1‑150 (MQ=255)
                               cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACTTGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTg                                                                                                                      <  1:1529797/150‑1 (MQ=255)
                                 tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGc                                                                                                                    >  2:650193/1‑150 (MQ=255)
                                    cgcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGt                                                                                                                 <  2:4711596/150‑1 (MQ=255)
                                     gcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGtt                                                                                                                >  2:23675/1‑150 (MQ=255)
                                      cATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                               <  1:3490986/150‑1 (MQ=255)
                                      cATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                               <  1:3738433/150‑1 (MQ=255)
                                      cATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                               <  1:3735169/150‑1 (MQ=255)
                                       aTCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttc                                                                                                              >  1:507121/1‑150 (MQ=255)
                                       aTCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttc                                                                                                              >  2:3117318/1‑150 (MQ=255)
                                       aTCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttc                                                                                                              <  1:805966/150‑1 (MQ=255)
                                        tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttct                                                                                                             <  2:4221700/150‑1 (MQ=255)
                                         cAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttctt                                                                                                            <  2:3898771/150‑1 (MQ=255)
                                          aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttcttc                                                                                                           <  2:5061551/150‑1 (MQ=255)
                                            gCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGc                                                                                                         <  2:3308017/150‑1 (MQ=255)
                                              aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                       >  1:3485798/1‑150 (MQ=255)
                                                 gAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAcc                                                                                                    <  2:3287529/150‑1 (MQ=255)
                                                     aCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAc                                                                                                >  1:1014435/1‑150 (MQ=255)
                                                       cTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                              <  1:550880/150‑1 (MQ=255)
                                                       cTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                              >  2:3684891/1‑150 (MQ=255)
                                                         aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAt                                                                                            >  1:3783788/1‑150 (MQ=255)
                                                         aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAt                                                                                            >  1:2232974/1‑150 (MQ=255)
                                                             ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTc                                                                                        <  2:4123587/150‑1 (MQ=255)
                                                              cGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCa                                                                                       >  1:2873096/1‑150 (MQ=255)
                                                               ggATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                      >  1:2010700/1‑150 (MQ=255)
                                                               ggATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                      >  1:2010851/1‑150 (MQ=255)
                                                                 aTGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGt                                                                                    <  1:3081668/150‑1 (MQ=255)
                                                                    cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                                 >  1:565320/1‑150 (MQ=255)
                                                                        gtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGc                                                                             >  2:2202167/1‑150 (MQ=255)
                                                                          gAACGCCTTATCCGGCCTACAAACGCTCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCGGCTGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCt                                                                           <  2:3813685/150‑1 (MQ=255)
                                                                              gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                       >  1:453129/1‑150 (MQ=255)
                                                                              gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                       >  2:1009510/1‑150 (MQ=255)
                                                                                    tCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGt                                                                 >  1:4192132/1‑150 (MQ=255)
                                                                                       ggCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATc                                                              >  2:4787458/1‑150 (MQ=255)
                                                                                           tACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAtt                                                          >  1:612103/1‑150 (MQ=255)
                                                                                           tACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTATTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAtt                                                          >  2:4678830/1‑150 (MQ=255)
                                                                                            aCAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTa                                                         >  1:1195375/1‑150 (MQ=255)
                                                                                                aCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCtt                                                     <  2:1809310/150‑1 (MQ=255)
                                                                                                 cGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtc                                                    >  2:3979411/1‑150 (MQ=255)
                                                                                                 cGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtc                                                    >  2:4264322/1‑150 (MQ=255)
                                                                                                   cACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcat                                                  >  2:3498828/1‑150 (MQ=255)
                                                                                                    aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                                 <  1:3665408/150‑1 (MQ=255)
                                                                                                      gAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCaa                                               >  1:1719757/1‑150 (MQ=255)
                                                                                                        aGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACg                                             >  2:2238765/1‑150 (MQ=255)
                                                                                                          aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                           >  2:4028957/1‑150 (MQ=255)
                                                                                                            tCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtctc                                         <  1:3558939/150‑1 (MQ=255)
                                                                                                             cAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCg                                        >  2:1225750/1‑150 (MQ=255)
                                                                                                                cATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTaa                                     <  1:3879113/150‑1 (MQ=255)
                                                                                                                cATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTaa                                     <  1:1392720/150‑1 (MQ=255)
                                                                                                                   gTTGCAATAGATTGTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                  <  2:263849/150‑1 (MQ=255)
                                                                                                                    ttGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGa                                 >  1:1096924/1‑150 (MQ=255)
                                                                                                                      gCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                               <  1:3032456/150‑1 (MQ=255)
                                                                                                                      gCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                               <  1:4296503/150‑1 (MQ=255)
                                                                                                                      gCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                               >  1:3990719/1‑150 (MQ=255)
                                                                                                                       cAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACgg                              <  1:4003328/150‑1 (MQ=255)
                                                                                                                           aGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                          >  1:113127/1‑150 (MQ=255)
                                                                                                                           aGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                          <  1:3583814/150‑1 (MQ=255)
                                                                                                                            gATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTg                         <  2:712334/150‑1 (MQ=255)
                                                                                                                              ttttGTAGGCCTGATGAGCGTAGCGCATCAGGCAATACTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTc                       >  2:742937/1‑150 (MQ=255)
                                                                                                                               tttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTcc                      <  2:3829701/150‑1 (MQ=255)
                                                                                                                               tttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTcc                      <  2:3830798/150‑1 (MQ=255)
                                                                                                                                   tAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTc                  >  2:2716890/1‑150 (MQ=255)
                                                                                                                                   tAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTc                  >  1:1293292/1‑150 (MQ=255)
                                                                                                                                    aGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCg                 >  2:2864445/1‑150 (MQ=255)
                                                                                                                                     ggCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCgg                >  2:4128091/1‑150 (MQ=255)
                                                                                                                                      gCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCGTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGga               <  1:752618/150‑1 (MQ=255)
                                                                                                                                       ccTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaa              <  2:3075812/150‑1 (MQ=255)
                                                                                                                                         tgatgaGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAGGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaaga            >  2:4135404/1‑150 (MQ=255)
                                                                                                                                          gatgaGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaagaa           >  1:775719/1‑150 (MQ=255)
                                                                                                                                            tgaGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAagaaga         >  1:2700668/1‑150 (MQ=255)
                                                                                                                                             gagcggagcgCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAt        <  1:4735260/150‑1 (MQ=255)
                                                                                                                                             gaGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAt        <  2:3533711/150‑1 (MQ=255)
                                                                                                                                             gaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAt        >  1:1915939/1‑150 (MQ=255)
                                                                                                                                              aGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtg       >  2:1342377/1‑150 (MQ=255)
                                                                                                                                                cGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGACGAAGAtggt     >  2:3949527/1‑150 (MQ=255)
                                                                                                                                                   aGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATggtggt  <  1:2206761/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGT  >  CP009273/3144009‑3144305

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: