Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,514,074 C→A 21.4% E252* (GAA→TAA)  yjhH ← putative lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,514,0740CA21.4% 513.5 / 90.9 195E252* (GAA→TAA) yjhHputative lyase/synthase
Reads supporting (aligned to +/- strand):  ref base C (93/60);  new base A (23/19);  total (116/79)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

GTAAAATGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAG  >  CP009273/4513930‑4514221
                                                                                                                                                |                                                                                                                                                   
gTAAAATGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAg                                                                                                                                                  <  1:885319/148‑1 (MQ=255)
 tAAAATGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAgg                                                                                                                                                 <  2:2343825/148‑1 (MQ=255)
  aaaaTGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAggg                                                                                                                                                <  1:3055569/148‑1 (MQ=255)
     aTGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGt                                                                                                                                             >  1:2545819/1‑148 (MQ=255)
      tGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTa                                                                                                                                            >  1:2046669/1‑148 (MQ=255)
            gcgcGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGc                                                                                                                                      >  1:1054322/1‑148 (MQ=255)
                   aaGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGt                                                                                                                               >  2:520150/1‑148 (MQ=255)
                   aaGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGt                                                                                                                               <  2:299222/148‑1 (MQ=255)
                    agcagcACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGtt                                                                                                                              >  1:19520/1‑148 (MQ=255)
                     gcagcaCGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTg                                                                                                                             >  2:374228/1‑148 (MQ=255)
                       agcaCGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTgta                                                                                                                           <  1:3018064/148‑1 (MQ=255)
                         caCGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTgtagt                                                                                                                         >  2:1139242/1‑148 (MQ=255)
                            gTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTagtagt                                                                                                                      <  2:1116878/148‑1 (MQ=255)
                             tGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtt                                                                                                                     <  2:1089400/148‑1 (MQ=255)
                               gACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                                                                                   <  2:1311788/148‑1 (MQ=255)
                               gACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                                                                                   <  2:2255486/148‑1 (MQ=255)
                                aCTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGttttt                                                                                                                  >  1:491384/1‑148 (MQ=255)
                                aCTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGttttt                                                                                                                  >  1:2154622/1‑148 (MQ=255)
                                 cTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttttt                                                                                                                 <  2:178397/148‑1 (MQ=255)
                                 cTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttttt                                                                                                                 <  2:2083560/148‑1 (MQ=255)
                                   ttATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAt                                                                                                               >  1:566882/1‑148 (MQ=255)
                                   ttATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAt                                                                                                               >  1:414701/1‑148 (MQ=255)
                                   ttATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAt                                                                                                               <  1:1736662/148‑1 (MQ=255)
                                      tCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCa                                                                                                            >  1:1921270/1‑148 (MQ=255)
                                       cTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGTCAGCCCGACACACTGCATGTTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGAGGGCAGTTGTAGTAGTTTTTTATTCAg                                                                                                           >  1:2582283/1‑148 (MQ=255)
                                       cTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGTTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAg                                                                                                           >  1:925950/1‑148 (MQ=255)
                                         tttGCTTCTTCAGATTCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCg                                                                                                         >  2:2079956/1‑148 (MQ=255)
                                          ttGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGt                                                                                                        >  2:3160147/1‑148 (MQ=255)
                                               tcttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAg                                                                                                   >  1:3300437/1‑148 (MQ=255)
                                               tcttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAg                                                                                                   >  1:3300644/1‑148 (MQ=255)
                                               tcttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAg                                                                                                   >  2:3146423/1‑148 (MQ=255)
                                                cttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGc                                                                                                  >  2:1767774/1‑148 (MQ=255)
                                                cttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGGTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGc                                                                                                  >  1:1857509/1‑148 (MQ=255)
                                                 ttcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACTCTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCg                                                                                                 <  2:665788/148‑1 (MQ=255)
                                                  tcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCgg                                                                                                >  2:762543/1‑148 (MQ=255)
                                                   cAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGt                                                                                               >  2:1997073/1‑148 (MQ=255)
                                                      aTGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGc                                                                                            >  1:1906532/1‑148 (MQ=255)
                                                      aTGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGc                                                                                            >  1:1395754/1‑148 (MQ=255)
                                                       tGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGcc                                                                                           >  2:2562697/1‑148 (MQ=255)
                                                       tGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGcc                                                                                           >  2:2967808/1‑148 (MQ=255)
                                                          ttCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGa                                                                                        <  1:789063/148‑1 (MQ=255)
                                                           tCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt                                                                                       >  2:2187336/1‑148 (MQ=255)
                                                                aaTCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCt                                                                                  >  1:2870685/1‑148 (MQ=255)
                                                                aaTCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCt                                                                                  >  2:2511474/1‑148 (MQ=255)
                                                                 aTCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCtt                                                                                 >  2:2242810/1‑148 (MQ=255)
                                                                 aTCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCtt                                                                                 >  1:768852/1‑148 (MQ=255)
                                                                  tCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTc                                                                                <  2:422317/148‑1 (MQ=255)
                                                                  tCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTc                                                                                >  1:1025932/1‑148 (MQ=255)
                                                                  tCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTc                                                                                >  1:2678748/1‑148 (MQ=255)
                                                                  tCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTc                                                                                >  1:1025929/1‑148 (MQ=255)
                                                                   cGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCa                                                                               >  2:2586611/1‑148 (MQ=255)
                                                                   cGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCa                                                                               >  2:2459385/1‑148 (MQ=255)
                                                                     gtggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACg                                                                             >  2:1039121/1‑148 (MQ=255)
                                                                      tggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGc                                                                            <  1:2742527/148‑1 (MQ=255)
                                                                         tAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAg                                                                         >  2:818394/1‑148 (MQ=255)
                                                                          aaGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAgg                                                                        >  1:2165004/1‑148 (MQ=255)
                                                                               atatGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg                                                                   >  1:1026268/1‑148 (MQ=255)
                                                                                tatGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGt                                                                  <  1:2002736/148‑1 (MQ=255)
                                                                                     ctctACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATg                                                             <  1:2192192/148‑1 (MQ=255)
                                                                                     ctctACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATg                                                             <  1:1133882/148‑1 (MQ=255)
                                                                                     ctctACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGGTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATg                                                             >  2:1095862/1‑148 (MQ=255)
                                                                                         aCAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGa                                                         <  1:1119743/148‑1 (MQ=255)
                                                                                          cAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGAc                                                        >  1:3114567/1‑148 (MQ=255)
                                                                                           aGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg                                                       <  1:34606/148‑1 (MQ=255)
                                                                                             gCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATTCCGACGga                                                     >  1:7325/1‑148 (MQ=255)
                                                                                             gCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGga                                                     <  2:2006188/148‑1 (MQ=255)
                                                                                              cAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGgag                                                    <  1:374228/148‑1 (MQ=255)
                                                                                                gCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGagag                                                  >  1:2604893/1‑148 (MQ=255)
                                                                                                gCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGagag                                                  >  1:2604078/1‑148 (MQ=255)
                                                                                                 cccGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGt                                                 <  1:3146423/148‑1 (MQ=255)
                                                                                                  ccGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGtt                                                >  2:1078744/1‑148 (MQ=255)
                                                                                                   cGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTc                                               >  1:1750713/1‑148 (MQ=255)
                                                                                                      cacacTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCgg                                            >  1:1644702/1‑148 (MQ=255)
                                                                                                      cacacTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCgg                                            >  2:2387242/1‑148 (MQ=255)
                                                                                                           tGCATGCTGTATTTGCTCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGGAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCaa                                       <  2:2700253/148‑1 (MQ=255)
                                                                                                           ggCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCaa                                       <  2:2271651/147‑1 (MQ=255)
                                                                                                            gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCaaa                                      <  2:3300479/148‑1 (MQ=255)
                                                                                                              aTGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGt                                    <  2:872568/148‑1 (MQ=255)
                                                                                                                gCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTa                                  >  1:2468938/1‑148 (MQ=255)
                                                                                                                gCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTa                                  >  2:1239087/1‑148 (MQ=255)
                                                                                                                 cTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAg                                 >  1:3111880/1‑148 (MQ=255)
                                                                                                                 cTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAg                                 <  1:430760/148‑1 (MQ=255)
                                                                                                                 cTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAg                                 <  1:77559/148‑1 (MQ=255)
                                                                                                                      tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTg                            <  1:762543/148‑1 (MQ=255)
                                                                                                                      tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTg                            <  2:161645/148‑1 (MQ=255)
                                                                                                                      tttGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTg                            <  1:1474231/148‑1 (MQ=255)
                                                                                                                       ttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTgg                           >  2:1908190/1‑148 (MQ=255)
                                                                                                                        tGATCAGTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGc                          <  2:611283/148‑1 (MQ=255)
                                                                                                                           tCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGt                       >  1:2311566/1‑148 (MQ=255)
                                                                                                                           tCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGt                       >  1:3090576/1‑148 (MQ=255)
                                                                                                                            cAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGtt                      <  2:709912/148‑1 (MQ=255)
                                                                                                                            cAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGtt                      <  1:401612/148‑1 (MQ=255)
                                                                                                                              gTGAGACAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTAt                    >  1:1324222/1‑148 (MQ=255)
                                                                                                                                 agaCAAACGGTGTTTAGAGGGCGTAAATAGCGGGCAGTAGTAGTAGTTTTTTATTCAGCGTCGCACCGGTCGCCAGAACCCCTTCACTCCAGGCGCGGTAGAAGCCGACGGAGAGTTCCGGAGCAAAGTTAGCCCTTGCGGTTATCGc                 >  1:2355210/1‑148 (MQ=255)
                                                                                                                                   aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCAc               >  1:1811119/1‑148 (MQ=255)
                                                                                                                                   aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCAc               >  2:1806525/1‑148 (MQ=255)
                                                                                                                                   aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCAc               >  2:1341548/1‑148 (MQ=255)
                                                                                                                                           tgtTTAGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTcgccg       <  1:855755/148‑1 (MQ=255)
                                                                                                                                             tttCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCgccgcc     >  1:169931/1‑148 (MQ=255)
                                                                                                                                             tttCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCgccgcc     <  2:491384/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTAAAATGCCCTGCGCGGTAAGCAGCACGTGGACTTTATCTTTTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAG  >  CP009273/4513930‑4514221

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: