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% 1264.5 / 42.3 328P272P (CCG→CCTdnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (123/183);  new base A (9/13);  total (132/196)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 333 total aligned reads displayed.

GCCACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGTTT  >  CP009273/3876125‑3876423
                                                                                                                                                     |                                                                                                                                                     
gCCACTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTc                                                                                                                                                       >  1:328255/1‑150 (MQ=255)
    cTTCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGAt                                                                                                                                                   <  1:3172772/150‑1 (MQ=255)
      tCGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAg                                                                                                                                                 >  1:2485175/1‑150 (MQ=255)
       cGCCCGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGc                                                                                                                                                <  1:4209922/150‑1 (MQ=255)
          ccGGCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGAt                                                                                                                                             <  2:431725/150‑1 (MQ=255)
            ggCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATcc                                                                                                                                           >  1:1073022/1‑150 (MQ=255)
            ggCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATcc                                                                                                                                           >  1:1447371/1‑150 (MQ=255)
            ggCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATcc                                                                                                                                           >  1:998089/1‑150 (MQ=255)
            ggCAAACGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATcc                                                                                                                                           <  2:2648385/150‑1 (MQ=255)
               aaaCGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAg                                                                                                                                        >  1:3081267/1‑150 (MQ=255)
                  cGAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTg                                                                                                                                     <  2:1515844/150‑1 (MQ=255)
                   gAATGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTga                                                                                                                                    >  1:4490971/1‑150 (MQ=255)
                    aaTGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTgag                                                                                                                                   >  2:3250741/1‑150 (MQ=255)
                    aaTGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTgag                                                                                                                                   >  2:736350/1‑150 (MQ=255)
                     aTGTCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTgaga                                                                                                                                  >  2:75947/1‑150 (MQ=255)
                        tCGTTTTCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATg                                                                                                                               >  2:1553804/1‑150 (MQ=255)
                           ttttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATc                                                                                                                            >  2:4836156/1‑150 (MQ=255)
                             ttCGTCGGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTg                                                                                                                          <  2:1618752/150‑1 (MQ=255)
                               cgtcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGtt                                                                                                                        <  2:3046340/150‑1 (MQ=255)
                               cgtcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGtt                                                                                                                        >  1:1528349/1‑150 (MQ=255)
                               cgtcgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGtt                                                                                                                        <  1:2993253/150‑1 (MQ=255)
                                  cgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGAt                                                                                                                     >  2:2193418/1‑150 (MQ=255)
                                  cgGCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGAt                                                                                                                     >  1:2205308/1‑150 (MQ=255)
                                    gCCTTTTTCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTa                                                                                                                   <  2:2425253/150‑1 (MQ=255)
                                       tttgtCATCCGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCACCCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCt                                                                                                                <  2:1912851/150‑1 (MQ=255)
                                         tttCATCAGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTc                                                                                                              <  1:2711749/150‑1 (MQ=255)
                                           tcatcaGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCca                                                                                                            <  1:2911711/150‑1 (MQ=255)
                                            catcaGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCcag                                                                                                           >  2:2847329/1‑150 (MQ=255)
                                               caGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCcagcag                                                                                                        <  2:4535741/150‑1 (MQ=255)
                                                aGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAgg                                                                                                       <  1:3529277/150‑1 (MQ=255)
                                                aGGATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAgg                                                                                                       <  1:2286956/150‑1 (MQ=255)
                                                 ggATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAggg                                                                                                      <  1:4575564/150‑1 (MQ=255)
                                                  gATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGc                                                                                                     <  2:1201494/150‑1 (MQ=255)
                                                  gATCGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGc                                                                                                     <  1:1156023/150‑1 (MQ=255)
                                                     cGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGtt                                                                                                  <  1:1770590/150‑1 (MQ=255)
                                                     cGCCACACGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGAGCTGGGGATTACCTTCCAGCAGGGCggt                                                                                                  >  1:4376550/1‑148 (MQ=255)
                                                        cacaCGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGaa                                                                                               >  2:177180/1‑150 (MQ=255)
                                                         acacGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAg                                                                                              >  2:92142/1‑150 (MQ=255)
                                                          cacGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGTGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAgg                                                                                             <  2:3005446/150‑1 (MQ=255)
                                                           acGGGTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgt                                                                                            <  1:1555518/150‑1 (MQ=255)
                                                             gggTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGgtgt                                                                                          <  1:1225963/150‑1 (MQ=255)
                                                              ggTTTCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGtgtg                                                                                         <  1:4076136/150‑1 (MQ=255)
                                                                tttCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGa                                                                                       >  2:1438703/1‑150 (MQ=255)
                                                                  tCCAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaa                                                                                     <  2:3127770/150‑1 (MQ=255)
                                                                    cAGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGaaaaa                                                                                   <  1:3635542/150‑1 (MQ=255)
                                                                     aGCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGTAGCGGGATTTCAATCGATCCTCTACGCAGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGTAGGTGTGGAAAAAc                                                                                  <  1:267172/150‑1 (MQ=255)
                                                                      gCTCTGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAAct                                                                                 <  2:288994/150‑1 (MQ=255)
                                                                        tctGGCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAActct                                                                               <  1:3737628/150‑1 (MQ=255)
                                                                            gCGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCt                                                                           <  2:1348597/150‑1 (MQ=255)
                                                                             cGGTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTAGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTg                                                                          <  1:708140/150‑1 (MQ=255)
                                                                              ggTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTga                                                                         >  1:3815884/1‑150 (MQ=255)
                                                                              ggTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTga                                                                         <  2:2734895/150‑1 (MQ=255)
                                                                               gTTCGATCGCCACAGTCAGTCCCCACCCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTgag                                                                        <  2:4948821/150‑1 (MQ=255)
                                                                               gTTCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTgag                                                                        <  2:3876458/150‑1 (MQ=255)
                                                                                ttCGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTgaga                                                                       <  2:572375/150‑1 (MQ=255)
                                                                                  cGATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATc                                                                     <  1:100888/150‑1 (MQ=255)
                                                                                   gATCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCg                                                                    >  1:659388/1‑150 (MQ=255)
                                                                                    aTCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGt                                                                   >  2:1502264/1‑150 (MQ=255)
                                                                                    aTCGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGt                                                                   >  2:3311873/1‑150 (MQ=255)
                                                                                      cGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTc                                                                 <  1:3315670/150‑1 (MQ=255)
                                                                                      cGCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTc                                                                 >  1:4672053/1‑150 (MQ=255)
                                                                                       gCCACAGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCt                                                                <  1:1829223/150‑1 (MQ=255)
                                                                                         cacaGTCAGTCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCttt                                                              <  2:4795997/150‑1 (MQ=255)
                                                                                           cagtcagtCCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTttat                                                            <  2:226069/150‑1 (MQ=255)
                                                                                            agtcagtcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTttatt                                                           >  1:3804794/1‑150 (MQ=255)
                                                                                            agtcagtcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTttatt                                                           >  2:373807/1‑150 (MQ=255)
                                                                                             gtcagtcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTttatta                                                          <  2:2404244/150‑1 (MQ=255)
                                                                                               cagtcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGc                                                        >  2:4261220/1‑150 (MQ=255)
                                                                                               cagtcCCCAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGc                                                        >  2:2078994/1‑150 (MQ=255)
                                                                                                   ccccAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCaaaa                                                    >  2:710512/1‑150 (MQ=255)
                                                                                                     ccAACCGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCaaaaaa                                                  <  2:4220811/150‑1 (MQ=255)
                                                                                                      cAACCCAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAAc                                                 <  2:1534951/150‑1 (MQ=255)
                                                                                                         ccGAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGa                                              >  1:4926152/1‑150 (MQ=255)
                                                                                                           gAAGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAat                                            <  1:4749086/150‑1 (MQ=255)
                                                                                                            aaGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAata                                           <  1:5072075/150‑1 (MQ=255)
                                                                                                             aGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAatat                                          <  2:464957/150‑1 (MQ=255)
                                                                                                             aGCGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAatat                                          <  1:3043725/150‑1 (MQ=255)
                                                                                                               cGGGATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcg                                        <  2:3531997/150‑1 (MQ=255)
                                                                                                                 ggATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtc                                      >  1:764639/1‑150 (MQ=255)
                                                                                                                 ggATTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATAtcgtc                                      <  2:347974/150‑1 (MQ=255)
                                                                                                                   aTTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGa                                    <  1:2043593/150‑1 (MQ=255)
                                                                                                                   aTTTCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGa                                    <  1:2111799/150‑1 (MQ=255)
                                                                                                                     ttCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATc                                  <  2:3950940/150‑1 (MQ=255)
                                                                                                                     ttCAAACGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATc                                  <  2:4031924/150‑1 (MQ=255)
                                                                                                                        aaaCGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATcagc                               >  1:1705859/1‑150 (MQ=255)
                                                                                                                         aaCGATCCTCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATcagca                              <  2:3643887/150‑1 (MQ=255)
                                                                                                                                 tCAACGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCt                      >  2:1302718/1‑150 (MQ=255)
                                                                                                                                   aaCGCCGTTGATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTAc                    <  2:3927962/150‑1 (MQ=255)
                                                                                                                                   aaCGCCGTTGATCTCTTTCGGATAGCGATCAGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCATGAGTTCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCTTCTAc                    >  2:3371718/1‑150 (MQ=255)
                                                                                                                                            gATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgta           <  1:4734388/150‑1 (MQ=255)
                                                                                                                                            gATCTCTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgta           >  1:2436590/1‑150 (MQ=255)
                                                                                                                                              tctcTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagt         >  1:1865434/1‑150 (MQ=255)
                                                                                                                                              tctcTTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagt         >  2:2619955/1‑150 (MQ=255)
                                                                                                                                               ctctTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGgtagta        <  2:5028549/150‑1 (MQ=255)
                                                                                                                                                tctTTCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGtagtag       <  2:4094457/150‑1 (MQ=255)
                                                                                                                                                    tCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGtt   <  1:2108132/150‑1 (MQ=255)
                                                                                                                                                    tCGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGtt   <  1:4544900/150‑1 (MQ=255)
                                                                                                                                                    tCGGAGAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGtt   <  2:2485175/150‑1 (MQ=255)
                                                                                                                                                     cGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCATTGCATCTACGGAACGGTAGTAGCGttt  <  2:5131844/150‑1 (MQ=255)
                                                                                                                                                     cGGATAGCGATCCGAGGTGAGAATGATCTGTTGATTACCTTCCAGCAGGGCGTTGAAGGTGTGGAAAAACTCTTCCTGAGATCGTTCTTTATTAGCAAAAAACTGAATATCGTCGATCAGCAGTGCATCTACGGAACGGTAGTAGCGttt  <  1:1901593/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: