Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,116 C→A 100% intergenic (‑103/+9) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,1160CA98.1% 906.0 / 6.9 267intergenic (‑103/+9)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (146/116);  total (149/118)
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.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 273 total aligned reads displayed.

GAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGC  >  CP009273/2304969‑2305262
                                                                                                                                                   |                                                                                                                                                  
gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATaaa                                                                                                                                                  >  1:1369041/1‑150 (MQ=255)
 agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAAt                                                                                                                                                 >  1:677057/1‑150 (MQ=255)
   cgcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCg                                                                                                                                               <  1:2035556/150‑1 (MQ=255)
   cgcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCg                                                                                                                                               <  1:1773867/150‑1 (MQ=255)
      ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCgaga                                                                                                                                            >  1:1662206/1‑150 (MQ=255)
         agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGATAAAAGGGCCCGCAGGCGCTTTGTTCGATATAAATCGAGATTa                                                                                                                                         >  2:1019170/1‑150 (MQ=255)
         agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTa                                                                                                                                         >  1:1431922/1‑150 (MQ=255)
         agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTa                                                                                                                                         >  1:1510261/1‑150 (MQ=255)
          gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAg                                                                                                                                        >  1:4157693/1‑150 (MQ=255)
           ccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGa                                                                                                                                       <  1:3396215/150‑1 (MQ=255)
            cagcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTTAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGTTATAAATCGAGATTAGaa                                                                                                                                      <  1:2374773/150‑1 (MQ=255)
            cagcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGaa                                                                                                                                      <  1:756977/150‑1 (MQ=255)
                agTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTgg                                                                                                                                  <  2:2195974/150‑1 (MQ=255)
                agTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTgg                                                                                                                                  >  2:3401138/1‑150 (MQ=255)
                           ggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAcc                                                                                                                       >  1:4369909/1‑150 (MQ=255)
                            gTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGATTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAccc                                                                                                                      <  1:915437/150‑1 (MQ=255)
                            gTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAccc                                                                                                                      >  2:3576805/1‑150 (MQ=255)
                              aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCag                                                                                                                    >  2:810256/1‑150 (MQ=255)
                              aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCGTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCag                                                                                                                    >  2:4045739/1‑150 (MQ=255)
                                aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  <  2:565784/150‑1 (MQ=255)
                                 gCTTATTTCCATACTAGAAAATCCTCATGCGAACGGTAGCAAGAGTACACCAAAAAACTCAGTTTCTACGCTGAAAACAATGAAAAAAGGGACCGCAGTCCCATAGTTCGATATACATCGAGATTAGAACTGGAAAACCAGACCCAGAGc                                                                                                                 >  2:3487258/1‑150 (MQ=255)
                                 gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGc                                                                                                                 >  1:1976147/1‑150 (MQ=255)
                                     aTTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACg                                                                                                             >  2:3450446/1‑150 (MQ=255)
                                            tACGTGATGATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTAt                                                                                                      <  2:3256002/150‑1 (MQ=255)
                                              cGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCa                                                                                                    <  1:2897134/150‑1 (MQ=255)
                                               gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                   >  2:2679015/1‑150 (MQ=255)
                                               gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                   >  1:3499810/1‑150 (MQ=255)
                                                tGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                                                  <  1:1129097/150‑1 (MQ=255)
                                                 gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                                                 >  2:4317657/1‑150 (MQ=255)
                                                 gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                                                 >  2:4317673/1‑150 (MQ=255)
                                                 gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                                                 >  1:1520144/1‑150 (MQ=255)
                                                  aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCCGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                                                >  2:453450/1‑150 (MQ=255)
                                                  aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                                                >  2:1932755/1‑150 (MQ=255)
                                                  aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATGAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                                                <  1:1264294/150‑1 (MQ=255)
                                                     aTCCTCATGCGAACGGTCGCAAGAGTACACCAAAAATCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTGCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGa                                                                                             >  2:1592819/1‑150 (MQ=255)
                                                      tCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                                            >  1:2057261/1‑150 (MQ=255)
                                                        cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                          >  1:4147676/1‑150 (MQ=255)
                                                        cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                          <  1:2224285/150‑1 (MQ=255)
                                                            tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCCGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                                                      <  1:972104/150‑1 (MQ=255)
                                                              cGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAATAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                                    >  2:1888384/1‑150 (MQ=255)
                                                              cGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                                    >  2:2265699/1‑150 (MQ=255)
                                                               gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTc                                                                                   >  1:263962/1‑150 (MQ=255)
                                                                aaCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCa                                                                                  <  2:2965109/150‑1 (MQ=255)
                                                                 aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCAc                                                                                 >  2:3086412/1‑150 (MQ=255)
                                                                 aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCAc                                                                                 >  2:841524/1‑150 (MQ=255)
                                                                 aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCAc                                                                                 <  2:2542611/150‑1 (MQ=255)
                                                                 aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCAc                                                                                 <  1:1299573/150‑1 (MQ=255)
                                                                  cGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACg                                                                                >  1:4475308/1‑150 (MQ=255)
                                                                   ggTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTAGCAGTGTTGATGCCAGCGTCACGa                                                                               >  1:3086003/1‑150 (MQ=255)
                                                                     tCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGt                                                                             <  2:454809/150‑1 (MQ=255)
                                                                     tCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGt                                                                             <  2:4162315/150‑1 (MQ=255)
                                                                        cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                                          >  1:2920838/1‑150 (MQ=255)
                                                                        cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                                          >  2:3962319/1‑150 (MQ=255)
                                                                          agagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACt                                                                        >  1:686642/1‑150 (MQ=255)
                                                                          agagAACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACt                                                                        <  1:1626709/150‑1 (MQ=255)
                                                                            agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                      <  2:2365299/150‑1 (MQ=255)
                                                                              tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                    >  1:1102631/1‑150 (MQ=255)
                                                                              tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                    >  1:2162316/1‑150 (MQ=255)
                                                                                  ccAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                                                <  1:995011/150‑1 (MQ=255)
                                                                                  ccAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                                                >  1:4025406/1‑150 (MQ=255)
                                                                                    aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGGTGATGCCAGCGTCACGAGTGAACTGGTTgtcgtc                                                              >  2:104311/1‑150 (MQ=255)
                                                                                     aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTcc                                                             >  1:1897455/1‑150 (MQ=255)
                                                                                      aaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCca                                                            <  2:1152252/150‑1 (MQ=255)
                                                                                      aaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCca                                                            <  2:601908/150‑1 (MQ=255)
                                                                                         aCTGAGTTGGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                         <  1:176981/150‑1 (MQ=255)
                                                                                          cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagcag                                                        >  2:1438913/1‑150 (MQ=255)
                                                                                           tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAgg                                                       >  2:2764835/1‑150 (MQ=255)
                                                                                             aGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                     <  1:1591135/150‑1 (MQ=255)
                                                                                             aGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                     <  2:1697707/150‑1 (MQ=255)
                                                                                               tttGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGa                                                   <  1:962472/150‑1 (MQ=255)
                                                                                                 tGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtt                                                 <  1:2391721/150‑1 (MQ=255)
                                                                                                     cGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGt                                             >  2:630952/1‑150 (MQ=255)
                                                                                                          aaaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCa                                        <  2:2633162/150‑1 (MQ=255)
                                                                                                           aaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCaa                                       >  1:4549888/1‑150 (MQ=255)
                                                                                                            aaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAAc                                      <  2:1047696/150‑1 (MQ=255)
                                                                                                             aCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg                                     >  2:125901/1‑150 (MQ=255)
                                                                                                             aCAATGAAAAAAGGGCCCGCAGGCCATTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg                                     >  1:310665/1‑150 (MQ=255)
                                                                                                               aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTa                                   <  1:874560/150‑1 (MQ=255)
                                                                                                                  gAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt                                <  2:585011/150‑1 (MQ=255)
                                                                                                                    aaaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTgg                              >  1:637924/1‑150 (MQ=255)
                                                                                                                    aaaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTgg                              >  1:2441703/1‑150 (MQ=255)
                                                                                                                     aaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                             >  1:3656873/1‑150 (MQ=255)
                                                                                                                     aaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                             >  1:2846233/1‑150 (MQ=255)
                                                                                                                     aaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                             >  1:3173182/1‑150 (MQ=255)
                                                                                                                      aaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAc                            <  1:1829104/150‑1 (MQ=255)
                                                                                                                       aaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACa                           <  1:1010517/150‑1 (MQ=255)
                                                                                                                        aGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACAt                          >  1:2730851/1‑150 (MQ=255)
                                                                                                                            cccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttt                      >  2:4433918/1‑150 (MQ=255)
                                                                                                                               gCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg                   >  2:3516290/1‑150 (MQ=255)
                                                                                                                                cAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgt                  <  2:2993654/150‑1 (MQ=255)
                                                                                                                                cAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgt                  <  1:887586/150‑1 (MQ=255)
                                                                                                                                  ggCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                <  2:174009/150‑1 (MQ=255)
                                                                                                                                  ggCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                <  2:173980/150‑1 (MQ=255)
                                                                                                                                    cccTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa              >  2:4257946/1‑150 (MQ=255)
                                                                                                                                      cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt            >  2:352100/1‑150 (MQ=255)
                                                                                                                                      cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt            <  2:2710154/150‑1 (MQ=255)
                                                                                                                                      cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt            <  2:1538855/150‑1 (MQ=255)
                                                                                                                                         tGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt         >  2:2885189/1‑150 (MQ=255)
                                                                                                                                             cGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggt     >  1:1814417/1‑150 (MQ=255)
                                                                                                                                              gATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta    >  1:3312150/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGC  >  CP009273/2304969‑2305262

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: