Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,624,581 C→A 5.6% E183* (GAA→TAA)  rob ← right oriC‑binding transcriptional activator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,624,5810CA 5.6% 1080.5 / 25.8 284E183* (GAA→TAA) robright oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base C (112/156);  new base A (10/6);  total (122/162)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 313 total aligned reads displayed.

CAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCA  >  CP009273/4624432‑4624730
                                                                                                                                                     |                                                                                                                                                     
cAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttc                                                                                                                                                       >  2:1155897/1‑150 (MQ=255)
 aGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttca                                                                                                                                                      >  1:3448399/1‑150 (MQ=255)
  gACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcat                                                                                                                                                     <  1:1173149/150‑1 (MQ=255)
        cATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGc                                                                                                                                               <  1:930454/150‑1 (MQ=255)
         aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGcc                                                                                                                                              <  1:3848791/150‑1 (MQ=255)
         aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGcc                                                                                                                                              >  2:1424815/1‑150 (MQ=255)
         aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGACAGTACATAGCCATCTGCCTGAACCTGGGCTAACGCGGTGGTAAAGAAAACCTCATGCTCGTCGTCTTTAACCTGACACGGACGCGATTCATTCAGGcc                                                                                                                                              >  2:2290850/1‑150 (MQ=255)
           aGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGt                                                                                                                                            >  2:3202942/1‑150 (MQ=255)
            ggTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTa                                                                                                                                           >  1:3090730/1‑150 (MQ=255)
            ggTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTa                                                                                                                                           >  2:1458355/1‑150 (MQ=255)
             gTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag                                                                                                                                          >  1:3105105/1‑150 (MQ=255)
               aaaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag                                                                                                                                        <  2:2386365/150‑1 (MQ=255)
               aaaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag                                                                                                                                        <  1:193971/150‑1 (MQ=255)
                 aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCa                                                                                                                                      <  2:1580165/150‑1 (MQ=255)
                 aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCa                                                                                                                                      <  2:4182907/150‑1 (MQ=255)
                    tCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcg                                                                                                                                   <  1:4299495/150‑1 (MQ=255)
                    tCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcg                                                                                                                                   >  2:4371337/1‑150 (MQ=255)
                     cacaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcgg                                                                                                                                  <  1:718480/150‑1 (MQ=255)
                      acaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggc                                                                                                                                 >  2:1858304/1‑150 (MQ=255)
                      acaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggc                                                                                                                                 >  1:2970461/1‑150 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGTCGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          <  1:751382/150‑1 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          >  2:2806785/1‑150 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          >  2:3272456/1‑150 (MQ=255)
                                 gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCgg                                                                                                                      <  2:2748246/150‑1 (MQ=255)
                                       cagcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGtt                                                                                                                <  2:1571/150‑1 (MQ=255)
                                       cagcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGtt                                                                                                                <  2:1727372/150‑1 (MQ=255)
                                       cagcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGtt                                                                                                                <  2:3785141/150‑1 (MQ=255)
                                       cagcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGtt                                                                                                                <  2:1413751/150‑1 (MQ=255)
                                             cACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCgag                                                                                                          >  2:971387/1‑150 (MQ=255)
                                               ccGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGaa                                                                                                        <  2:3896432/150‑1 (MQ=255)
                                               ccGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGaa                                                                                                        <  2:489854/150‑1 (MQ=255)
                                               ccGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGaa                                                                                                        <  2:491904/150‑1 (MQ=255)
                                                 ggATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGaaaa                                                                                                      >  2:2201300/1‑150 (MQ=255)
                                                   aTGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATc                                                                                                    <  1:1496275/150‑1 (MQ=255)
                                                     gCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGTCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGt                                                                                                  <  1:3695711/150‑1 (MQ=255)
                                                     gCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGt                                                                                                  >  2:695920/1‑150 (MQ=255)
                                                         cGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCa                                                                                              >  1:3570870/1‑150 (MQ=255)
                                                          gTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAg                                                                                             <  2:524906/150‑1 (MQ=255)
                                                          gTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAg                                                                                             <  2:3548755/150‑1 (MQ=255)
                                                          gTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAg                                                                                             <  2:1277088/150‑1 (MQ=255)
                                                             aGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAAc                                                                                          >  2:3486091/1‑150 (MQ=255)
                                                              gTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACt                                                                                         <  2:2777696/150‑1 (MQ=255)
                                                               tACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTg                                                                                        <  1:2534316/150‑1 (MQ=255)
                                                                 cATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGAt                                                                                      >  1:3291316/1‑150 (MQ=255)
                                                                    aGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAAc                                                                                   <  1:3402336/150‑1 (MQ=255)
                                                                    aGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAAc                                                                                   >  1:173859/1‑150 (MQ=255)
                                                                      ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGc                                                                                 >  2:417297/1‑150 (MQ=255)
                                                                        aTCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCAt                                                                               >  2:1504408/1‑150 (MQ=255)
                                                                             ccTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCAt                                                                          <  1:822518/150‑1 (MQ=255)
                                                                              cTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATg                                                                         >  1:3402433/1‑150 (MQ=255)
                                                                               tGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATgg                                                                        <  1:3530940/150‑1 (MQ=255)
                                                                               tGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATgg                                                                        >  1:3181368/1‑150 (MQ=255)
                                                                                gATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGc                                                                       <  1:1082282/150‑1 (MQ=255)
                                                                                gATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGc                                                                       >  1:1398256/1‑150 (MQ=255)
                                                                                  tCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCgg                                                                     >  1:1203321/1‑150 (MQ=255)
                                                                                    cTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGaa                                                                   >  1:2227954/1‑150 (MQ=255)
                                                                                    cTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGaa                                                                   >  1:405023/1‑150 (MQ=255)
                                                                                          tAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCaga                                                             <  1:36250/150‑1 (MQ=255)
                                                                                           aaCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCagag                                                            <  1:3774176/150‑1 (MQ=255)
                                                                                              gcgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAtt                                                         <  1:1623832/150‑1 (MQ=255)
                                                                                              gcgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAtt                                                         <  1:2063964/150‑1 (MQ=255)
                                                                                                      ataGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAg                                                 <  2:2197485/150‑1 (MQ=255)
                                                                                                        aGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCg                                               <  1:2855931/150‑1 (MQ=255)
                                                                                                        aGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCg                                               <  2:2454548/150‑1 (MQ=255)
                                                                                                            tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACa                                           >  2:770554/1‑150 (MQ=255)
                                                                                                            tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACa                                           <  2:2813834/150‑1 (MQ=255)
                                                                                                              ccTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAgg                                         >  2:1540746/1‑150 (MQ=255)
                                                                                                               ctctTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGa                                        <  2:1657686/150‑1 (MQ=255)
                                                                                                               ctctTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGa                                        <  1:458175/150‑1 (MQ=255)
                                                                                                                  ttGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTa                                     >  1:2481930/1‑150 (MQ=255)
                                                                                                                   tGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAg                                    >  2:3925861/1‑150 (MQ=255)
                                                                                                                   tGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAg                                    <  1:2449949/150‑1 (MQ=255)
                                                                                                                   tGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAg                                    <  1:2649849/150‑1 (MQ=255)
                                                                                                                     cTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGct                                  <  2:2374002/150‑1 (MQ=255)
                                                                                                                     cTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGct                                  >  2:4440886/1‑150 (MQ=255)
                                                                                                                        gtcgtcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTg                               <  1:495823/150‑1 (MQ=255)
                                                                                                                        gtcgtcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTg                               <  2:310040/150‑1 (MQ=255)
                                                                                                                          cgtcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTggg                             >  1:4180321/1‑150 (MQ=255)
                                                                                                                           gtcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGt                            >  1:2120563/1‑150 (MQ=255)
                                                                                                                            tcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTCGCTCTGGGTa                           >  1:3991795/1‑150 (MQ=255)
                                                                                                                            tcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTa                           <  1:1138695/150‑1 (MQ=255)
                                                                                                                            tcTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTa                           <  1:1138785/150‑1 (MQ=255)
                                                                                                                               taaTCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCAGGGTAaca                        >  2:2650202/3‑150 (MQ=255)
                                                                                                                                 aTCCTGACTCGGACGCGTTTCATTCAGGCCTTAGAGCAACGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACAcc                      <  1:4371780/150‑1 (MQ=255)
                                                                                                                                 aTCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACAcc                      <  1:2215391/150‑1 (MQ=255)
                                                                                                                                 aTCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACAcc                      >  1:3014851/1‑150 (MQ=255)
                                                                                                                                   ccTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCaa                    >  2:927769/1‑150 (MQ=255)
                                                                                                                                    cTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAt                   <  1:2290706/150‑1 (MQ=255)
                                                                                                                                     tGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATc                  <  1:3081960/150‑1 (MQ=255)
                                                                                                                                         tCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcg              >  2:4395664/1‑150 (MQ=255)
                                                                                                                                           ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggc            >  2:283706/1‑150 (MQ=255)
                                                                                                                                           ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggc            <  1:1838222/150‑1 (MQ=255)
                                                                                                                                           ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggc            <  2:2069801/150‑1 (MQ=255)
                                                                                                                                            gACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggcg           >  2:2866881/1‑150 (MQ=255)
                                                                                                                                                 gTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCt      <  2:3108522/150‑1 (MQ=255)
                                                                                                                                                    tcattcaGGCCGTAGAGCACCGGCGGATTGGTCGGCGCTTTGCCGAGTAAATCGTGCCAGATCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTcc   <  1:1424815/150‑1 (MQ=255)
                                                                                                                                                    tcattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTcc   <  1:4333489/150‑1 (MQ=255)
                                                                                                                                                    tcattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTcc   <  1:2051543/150‑1 (MQ=255)
                                                                                                                                                     cattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCa  <  2:3422765/150‑1 (MQ=255)
                                                                                                                                                     cattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCa  <  2:2962182/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCA  >  CP009273/4624432‑4624730

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: