Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,358,920 G→T 100% R147L (CGC→CTC)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,358,9200GT97.6% 527.1 / 3.4 170R147L (CGC→CTC) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base T (73/92);  total (78/92)
Fisher's exact test for biased strand distribution p-value = 4.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 171 total aligned reads displayed.

CTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCC  >  CP009273/4358775‑4359061
                                                                                                                                                 |                                                                                                                                             
     aTTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTGTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGc                                                                                                                                        <  1:2483653/148‑1 (MQ=255)
     aTTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGc                                                                                                                                        <  1:2621223/148‑1 (MQ=255)
     aTTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGc                                                                                                                                        <  2:2420405/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  2:774735/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  2:2134385/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  2:342258/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  2:2556669/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  2:2831450/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                       <  1:2213543/148‑1 (MQ=255)
       ggccgccGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCt                                                                                                                                      <  2:2909390/147‑1 (MQ=255)
           gccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGAt                                                                                                                                  >  1:2935088/1‑148 (MQ=255)
             cgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCa                                                                                                                                >  2:870835/1‑148 (MQ=255)
                tGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATaa                                                                                                                             <  2:3051445/148‑1 (MQ=255)
                 gCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATaaa                                                                                                                            >  2:498077/1‑148 (MQ=255)
                   aaaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGa                                                                                                                          <  1:2883579/148‑1 (MQ=255)
                    aaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGAc                                                                                                                         <  1:344835/148‑1 (MQ=255)
                     aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                        >  1:468802/1‑148 (MQ=255)
                     aaGCATCTAACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                        <  1:2088383/148‑1 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                       >  1:1428267/1‑148 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                       >  2:421145/1‑148 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                       <  1:452630/148‑1 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGa                                                                                                                       >  1:3277495/1‑148 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGa                                                                                                                       >  1:3277349/1‑148 (MQ=255)
                          tCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGc                                                                                                                   >  1:1424974/1‑148 (MQ=255)
                          tCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGc                                                                                                                   >  1:3192955/1‑148 (MQ=255)
                               cAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                              <  2:1949866/148‑1 (MQ=255)
                                aGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGAt                                                                                                             <  2:1206984/148‑1 (MQ=255)
                                 gTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGAtt                                                                                                            <  2:3224066/148‑1 (MQ=255)
                                  tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                           <  2:1679271/148‑1 (MQ=255)
                                  tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                           <  2:1102465/148‑1 (MQ=255)
                                   gAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTaa                                                                                                          <  1:3535828/148‑1 (MQ=255)
                                    aaGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAt                                                                                                         >  2:815081/1‑148 (MQ=255)
                                      gTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGt                                                                                                       >  2:2245432/1‑148 (MQ=255)
                                      gTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGt                                                                                                       >  2:2245364/1‑148 (MQ=255)
                                      gTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGt                                                                                                       >  2:2245356/1‑148 (MQ=255)
                                          aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCTGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                   >  1:3086080/1‑148 (MQ=255)
                                           tGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                  <  2:1292072/148‑1 (MQ=255)
                                           tGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                  <  2:1291659/148‑1 (MQ=255)
                                            gCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                                                                                                 >  1:1632621/1‑148 (MQ=255)
                                            gCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                                                                                                 >  2:1520744/1‑148 (MQ=255)
                                             ccGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                                                                                >  2:2833029/1‑148 (MQ=255)
                                             ccGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                                                                                <  1:3012394/148‑1 (MQ=255)
                                                cATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCg                                                                                             >  2:2835596/1‑148 (MQ=255)
                                                   ttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                          >  1:667554/1‑148 (MQ=255)
                                                   ttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                          >  2:1671883/1‑148 (MQ=255)
                                                   ttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                          >  1:2887015/1‑148 (MQ=255)
                                                      cgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAt                                                                                       <  2:191363/148‑1 (MQ=255)
                                                      cgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAt                                                                                       <  2:3496117/148‑1 (MQ=255)
                                                      cgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAt                                                                                       <  2:529971/148‑1 (MQ=255)
                                                            tctcGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGtt                                                                                 >  1:2424591/1‑148 (MQ=255)
                                                                 cATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGcc                                                                            >  1:3423568/1‑148 (MQ=255)
                                                                  aTGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCa                                                                           <  2:2153548/148‑1 (MQ=255)
                                                                     ccTGGCGGCGGTTTTAGCGCCGGGACGGGTGGCGGTATTTCTTTTGATGGTGAGTACCTGCGAATGGTTGATGAGCTCGACCGCATTGATCATAAAGACGATCGCCAGGATTATTGTCGAAACGCCGGATTATGTGGTTATGCCAttt                                                                        <  1:515895/148‑1 (MQ=255)
                                                                     ccTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAttt                                                                        <  1:516078/148‑1 (MQ=255)
                                                                            gcggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCg                                                                 <  1:1141981/148‑1 (MQ=255)
                                                                            gcggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCg                                                                 <  1:3404034/148‑1 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGCTTATGTGGTTATGCCATTTTCCGGCGt                                                                >  2:3094944/1‑148 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                >  1:1316684/1‑148 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                >  2:2822024/1‑148 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                <  2:1978612/148‑1 (MQ=255)
                                                                               gTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttt                                                              >  2:2489411/1‑148 (MQ=255)
                                                                                ttttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGtttt                                                             >  2:813077/1‑148 (MQ=255)
                                                                                 tttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                            <  1:866999/148‑1 (MQ=255)
                                                                                    aGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGt                                                         <  2:2935088/148‑1 (MQ=255)
                                                                                       gcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGtttt                                                      <  2:846419/148‑1 (MQ=255)
                                                                                         cGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTgg                                                    <  1:3008977/148‑1 (MQ=255)
                                                                                          gggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgc                                                   <  1:2078161/148‑1 (MQ=255)
                                                                                          gggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgc                                                   >  2:2982812/1‑148 (MQ=255)
                                                                                           ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgca                                                  >  2:3295789/1‑148 (MQ=255)
                                                                                                   ggCGGTAATACTTTTGAGGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGc                                          <  1:2399123/148‑1 (MQ=255)
                                                                                                    gCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGcc                                         <  2:1817502/148‑1 (MQ=255)
                                                                                                     cGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTGTTTCGTTTTGGCAGCAACCGCCa                                        >  1:1207776/1‑148 (MQ=255)
                                                                                                       gTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATa                                      >  2:2011998/1‑148 (MQ=255)
                                                                                                          aTACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAcc                                   >  1:2950201/1‑148 (MQ=255)
                                                                                                             cTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCg                                <  2:1434481/148‑1 (MQ=255)
                                                                                                             cTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACACCCATAACCCCg                                >  1:988190/1‑148 (MQ=255)
                                                                                                               tttGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGcc                              <  2:258153/148‑1 (MQ=255)
                                                                                                                   aTGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCa                          >  2:3475538/1‑148 (MQ=255)
                                                                                                                   aTGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCa                          >  1:974846/1‑148 (MQ=255)
                                                                                                                     ggTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATg                        >  2:696096/1‑148 (MQ=255)
                                                                                                                     ggTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATg                        >  1:3241045/1‑148 (MQ=255)
                                                                                                                      gTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGa                       <  1:1011690/148‑1 (MQ=255)
                                                                                                                       tGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGAt                      >  2:182150/1‑148 (MQ=255)
                                                                                                                         aGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAATCGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCa                    <  2:5619/148‑1 (MQ=255)
                                                                                                                         aGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCa                    <  2:667554/148‑1 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATTAGCTCGACCGCCTTGATCAAAAAGACGATCGCCAGGATAAATGTCGAAACGCAGGATTATCTGGTTATGCCATTTTCCGGCTTTTTTCGTTTTGGCATAAACAGCCAAAACCCACCCACCATGATCat                   >  1:3166016/1‑147 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                   <  1:3215228/148‑1 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                   <  1:525797/148‑1 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                   >  1:2967078/1‑148 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                   >  1:3164337/1‑148 (MQ=255)
                                                                                                                           tACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgc                  >  1:1574829/1‑148 (MQ=255)
                                                                                                                                gCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATa             >  2:2053294/1‑148 (MQ=255)
                                                                                                                                 cGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATag            >  2:1631667/1‑148 (MQ=255)
                                                                                                                                 cGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATag            >  1:1125796/1‑148 (MQ=255)
                                                                                                                                    aaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAgaga         <  2:692630/148‑1 (MQ=255)
                                                                                                                                       gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt      <  1:3475538/148‑1 (MQ=255)
                                                                                                                                       gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt      <  1:3418034/148‑1 (MQ=255)
                                                                                                                                       gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt      <  1:3415960/148‑1 (MQ=255)
                                                                                                                                        atgatgAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACttt     >  1:1392399/1‑148 (MQ=255)
                                                                                                                                           atgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTcc  <  1:589312/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCC  >  CP009273/4358775‑4359061

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: