Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,624,581 C→A 17.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,5810CA17.6% 863.5 / 124.7 289E183* (GAA→TAA) robright oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base C (125/113);  new base A (28/23);  total (153/136)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

CAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCA  >  CP009273/4624432‑4624730
                                                                                                                                                     |                                                                                                                                                     
   aCCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcatt                                                                                                                                                    <  1:3141782/150‑1 (MQ=255)
    ccTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcattc                                                                                                                                                   >  2:2859396/1‑150 (MQ=255)
    cattCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATTGTATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcattc                                                                                                                                                   <  2:3019957/148‑1 (MQ=255)
       tCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAgg                                                                                                                                                >  2:4170819/1‑150 (MQ=255)
         aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGcc                                                                                                                                              <  1:636927/150‑1 (MQ=255)
         aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGcc                                                                                                                                              >  2:3456898/1‑150 (MQ=255)
             gTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag                                                                                                                                          <  1:2025047/150‑1 (MQ=255)
               aaaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTGTCCTGACTCGGACGCGTTTCATTCAGGCCGTagag                                                                                                                                        >  2:544114/1‑150 (MQ=255)
                aaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGc                                                                                                                                       >  1:1117647/1‑150 (MQ=255)
                 aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCa                                                                                                                                      >  1:3360654/1‑150 (MQ=255)
                 aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCa                                                                                                                                      <  1:2765789/150‑1 (MQ=255)
                 aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCa                                                                                                                                      >  2:1087491/1‑150 (MQ=255)
                  cATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCAc                                                                                                                                     <  2:4838779/150‑1 (MQ=255)
                   aTCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCAcc                                                                                                                                    <  1:2006833/150‑1 (MQ=255)
                    tCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcg                                                                                                                                   <  1:1560919/150‑1 (MQ=255)
                     cacaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcgg                                                                                                                                  >  2:2368650/1‑150 (MQ=255)
                      acaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggc                                                                                                                                 >  1:3774914/1‑150 (MQ=255)
                       caTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACcggcg                                                                                                                                >  1:10741/1‑150 (MQ=255)
                        atatTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACcggcgg                                                                                                                               >  1:3605477/1‑150 (MQ=255)
                         tatTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGa                                                                                                                              <  1:2559494/150‑1 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          >  1:326246/1‑150 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          <  1:3433995/150‑1 (MQ=255)
                             cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATgg                                                                                                                          <  1:3641653/150‑1 (MQ=255)
                              gccgccCTGCAGCATCACCGGCTGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGt                                                                                                                         <  1:441025/150‑1 (MQ=255)
                              gccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGt                                                                                                                         >  1:1006502/1‑150 (MQ=255)
                               ccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTc                                                                                                                        <  2:547286/150‑1 (MQ=255)
                               ccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTc                                                                                                                        >  2:1484299/1‑150 (MQ=255)
                               ccgacCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTATCTCTGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCTGACGCGTTTCATTCATGCCGTAGTGCACCGTCGGAATTGTc                                                                                                                        <  2:3258592/150‑1 (MQ=255)
                                 gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCgg                                                                                                                      >  2:1857858/1‑150 (MQ=255)
                                 gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATTGTCgg                                                                                                                      >  1:800084/1‑150 (MQ=255)
                                  cccTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGgc                                                                                                                     <  2:731711/150‑1 (MQ=255)
                                        agcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTg                                                                                                               >  1:4699431/1‑150 (MQ=255)
                                        agcaTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTg                                                                                                               <  2:1066712/150‑1 (MQ=255)
                                          caTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGcc                                                                                                             <  2:1554343/150‑1 (MQ=255)
                                          caTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGcc                                                                                                             <  1:93661/150‑1 (MQ=255)
                                           aTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCg                                                                                                            <  2:2992209/150‑1 (MQ=255)
                                           aTCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCg                                                                                                            >  1:1040350/1‑150 (MQ=255)
                                            tCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCga                                                                                                           >  1:4991512/1‑150 (MQ=255)
                                              aCCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCTTATTTATCCTGACTCGGACGCGTTTCATTAAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCgaga                                                                                                         >  2:829737/1‑150 (MQ=255)
                                                cGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGaaa                                                                                                       <  2:2725968/150‑1 (MQ=255)
                                                   aTGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATc                                                                                                    <  1:743219/150‑1 (MQ=255)
                                                   aTGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATc                                                                                                    >  2:2171991/1‑150 (MQ=255)
                                                    tGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCAGTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCg                                                                                                   >  2:3870163/1‑150 (MQ=255)
                                                     gCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGt                                                                                                  >  2:4158909/1‑150 (MQ=255)
                                                      ccccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTg                                                                                                 >  1:3305011/1‑150 (MQ=255)
                                                       cccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGc                                                                                                >  2:3110480/1‑150 (MQ=255)
                                                        ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGcc                                                                                               >  2:3439586/1‑150 (MQ=255)
                                                        ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGcc                                                                                               >  2:2077031/1‑150 (MQ=255)
                                                              gTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACt                                                                                         <  2:101106/150‑1 (MQ=255)
                                                               tACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTg                                                                                        <  2:1537544/150‑1 (MQ=255)
                                                                aCATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGa                                                                                       <  2:42108/150‑1 (MQ=255)
                                                                aCATAGCCATCTGACTGATCCTGGGATAACGCGGTGGTTTAGAATACATCTTGCTCGTAGTCTTTATCCTGACTCTGACGAGTTTCATTCTGGCCGTAGAGCACCGGCGGAATTGTCGGCGCGTTGCCGTGAAAATCGTGCCAGAACTGa                                                                                       <  1:1740112/150‑1 (MQ=255)
                                                                 cATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGAt                                                                                      <  1:4735321/150‑1 (MQ=255)
                                                                   tAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATaa                                                                                    <  2:2248777/150‑1 (MQ=255)
                                                                    aGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAAc                                                                                   <  1:1330464/150‑1 (MQ=255)
                                                                    aGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAAc                                                                                   <  1:593711/150‑1 (MQ=255)
                                                                    aGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAAc                                                                                   <  1:2965694/150‑1 (MQ=255)
                                                                         tCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCAtt                                                                              >  2:2120299/1‑150 (MQ=255)
                                                                         tCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCAtt                                                                              <  1:3062000/150‑1 (MQ=255)
                                                                           tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTc                                                                            <  2:4974056/150‑1 (MQ=255)
                                                                            gCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCa                                                                           <  2:4730730/150‑1 (MQ=255)
                                                                             ccTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCAt                                                                          <  2:1659958/150‑1 (MQ=255)
                                                                              cTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATg                                                                         <  2:3977402/150‑1 (MQ=255)
                                                                               tGATCCTGGGCTAACGCGGTGGTATAGAAGACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATgg                                                                        <  2:2282355/150‑1 (MQ=255)
                                                                                gATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGc                                                                       >  2:361106/1‑150 (MQ=255)
                                                                                     tGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACAGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGaaa                                                                  >  2:180547/1‑150 (MQ=255)
                                                                                       ggCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAACTc                                                                >  2:4566594/1‑150 (MQ=255)
                                                                                       ggCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATc                                                                <  2:4417196/150‑1 (MQ=255)
                                                                                          tAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCaga                                                             <  1:4874571/150‑1 (MQ=255)
                                                                                           aaCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCagag                                                            <  2:591522/150‑1 (MQ=255)
                                                                                            aCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAgaga                                                           <  2:2630506/150‑1 (MQ=255)
                                                                                             cgcgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAt                                                          <  1:715858/150‑1 (MQ=255)
                                                                                              gcgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAtt                                                         <  2:2566916/150‑1 (MQ=255)
                                                                                               cgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAttt                                                        >  1:4300131/1‑150 (MQ=255)
                                                                                               cgGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGAttt                                                        <  1:408599/150‑1 (MQ=255)
                                                                                                  tggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCt                                                     >  1:4570541/1‑150 (MQ=255)
                                                                                                     tataGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCa                                                  >  2:4969902/1‑150 (MQ=255)
                                                                                                           aTACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAAc                                            <  1:664897/150‑1 (MQ=255)
                                                                                                            tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACa                                           <  2:5105332/150‑1 (MQ=255)
                                                                                                             aCCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGACTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAg                                          >  2:610433/1‑150 (MQ=255)
                                                                                                                tctTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAg                                       <  2:4860712/150‑1 (MQ=255)
                                                                                                                 ctTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGt                                      <  1:1891109/150‑1 (MQ=255)
                                                                                                                    gCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGc                                   <  1:3858037/150‑1 (MQ=255)
                                                                                                                     cTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGct                                  <  2:3635331/150‑1 (MQ=255)
                                                                                                                      tcgtcgTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGctc                                 >  1:2287260/1‑150 (MQ=255)
                                                                                                                      tcgtcgTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGctc                                 >  2:1848625/1‑150 (MQ=255)
                                                                                                                           gtcTTTATCCTGACTCGGACGCGTTTAATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGt                            >  2:4941062/1‑150 (MQ=255)
                                                                                                                             cTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTaa                          >  2:4367796/1‑150 (MQ=255)
                                                                                                                               ttATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAaca                        >  1:4526633/1‑150 (MQ=255)
                                                                                                                                       aCTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAg                <  1:748788/150‑1 (MQ=255)
                                                                                                                                        cTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgc               >  2:355755/1‑150 (MQ=255)
                                                                                                                                           ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggc            <  1:3528848/150‑1 (MQ=255)
                                                                                                                                            gACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggcg           >  2:1292300/1‑150 (MQ=255)
                                                                                                                                            gACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAgcggcg           >  2:2907240/1‑150 (MQ=255)
                                                                                                                                             aCGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGt          <  2:3798314/150‑1 (MQ=255)
                                                                                                                                             aCGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGt          <  1:544233/150‑1 (MQ=255)
                                                                                                                                              cgcgTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTa         >  1:243550/1‑150 (MQ=255)
                                                                                                                                                 gTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCt      <  2:1368419/150‑1 (MQ=255)
                                                                                                                                                   ttcattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTc    <  1:2256247/150‑1 (MQ=255)
                                                                                                                                                    tcattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTcc   <  2:3766717/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: