Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,276,667 C→A 7.8% E254* (GAA→TAA)  yejK ← spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,276,6670CA 7.8% 1048.8 / 40.1 282E254* (GAA→TAA) yejKspermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (111/149);  new base A (8/14);  total (119/163)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 285 total aligned reads displayed.

TTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCG  >  CP009273/2276518‑2276816
                                                                                                                                                     |                                                                                                                                                     
     cTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTa                                                                                                                                                  >  1:1434684/1‑150 (MQ=255)
     cTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTa                                                                                                                                                  <  1:75683/150‑1 (MQ=255)
     cTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTa                                                                                                                                                  <  1:2673539/150‑1 (MQ=255)
       aaCAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACg                                                                                                                                                >  2:4149465/1‑150 (MQ=255)
        aCAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGc                                                                                                                                               >  1:4061530/1‑150 (MQ=255)
            catcgcatcAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGc                                                                                                                                           >  1:2091869/1‑150 (MQ=255)
             atcgcatcAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGcc                                                                                                                                          >  1:2138886/1‑150 (MQ=255)
                gcatcAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGc                                                                                                                                       >  2:3921140/1‑150 (MQ=255)
                    cAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCtt                                                                                                                                   >  2:2497720/1‑150 (MQ=255)
                      aaGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                                                                                                                                 >  1:185427/1‑150 (MQ=255)
                        gTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGa                                                                                                                               >  2:891853/1‑150 (MQ=255)
                         ttAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGAc                                                                                                                              <  2:3865753/150‑1 (MQ=255)
                         ttAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGAc                                                                                                                              >  1:776378/1‑150 (MQ=255)
                           aaTAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAg                                                                                                                            >  1:241833/1‑150 (MQ=255)
                             tAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCg                                                                                                                          >  2:2154150/1‑150 (MQ=255)
                              aGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGa                                                                                                                         >  2:3811776/1‑150 (MQ=255)
                               gTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAt                                                                                                                        >  1:3865032/1‑150 (MQ=255)
                                   aGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCa                                                                                                                    <  1:2059482/150‑1 (MQ=255)
                                   aGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCa                                                                                                                    <  1:3380273/150‑1 (MQ=255)
                                    gCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAg                                                                                                                   >  2:1814796/1‑150 (MQ=255)
                                    gCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAg                                                                                                                   <  2:1215755/150‑1 (MQ=255)
                                      cACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGACGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGtt                                                                                                                 >  2:3840639/1‑150 (MQ=255)
                                          gccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGAt                                                                                                             >  1:4120171/1‑150 (MQ=255)
                                          gccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGAt                                                                                                             >  1:1869930/1‑150 (MQ=255)
                                                aCCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTc                                                                                                       >  1:462264/1‑150 (MQ=255)
                                                aCCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTc                                                                                                       >  2:553087/1‑150 (MQ=255)
                                                 ccGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCg                                                                                                      <  1:3776424/150‑1 (MQ=255)
                                                    gCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCgcct                                                                                                   <  1:1175089/150‑1 (MQ=255)
                                                    gCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCgcct                                                                                                   <  1:2444042/150‑1 (MQ=255)
                                                    gCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCgcct                                                                                                   <  2:1566783/150‑1 (MQ=255)
                                                     cAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCgcctg                                                                                                  <  1:1622390/150‑1 (MQ=255)
                                                         cTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGcctgcctg                                                                                              >  2:2356269/1‑150 (MQ=255)
                                                           tCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCa                                                                                            <  1:753859/150‑1 (MQ=255)
                                                           tCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCa                                                                                            <  1:753957/150‑1 (MQ=255)
                                                              tCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgct                                                                                         <  2:1165775/150‑1 (MQ=255)
                                                               cAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgctg                                                                                        >  1:1666970/1‑150 (MQ=255)
                                                                 aCTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgctgct                                                                                      >  1:2745473/1‑150 (MQ=255)
                                                                 aCTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgctgct                                                                                      <  2:4477632/150‑1 (MQ=255)
                                                                  cTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTc                                                                                     >  2:1990548/1‑150 (MQ=255)
                                                                   tGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCa                                                                                    <  1:3999139/150‑1 (MQ=255)
                                                                     aCGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCAtt                                                                                  <  1:603772/150‑1 (MQ=255)
                                                                      cGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                 <  2:344468/150‑1 (MQ=255)
                                                                      cGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                 <  2:2641282/150‑1 (MQ=255)
                                                                      cGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                 <  1:1000012/150‑1 (MQ=255)
                                                                      cGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                 <  2:4314886/150‑1 (MQ=255)
                                                                      cGCAGTGTGCTACGGTCTGCCGGGAAGCTCACTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                 <  2:4343653/150‑1 (MQ=255)
                                                                       gCAGTGTGCTACGGTCTGCCGGTAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTAc                                                                                <  2:4164973/150‑1 (MQ=255)
                                                                        cAGTGTGCTACTGTCTGCCGGGAAGCTCTCGTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACa                                                                               <  1:4343286/150‑1 (MQ=255)
                                                                            gtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAg                                                                           <  2:1669676/150‑1 (MQ=255)
                                                                            gtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGGTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAg                                                                           >  2:1906120/1‑150 (MQ=255)
                                                                             tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGc                                                                          >  2:2941872/1‑150 (MQ=255)
                                                                             tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGc                                                                          <  1:4374055/150‑1 (MQ=255)
                                                                               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCta                                                                        <  2:3163091/150‑1 (MQ=255)
                                                                               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCta                                                                        <  1:3921140/150‑1 (MQ=255)
                                                                                tACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACa                                                                               <  2:3066750/142‑1 (MQ=255)
                                                                                tACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACa                                                                               >  1:3066750/1‑142 (MQ=255)
                                                                                tACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCtat                                                                       >  1:2010190/1‑150 (MQ=255)
                                                                                    gTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCTCTTCGCTTTCGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGGTGCTCATTACAGTAGCTATAAAc                                                                   <  1:3699432/150‑1 (MQ=255)
                                                                                      cTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACct                                                                 <  1:3550161/150‑1 (MQ=255)
                                                                                         ccGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACctgct                                                              <  1:1814796/150‑1 (MQ=255)
                                                                                         ccGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACctgct                                                              <  2:233622/150‑1 (MQ=255)
                                                                                          cGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACcagcta                                                             >  2:298240/1‑149 (MQ=255)
                                                                                            ggAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGgc                                                           <  1:2538275/150‑1 (MQ=255)
                                                                                              aaGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGgcgc                                                         >  2:4218989/1‑150 (MQ=255)
                                                                                              aaGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGgcgc                                                         >  2:4218939/1‑150 (MQ=255)
                                                                                               aGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCa                                                        <  2:1142677/150‑1 (MQ=255)
                                                                                                 ctctctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACg                                                      <  2:4508351/150‑1 (MQ=255)
                                                                                                  tctctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGt                                                     >  2:1957984/1‑150 (MQ=255)
                                                                                                     ctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCt                                                  <  1:1336389/150‑1 (MQ=255)
                                                                                                     ctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCt                                                  <  2:1759408/150‑1 (MQ=255)
                                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTg                                                 >  2:627720/1‑150 (MQ=255)
                                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTg                                                 <  1:769316/150‑1 (MQ=255)
                                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTg                                                 <  1:2286400/150‑1 (MQ=255)
                                                                                                       tCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGt                                                <  2:1843742/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTc                                               <  2:2745473/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTc                                               <  1:3599616/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTc                                               <  1:314819/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTc                                               <  2:490168/150‑1 (MQ=255)
                                                                                                         cAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCg                                              <  2:4106046/150‑1 (MQ=255)
                                                                                                         cAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCg                                              <  2:925095/150‑1 (MQ=255)
                                                                                                               gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGc                                        >  1:1204066/1‑150 (MQ=255)
                                                                                                               gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGc                                        >  1:3476028/1‑150 (MQ=255)
                                                                                                                   ccTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTa                                    >  2:3536547/1‑150 (MQ=255)
                                                                                                                   ccTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTa                                    >  2:1451103/1‑150 (MQ=255)
                                                                                                                      ttttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATcc                                 <  2:372871/150‑1 (MQ=255)
                                                                                                                      ttttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATcc                                 <  1:2410112/150‑1 (MQ=255)
                                                                                                                      ttttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATcc                                 <  1:2410400/150‑1 (MQ=255)
                                                                                                                          cAGCCGCAAACTCCGTGAAGCTCACTTAGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGtt                             >  1:2871609/1‑150 (MQ=255)
                                                                                                                              cGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGc                         <  1:2962214/150‑1 (MQ=255)
                                                                                                                                 aaaCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTc                      >  2:743362/1‑150 (MQ=255)
                                                                                                                                  aaCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCg                     <  1:105874/150‑1 (MQ=255)
                                                                                                                                  aaCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCg                     >  1:1554102/1‑150 (MQ=255)
                                                                                                                                   aCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCgg                    >  1:3011749/1‑150 (MQ=255)
                                                                                                                                        gTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTg               <  1:741571/150‑1 (MQ=255)
                                                                                                                                               tCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAatcat        >  2:28256/1‑150 (MQ=255)
                                                                                                                                               tCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAatcat        >  2:3275338/1‑150 (MQ=255)
                                                                                                                                                 aCTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAAtcatca      <  1:3357650/150‑1 (MQ=255)
                                                                                                                                                    tCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAAcc   <  2:1015120/150‑1 (MQ=255)
                                                                                                                                                     cGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCg  >  1:1764782/1‑150 (MQ=255)
                                                                                                                                                     cGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCg  <  2:2047620/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCG  >  CP009273/2276518‑2276816

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: