Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,876,274 C→A 6.7% P272P (CCG→CCT dnaA ← chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,876,2740CA 6.7% 970.0 / 29.4 261P272P (CCG→CCTdnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (116/127);  new base A (10/8);  total (126/135)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

GCCACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGTTT  >  CP009273/3876125‑3876423
                                                                                                                                                     |                                                                                                                                                     
gCCACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTc                                                                                                                                                       <  2:3916891/150‑1 (MQ=255)
 ccACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCg                                                                                                                                                      >  1:3934533/1‑150 (MQ=255)
 ccACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCg                                                                                                                                                      >  1:4135879/1‑150 (MQ=255)
   aCTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGa                                                                                                                                                    >  1:3635529/1‑150 (MQ=255)
     ttCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATa                                                                                                                                                  <  1:3817488/150‑1 (MQ=255)
     ttCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATa                                                                                                                                                  <  1:3817406/150‑1 (MQ=255)
     ttCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATa                                                                                                                                                  <  1:909568/150‑1 (MQ=255)
     ttCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATa                                                                                                                                                  <  1:3817277/150‑1 (MQ=255)
      tCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAg                                                                                                                                                 >  2:2416642/1‑150 (MQ=255)
        gCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCg                                                                                                                                               >  1:2455281/1‑150 (MQ=255)
              cAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGa                                                                                                                                         <  1:1738554/150‑1 (MQ=255)
                  cGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTg                                                                                                                                     <  1:962283/150‑1 (MQ=255)
                  cGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTg                                                                                                                                     >  1:3086702/1‑150 (MQ=255)
                   gAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTga                                                                                                                                    >  1:2775325/1‑150 (MQ=255)
                     aTGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTgaga                                                                                                                                  >  2:3295388/1‑150 (MQ=255)
                     aTGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTgaga                                                                                                                                  >  2:1016421/1‑150 (MQ=255)
                       gTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAAt                                                                                                                                >  1:1835791/1‑150 (MQ=255)
                         cGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACTGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTAAGGATAGCGATCCGAGGTGAGAATGa                                                                                                                              >  2:485133/1‑150 (MQ=255)
                          gTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGAt                                                                                                                             >  1:3503520/1‑150 (MQ=255)
                          gTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGAt                                                                                                                             <  1:450048/150‑1 (MQ=255)
                           ttttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAATCGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATc                                                                                                                            <  2:1549025/150‑1 (MQ=255)
                            tttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCt                                                                                                                           <  1:3928599/150‑1 (MQ=255)
                             ttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTg                                                                                                                          <  1:2750795/150‑1 (MQ=255)
                             ttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTg                                                                                                                          >  2:3489926/1‑150 (MQ=255)
                               cgtcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGtt                                                                                                                        <  2:2477995/150‑1 (MQ=255)
                               cgtcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGtt                                                                                                                        <  2:103756/150‑1 (MQ=255)
                                 tcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGa                                                                                                                      >  2:4538240/1‑150 (MQ=255)
                                    gCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCTATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTa                                                                                                                   <  2:4037867/150‑1 (MQ=255)
                                    gCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTa                                                                                                                   <  2:3060509/150‑1 (MQ=255)
                                       tttttCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCt                                                                                                                >  2:1469057/1‑150 (MQ=255)
                                       tttttCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCt                                                                                                                <  2:1847067/150‑1 (MQ=255)
                                        ttttCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCtt                                                                                                               <  1:2574931/150‑1 (MQ=255)
                                        ttttCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCtt                                                                                                               <  2:355648/150‑1 (MQ=255)
                                           tcatcaGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCca                                                                                                            >  1:673267/1‑150 (MQ=255)
                                               caGGATCGCCACACGGGTTACCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGAACCTCAACGCCGATGATCTCATTCGGATAGCGATCCTAGGTGAGACTGATCTGTTGATTACCTTCcagcag                                                                                                        >  2:311545/1‑150 (MQ=255)
                                                aGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAgg                                                                                                       >  2:2908475/1‑150 (MQ=255)
                                                   aTCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCg                                                                                                    <  1:2541707/150‑1 (MQ=255)
                                                     cGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGtt                                                                                                  <  2:3712651/150‑1 (MQ=255)
                                                     cGCCACACGGGTGTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGtt                                                                                                  <  2:3712702/150‑1 (MQ=255)
                                                      gCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTg                                                                                                 <  1:1236743/150‑1 (MQ=255)
                                                       ccACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGa                                                                                                >  1:4434434/1‑150 (MQ=255)
                                                         acacGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAg                                                                                              <  1:1193674/150‑1 (MQ=255)
                                                          cacGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAgg                                                                                             >  2:3998509/1‑150 (MQ=255)
                                                          cacGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAgg                                                                                             >  2:4320765/1‑150 (MQ=255)
                                                           acGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgt                                                                                            <  2:358401/150‑1 (MQ=255)
                                                           acGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgt                                                                                            >  1:355375/1‑150 (MQ=255)
                                                           acGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgt                                                                                            <  2:1667889/150‑1 (MQ=255)
                                                            cGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgtg                                                                                           <  2:2546570/150‑1 (MQ=255)
                                                             gggTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgtgt                                                                                          >  1:2689371/1‑150 (MQ=255)
                                                                 ttCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaa                                                                                      <  1:3567713/150‑1 (MQ=255)
                                                                 ttCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaa                                                                                      <  1:3567184/150‑1 (MQ=255)
                                                                  tCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaa                                                                                     <  1:323157/150‑1 (MQ=255)
                                                                   ccAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaaa                                                                                    >  2:2820394/1‑150 (MQ=255)
                                                                   ccAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaaa                                                                                    <  1:2407461/150‑1 (MQ=255)
                                                                   ccAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaaa                                                                                    <  1:1807224/150‑1 (MQ=255)
                                                                     aGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGACTGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAAc                                                                                  >  2:197505/1‑150 (MQ=255)
                                                                         ctGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCtt                                                                              >  2:2671267/1‑150 (MQ=255)
                                                                          tGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTc                                                                             >  1:1951553/1‑150 (MQ=255)
                                                                            gCGGTTCGCTCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGATGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGATAAACTCTTCCt                                                                           <  2:707390/150‑1 (MQ=255)
                                                                             cGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTg                                                                          <  1:2619900/150‑1 (MQ=255)
                                                                             cGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTg                                                                          >  1:129786/1‑150 (MQ=255)
                                                                                 tCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGAt                                                                      >  2:4511197/1‑150 (MQ=255)
                                                                                  cGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATc                                                                     <  1:2099762/150‑1 (MQ=255)
                                                                                  cGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATc                                                                     <  2:413818/150‑1 (MQ=255)
                                                                                  cGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATc                                                                     >  2:3691391/1‑150 (MQ=255)
                                                                                   gATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCg                                                                    <  1:2399294/150‑1 (MQ=255)
                                                                                         cacaGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCttt                                                              <  2:2713539/150‑1 (MQ=255)
                                                                                         cacaGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCttt                                                              <  2:2576196/150‑1 (MQ=255)
                                                                                           cagtcagtCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTttat                                                            >  2:4251722/1‑150 (MQ=255)
                                                                                                  tcCCCTACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGTTCCGAGGTGAGAATGATCTGTTGATTACCTTCCATCTTGGCGTTGATGGTGTGGAATAACTCTTCCTGAGATCGTTCTTTATTAGCaaa                                                     <  1:1206200/150‑1 (MQ=255)
                                                                                                  tcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCaaa                                                     >  2:1065030/1‑150 (MQ=255)
                                                                                                   ccccAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCaaaa                                                    >  1:3997953/1‑150 (MQ=255)
                                                                                                      cAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAAc                                                 <  1:2978966/150‑1 (MQ=255)
                                                                                                      cAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAAc                                                 >  2:1487471/1‑150 (MQ=255)
                                                                                                      cAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAAc                                                 <  1:1565847/150‑1 (MQ=255)
                                                                                                       aaCCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACt                                                <  2:4499338/150‑1 (MQ=255)
                                                                                                          cGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGaa                                             >  1:3392980/1‑150 (MQ=255)
                                                                                                           gAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAat                                            >  2:1808630/1‑150 (MQ=255)
                                                                                                             aGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAatat                                          <  1:3564393/150‑1 (MQ=255)
                                                                                                              gCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGTTAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtc                                         <  2:4038886/150‑1 (MQ=255)
                                                                                                              gCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtc                                         <  1:2839236/150‑1 (MQ=255)
                                                                                                                 ggATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtc                                      >  2:314278/1‑150 (MQ=255)
                                                                                                                 ggATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtc                                      >  2:1199002/1‑150 (MQ=255)
                                                                                                                 ggATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGGTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtc                                      >  2:1198899/1‑150 (MQ=255)
                                                                                                                  gATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtcg                                     >  1:2310912/1‑150 (MQ=255)
                                                                                                                  gATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtcg                                     >  1:220857/1‑150 (MQ=255)
                                                                                                                   aTTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGa                                    >  2:3415710/1‑150 (MQ=255)
                                                                                                                         aaCGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATcagca                              <  1:4251722/150‑1 (MQ=255)
                                                                                                                          aCGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATcagcag                             <  1:754787/150‑1 (MQ=255)
                                                                                                                            gATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTg                           <  1:516792/150‑1 (MQ=255)
                                                                                                                             aTCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGc                          <  1:1374314/150‑1 (MQ=255)
                                                                                                                                  cAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTa                     <  1:2725580/150‑1 (MQ=255)
                                                                                                                                   aaCGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTAc                    >  1:3328174/1‑150 (MQ=255)
                                                                                                                                         gTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACg              <  2:893735/150‑1 (MQ=255)
                                                                                                                                              tctcTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagt         >  2:3885269/1‑150 (MQ=255)
                                                                                                                                              tctcTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagt         <  1:1613356/150‑1 (MQ=255)
                                                                                                                                               ctctTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagta        <  2:673267/150‑1 (MQ=255)
                                                                                                                                                tctTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGtagtag       >  2:3758610/1‑150 (MQ=255)
                                                                                                                                                tctTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGtagtag       <  1:4077196/150‑1 (MQ=255)
                                                                                                                                                 ctTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGc      <  2:2215942/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGTTT  >  CP009273/3876125‑3876423

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: