Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 918,501 G→A 18.8% D45N (GAC→AAC)  clpS → regulatory protein for ClpA substrate specificity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273918,5010GA18.8% 500.8 / 102.6 202D45N (GAC→AAC) clpSregulatory protein for ClpA substrate specificity
Reads supporting (aligned to +/- strand):  ref base G (81/83);  new base A (17/21);  total (98/104)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 206 total aligned reads displayed.

GAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGA  >  CP009273/918352‑918649
                                                                                                                                                       |                                                                                                                                                    
gAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTg                                                                                                                                                      <  2:292989/150‑1 (MQ=255)
gAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTg                                                                                                                                                      >  1:1236497/1‑150 (MQ=255)
     ttatAACTGACATTGGTAAAAACTAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgtg                                                                                                                                                 <  1:1439846/148‑1 (MQ=255)
     tgatAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgtg                                                                                                                                                 >  1:3449309/1‑150 (MQ=255)
      gatAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgtgt                                                                                                                                                <  2:1657938/150‑1 (MQ=255)
      gagaaCTGACCCGGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgtgt                                                                                                                                                <  2:1656516/147‑1 (MQ=255)
       atACCTGACATTGGTTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGtt                                                                                                                                               <  2:4676253/150‑1 (MQ=255)
       atAACTGACAATGGTTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGtt                                                                                                                                               <  2:4676171/150‑1 (MQ=255)
        tAACTGACATTGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTa                                                                                                                                              <  2:880247/150‑1 (MQ=255)
        tAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTa                                                                                                                                              >  1:4446838/1‑150 (MQ=255)
        tAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGATATTGACGTGTTa                                                                                                                                              <  2:1573746/150‑1 (MQ=255)
        tAACTGAAATTGGGTAAACCGCACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTa                                                                                                                                              <  1:2926451/150‑1 (MQ=255)
        tAACTCAAAATGGTAAACACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTa                                                                                                                                              <  1:4196435/150‑1 (MQ=255)
        aaaCTGACATTGGGAAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTa                                                                                                                                              <  1:199185/149‑1 (MQ=255)
                aaagggTAAAACCCAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAt                                                                                                                                       <  1:3721215/147‑1 (MQ=255)
                 aTGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                                                                                                                                     >  2:1239857/1‑150 (MQ=255)
                  tGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttct                                                                                                                                    <  2:815286/150‑1 (MQ=255)
                    ggTAAACACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctt                                                                                                                                   <  1:4116528/150‑1 (MQ=255)
                     gTAAACCGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctttt                                                                                                                                 <  2:2675937/150‑1 (MQ=255)
                     gTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctttt                                                                                                                                 <  2:2675839/150‑1 (MQ=255)
                      tAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttc                                                                                                                                >  2:3965221/1‑150 (MQ=255)
                       aaaaCGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttct                                                                                                                               >  1:3580674/1‑150 (MQ=255)
                       aaaaCGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTAAACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttct                                                                                                                               >  2:2328680/1‑150 (MQ=255)
                       aaaaCCAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAAttcttttct                                                                                                                               <  1:2498476/150‑1 (MQ=255)
                        aaaCGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTCTTATCGTGTTACAAAAAttctttgctt                                                                                                                              >  2:3335620/1‑150 (MQ=255)
                        aaaCGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttctt                                                                                                                              <  2:4211062/150‑1 (MQ=255)
                             cgaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatg                                                                                                                           >  1:2642212/1‑150 (MQ=255)
                             cgaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatg                                                                                                                           >  1:4634226/1‑150 (MQ=255)
                              gaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatga                                                                                                                          <  1:1588917/150‑1 (MQ=255)
                              gaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatga                                                                                                                          >  2:720613/1‑150 (MQ=255)
                              gaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTatga                                                                                                                          <  2:3254028/150‑1 (MQ=255)
                               aacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgat                                                                                                                         <  1:1142746/150‑1 (MQ=255)
                                acgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTatgatg                                                                                                                        <  2:2907266/150‑1 (MQ=255)
                                 cgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGt                                                                                                                       <  2:2635367/150‑1 (MQ=255)
                                  gaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTa                                                                                                                      <  1:3846741/150‑1 (MQ=255)
                                    ctggctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGa                                                                                                                    <  1:2045770/150‑1 (MQ=255)
                                    ctggctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGa                                                                                                                    <  1:1553121/150‑1 (MQ=255)
                                    ctggctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGa                                                                                                                    >  2:3491362/1‑150 (MQ=255)
                                      ggctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACTTGTTACAAAAATTCTTTTCTTATGATGTAGAAc                                                                                                                  <  1:3046686/150‑1 (MQ=255)
                                      ggctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAAc                                                                                                                  >  1:51263/1‑150 (MQ=255)
                                        ctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGt                                                                                                                >  1:2614268/1‑150 (MQ=255)
                                           gACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCa                                                                                                             <  1:750588/150‑1 (MQ=255)
                                                 gATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCaa                                                                                                       >  1:3813227/1‑150 (MQ=255)
                                                 gATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCaa                                                                                                       >  2:2790177/1‑150 (MQ=255)
                                                  aTCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAAt                                                                                                      >  1:2097451/1‑150 (MQ=255)
                                                     aaCTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGa                                                                                                   >  1:1206599/1‑150 (MQ=255)
                                                           cGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCg                                                                                             >  1:2484605/1‑150 (MQ=255)
                                                             gaagaaAAAGTTCGCGACGCGCTCACACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAGTCTTTTCGTATGATGTAGAACGTGCAACGCAATTGATGCTCGCt                                                                                           <  2:2985446/150‑1 (MQ=255)
                                                             gaagaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCt                                                                                           >  1:631980/1‑150 (MQ=255)
                                                              aagaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTg                                                                                          >  2:2711999/1‑150 (MQ=255)
                                                              aagaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTg                                                                                          >  2:2786299/1‑150 (MQ=255)
                                                              aagaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTg                                                                                          >  1:4343402/1‑150 (MQ=255)
                                                                gaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGtt                                                                                        >  1:108922/1‑150 (MQ=255)
                                                                    aaGTTCGCGACGCGCTAAAACCGCCATCTATGTATACAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACt                                                                                    <  2:4462288/150‑1 (MQ=255)
                                                                     aGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTa                                                                                   >  1:3403007/1‑150 (MQ=255)
                                                                       ttCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTAcc                                                                                 <  1:4608547/150‑1 (MQ=255)
                                                                       ttCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTAcc                                                                                 <  1:440204/150‑1 (MQ=255)
                                                                       ttCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTAcc                                                                                 >  1:4273507/1‑150 (MQ=255)
                                                                          gcgACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAgg                                                                              <  2:3580674/150‑1 (MQ=255)
                                                                             aCGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGa                                                                           <  1:781448/150‑1 (MQ=255)
                                                                                 gcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGc                                                                       >  2:2498125/1‑150 (MQ=255)
                                                                                  cTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGcc                                                                      >  1:1169890/1‑150 (MQ=255)
                                                                                   tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     >  2:1394703/1‑150 (MQ=255)
                                                                                   tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     >  1:3301226/1‑150 (MQ=255)
                                                                                          gCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCgg                                                              >  2:477253/1‑150 (MQ=255)
                                                                                          gCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCgg                                                              >  1:3019791/1‑150 (MQ=255)
                                                                                          gCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGATATTGACGTGATACAAAAATTCATTACTTATGATGAAGAACGTGCAAAGCAATTGATGCTCGCTGTTCACTACAAGGGGAAGCCCATATGCgg                                                              >  2:4007307/1‑150 (MQ=255)
                                                                                           ccATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGa                                                             >  2:1178256/1‑150 (MQ=255)
                                                                                             aTCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGt                                                           <  1:1186111/150‑1 (MQ=255)
                                                                                               cTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCt                                                         <  2:316480/150‑1 (MQ=255)
                                                                                                tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCtt                                                        >  2:2821134/1‑150 (MQ=255)
                                                                                                tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCtt                                                        <  1:1087554/150‑1 (MQ=255)
                                                                                                 aTGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCttt                                                       <  2:2958250/150‑1 (MQ=255)
                                                                                                 aTGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCttt                                                       <  2:1818170/150‑1 (MQ=255)
                                                                                                  tGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTa                                                      <  1:2300454/150‑1 (MQ=255)
                                                                                                   gTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTAc                                                     >  1:3216054/1‑150 (MQ=255)
                                                                                                     ataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTAccg                                                   <  1:350784/150‑1 (MQ=255)
                                                                                                         aGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGa                                               <  2:4617521/150‑1 (MQ=255)
                                                                                                             atatTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGtt                                           <  2:4038579/150‑1 (MQ=255)
                                                                                                                ttAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGACGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCa                                        >  2:3038733/1‑150 (MQ=255)
                                                                                                                   gTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACAAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGaa                                     >  2:4697896/1‑150 (MQ=255)
                                                                                                                    tCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGaaa                                    >  2:1976840/1‑150 (MQ=255)
                                                                                                                     cAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAAc                                   <  2:3449309/150‑1 (MQ=255)
                                                                                                                     cAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAAc                                   <  1:3992537/150‑1 (MQ=255)
                                                                                                                      aaTGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAAcc                                  >  1:2823250/1‑150 (MQ=255)
                                                                                                                       atgatgATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCa                                 >  1:2287690/1‑150 (MQ=255)
                                                                                                                        tgatgaTTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCaa                                >  1:3069413/1‑150 (MQ=255)
                                                                                                                               tACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCg                         <  2:2202019/150‑1 (MQ=255)
                                                                                                                               tACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCg                         <  2:4652246/150‑1 (MQ=255)
                                                                                                                                 cacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGAt                       >  1:479552/1‑150 (MQ=255)
                                                                                                                                  acTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTACAGAAACCAAAGTGGCGATg                      >  1:2566894/1‑150 (MQ=255)
                                                                                                                                   cTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATgg                     <  2:4663824/150‑1 (MQ=255)
                                                                                                                                     ccGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTg                   <  2:108922/150‑1 (MQ=255)
                                                                                                                                      cGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGa                  <  1:2761258/150‑1 (MQ=255)
                                                                                                                                       gATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGaa                 >  1:1291972/1‑150 (MQ=255)
                                                                                                                                          ggAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACaa              <  2:1315819/150‑1 (MQ=255)
                                                                                                                                          ggAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACaa              >  1:1006315/1‑150 (MQ=255)
                                                                                                                                               ttGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTAcg         <  1:3633902/150‑1 (MQ=255)
                                                                                                                                                tGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTAcgc        >  1:2850679/1‑150 (MQ=255)
                                                                                                                                                      tGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGga  <  1:2417787/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGA  >  CP009273/918352‑918649

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: