Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,546,672 G→T 26.0% intergenic (‑535/‑138) yjiC ← / → iraD uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,546,6720GT26.0% 664.2 / 179.7 308intergenic (‑535/‑138)yjiC/iraDuncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base G (117/110);  new base T (46/35);  total (163/145)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

GTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCG  >  CP009273/4546523‑4546817
                                                                                                                                                      |                                                                                                                                                 
gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTg                                                                                                                                                   <  1:2083693/150‑1 (MQ=255)
  gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGAc                                                                                                                                                 >  2:3854510/1‑150 (MQ=255)
    cAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTTTACTTTGCTGACAg                                                                                                                                               >  2:1615896/1‑150 (MQ=255)
    cAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATAATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAg                                                                                                                                               >  2:3114184/1‑150 (MQ=255)
     aaGTATTTTCATCTACGTAACTATTAGGAATATTAAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAgg                                                                                                                                              <  1:1761628/150‑1 (MQ=255)
      aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGa                                                                                                                                             <  1:1482069/150‑1 (MQ=255)
      aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGa                                                                                                                                             >  1:4745352/1‑150 (MQ=255)
       gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAt                                                                                                                                            >  2:368871/1‑150 (MQ=255)
       gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTCAGTAGATATTCATATGCAAGTAAGGTTTTATACTTTGATGACAGGAt                                                                                                                                            >  1:2283159/1‑150 (MQ=255)
          ttttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATAATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAATGTTTTATACTTTGCTGACAGGATTCa                                                                                                                                         >  1:942057/1‑150 (MQ=255)
                     gTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCa                                                                                                                              >  2:1705458/1‑150 (MQ=255)
                      tAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAg                                                                                                                             >  2:4587805/1‑150 (MQ=255)
                         cTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGATTTTGATTATTTTTATAAAATAATTACTTTTTTACTGTAAATTAATTAGATATTAATATGCAAGTAAGGTTTTATAATTTGCTGACAGGATTCAGGCCTGTCTCAGACt                                                                                                                          >  1:1467668/1‑150 (MQ=255)
                           aTTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATTATTTTTATAAAATAATTACATTTTTAATGAAAATTAAGTAGATATTAATATGCAAGTTAGGTTTTTTAAATTTCTTAAAGGATTCAGGCCTTTATCAGACTg                                                                                                                         >  1:1022994/1‑150 (MQ=255)
                             tCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACa                                                                                                                      <  1:1402246/150‑1 (MQ=255)
                              cGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACAt                                                                                                                     <  1:4797036/150‑1 (MQ=255)
                                gAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATAATTTGCTGACAGGATTCAGGCCTGTCTCATACTGACATgg                                                                                                                   >  2:4904962/1‑150 (MQ=255)
                                   tatTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATg                                                                                                                <  2:2150180/150‑1 (MQ=255)
                                    atTCAGAATATACTAGCGGTTTTTTTTTGTGTTTTATGATTTTTATAAAATAAATAATTTTTTACAGAAAAGTAAGTATATATTAATAGGCAAATAATGTTTTATACATTTCTGACAGGATTAATGCCAGTCTAAGACTGACATTGAAGt                                                                                                               >  1:2303928/1‑150 (MQ=255)
                                         gAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATg                                                                                                           <  2:2618675/150‑1 (MQ=255)
                                           atatACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATTACTTTTTTACTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATAATTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATTTAATGAAc                                                                                                        >  2:2418639/1‑150 (MQ=255)
                                              taCTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATAAAATCATTACTTTTTTCCTGAAAATTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTTATTAAAGGATTCATGCCTGTCTCAGACAGACATGGATGTAATGaaaaaa                                                                                                     >  1:2868244/1‑150 (MQ=255)
                                                cTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAg                                                                                                   <  2:3029102/150‑1 (MQ=255)
                                                cTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAg                                                                                                   >  2:4534479/1‑150 (MQ=255)
                                                cTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAg                                                                                                   <  2:1176673/150‑1 (MQ=255)
                                                 tAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTANTGAACAAAAgg                                                                                                  <  1:2762856/150‑1 (MQ=255)
                                                  aGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTATATATTCATAGGCAAGTAAGGTTTTATAATTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAggg                                                                                                 >  1:22690/1‑150 (MQ=255)
                                                  aGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAggg                                                                                                 <  2:64154/150‑1 (MQ=255)
                                                    cGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACATTGCTTACAGGATTAAGGCCTGTATCAGACTGACATGGATGTAATGAAAAAAAGGGaa                                                                                               >  1:2468118/1‑150 (MQ=255)
                                                     ggTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAAt                                                                                              <  2:5146567/150‑1 (MQ=255)
                                                      gTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATg                                                                                             <  2:842873/150‑1 (MQ=255)
                                                       ggttttttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATg                                                                                             >  2:3802265/3‑150 (MQ=255)
                                                        tttttttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATgg                                                                                            <  2:3547533/150‑1 (MQ=255)
                                                          tttttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCt                                                                                          <  1:4198267/150‑1 (MQ=255)
                                                          tttttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCt                                                                                          <  1:724085/150‑1 (MQ=255)
                                                            tttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTAt                                                                                        <  1:3535958/150‑1 (MQ=255)
                                                              tttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATgg                                                                                      <  1:36888/150‑1 (MQ=255)
                                                                tGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTAGACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGaa                                                                                    <  2:4743117/150‑1 (MQ=255)
                                                                  aGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGaaaa                                                                                  >  2:4666995/1‑150 (MQ=255)
                                                                   gTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAAt                                                                                 <  1:2153101/150‑1 (MQ=255)
                                                                    tttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATg                                                                                >  2:544943/1‑150 (MQ=255)
                                                                    tttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATg                                                                                >  2:3677577/1‑150 (MQ=255)
                                                                    tttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATg                                                                                >  1:1241149/1‑150 (MQ=255)
                                                                     ttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGa                                                                               <  1:2365271/150‑1 (MQ=255)
                                                                      tgatgaTTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAg                                                                              >  2:2935760/1‑150 (MQ=255)
                                                                       gatgatTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGc                                                                             <  2:882252/150‑1 (MQ=255)
                                                                       gatgatTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGc                                                                             <  2:2563293/150‑1 (MQ=255)
                                                                       gatgatTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGc                                                                             <  1:2577881/150‑1 (MQ=255)
                                                                           atTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATTGCTATGGAAAATGAGCATCa                                                                         <  1:3870720/150‑1 (MQ=255)
                                                                                tATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACa                                                                    <  1:4543387/150‑1 (MQ=255)
                                                                                tATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACa                                                                    >  2:4050964/1‑150 (MQ=255)
                                                                                  tCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAgt                                                                  <  2:1328451/150‑1 (MQ=255)
                                                                                   cAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAgtg                                                                 <  2:530145/150‑1 (MQ=255)
                                                                                    aaaTCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAgtgg                                                                >  1:1688959/1‑150 (MQ=255)
                                                                                     aaTCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAgtggt                                                               <  2:1974243/150‑1 (MQ=255)
                                                                                      aTCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAgtggtg                                                              >  1:1367671/1‑150 (MQ=255)
                                                                                             cTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGgtg                                                       >  2:11092/1‑150 (MQ=255)
                                                                                             cTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGgtg                                                       <  2:3845650/150‑1 (MQ=255)
                                                                                                ttgtCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTc                                                    <  1:5113079/147‑1 (MQ=255)
                                                                                                 tttCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCa                                                   <  1:2646927/150‑1 (MQ=255)
                                                                                                 tttCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCa                                                   <  2:4437192/150‑1 (MQ=255)
                                                                                                   tCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAgg                                                 <  1:776624/150‑1 (MQ=255)
                                                                                                    ccTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAggg                                                <  1:4771519/150‑1 (MQ=255)
                                                                                                     cTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGc                                               <  1:2936341/150‑1 (MQ=255)
                                                                                                        aaaaGTCAGTAGATATTCATAGGCAAGTAATGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAg                                            <  1:518989/150‑1 (MQ=255)
                                                                                                          aaGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGcc                                          <  1:3625970/150‑1 (MQ=255)
                                                                                                            gTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGc                                        >  2:5234935/1‑150 (MQ=255)
                                                                                                             tCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCa                                       >  2:2465044/1‑150 (MQ=255)
                                                                                                               aGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCata                                     >  2:810375/1‑150 (MQ=255)
                                                                                                               aGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTGATGAACAAAAGGGAATGGCGATGGAAAATGAGCCTCAATACAGGGGTGCCCGGTGTTCAGGGCGAGCCGctgc                                     >  2:5113661/1‑147 (MQ=255)
                                                                                                                 tAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATg                                   <  1:1890665/150‑1 (MQ=255)
                                                                                                                  aGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGt                                  >  2:3385417/1‑150 (MQ=255)
                                                                                                                  aGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGt                                  >  1:1809893/1‑150 (MQ=255)
                                                                                                                    atatTCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTg                                >  2:4485614/1‑150 (MQ=255)
                                                                                                                        tCATAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTaa                            <  2:3284845/150‑1 (MQ=255)
                                                                                                                         cATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTaaa                           <  1:3223153/150‑1 (MQ=255)
                                                                                                                          aTAGGCAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCACGCCGCATATGTTGCTAAAc                          >  2:2388006/1‑150 (MQ=255)
                                                                                                                           gaGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCCGATGTTGCTAAACg                         <  2:1532798/149‑1 (MQ=255)
                                                                                                                               cAAGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAg                     >  1:5277870/1‑150 (MQ=255)
                                                                                                                               cAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAg                     <  2:2009193/150‑1 (MQ=255)
                                                                                                                                aaGTAAGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAgg                    >  1:4396213/1‑150 (MQ=255)
                                                                                                                                aaGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAgg                    >  1:2994192/1‑150 (MQ=255)
                                                                                                                                 aGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGa                   >  2:684164/1‑150 (MQ=255)
                                                                                                                                 aGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGa                   <  2:4796024/150‑1 (MQ=255)
                                                                                                                                  gTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAg                  <  2:2761572/150‑1 (MQ=255)
                                                                                                                                    aaGGTTTTATACTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTg                >  1:4554533/1‑150 (MQ=255)
                                                                                                                                    aaGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTg                >  2:192233/1‑150 (MQ=255)
                                                                                                                                      ggTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTgcg              >  2:4999575/1‑150 (MQ=255)
                                                                                                                                           tataCTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCaaaa         >  1:3042296/1‑150 (MQ=255)
                                                                                                                                           tataCTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCaaaa         >  2:2348436/1‑150 (MQ=255)
                                                                                                                                            ataCTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAat        >  1:465475/1‑150 (MQ=255)
                                                                                                                                            ataCTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAat        >  2:4824266/1‑150 (MQ=255)
                                                                                                                                            ataCTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAat        <  1:3253864/150‑1 (MQ=255)
                                                                                                                                             taCTTTGCTTACAGGATTCATGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAatg       >  1:154221/1‑150 (MQ=255)
                                                                                                                                             taCTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAatg       >  1:1321852/1‑150 (MQ=255)
                                                                                                                                               cTTTGCTTACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAatgat     >  2:5166415/1‑150 (MQ=255)
                                                                                                                                               cTTTGCTGACAGGATTCAGTCCTGTCTAAGACTGACATGGTTGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCTAAatgat     <  1:2503480/150‑1 (MQ=255)
                                                                                                                                               cTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAatgat     >  1:2248243/1‑150 (MQ=255)
                                                                                                                                               cTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAatgat     <  2:942057/150‑1 (MQ=255)
                                                                                                                                                  tGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCg  <  1:4999930/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCG  >  CP009273/4546523‑4546817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: