Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,546,198 T→G 28.8% L163R (CTG→CGG)  amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,546,1980TG28.8% 386.1 / 308.3 339L163R (CTG→CGG) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base T (112/129);  new base G (46/52);  total (158/181)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 342 total aligned reads displayed.

ACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTGAAGAAGATTA  >  CP009273/2546049‑2546347
                                                                                                                                                     |                                                                                                                                                     
aCTTTACGATCGCGTTGTAATCGCCCATAATCATGGCGCAGTTCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCt                                                                                                                                                       <  2:2363644/150‑1 (MQ=255)
 cTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTg                                                                                                                                                      <  1:2474459/150‑1 (MQ=255)
   ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTc                                                                                                                                                    >  2:4954624/1‑150 (MQ=255)
     aCGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTg                                                                                                                                                  <  2:4838378/150‑1 (MQ=255)
      cGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTCACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGa                                                                                                                                                 >  2:1965232/1‑150 (MQ=255)
      cGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGa                                                                                                                                                 <  1:4104465/150‑1 (MQ=255)
      cGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGa                                                                                                                                                 >  1:879363/1‑150 (MQ=255)
      cGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGa                                                                                                                                                 >  2:1120532/1‑150 (MQ=255)
       gATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGaa                                                                                                                                                >  2:5250414/1‑150 (MQ=255)
        aTCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAAc                                                                                                                                               <  2:907804/150‑1 (MQ=255)
        aTCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAAc                                                                                                                                               <  2:2547672/150‑1 (MQ=255)
         tCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAAcg                                                                                                                                              <  1:4478749/150‑1 (MQ=255)
          cgcgTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAAcgc                                                                                                                                             >  2:2845224/1‑150 (MQ=255)
             gTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGaa                                                                                                                                          >  2:4582795/1‑150 (MQ=255)
             gTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGaa                                                                                                                                          >  2:2200852/1‑150 (MQ=255)
                  aaTCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATTCCTGTCTGATCGAGAAAACcg                                                                                                                                     <  2:2636790/150‑1 (MQ=255)
                     cGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCgcgc                                                                                                                                  >  2:5182281/1‑150 (MQ=255)
                      gCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAACACCGCGcc                                                                                                                                 >  2:3369518/1‑150 (MQ=255)
                           tAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGa                                                                                                                            <  2:3189492/150‑1 (MQ=255)
                             aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAg                                                                                                                          >  2:2877643/1‑150 (MQ=255)
                              aCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGt                                                                                                                         >  2:1229134/1‑150 (MQ=255)
                               cATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGtt                                                                                                                        >  1:5110820/1‑150 (MQ=255)
                                   gcgcAGATCTGTTTATGTCAAGTCATGCCGATGGCTTTACCAACCCGAACGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCg                                                                                                                    <  2:1523009/150‑1 (MQ=255)
                                       aGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaa                                                                                                                <  2:680418/150‑1 (MQ=255)
                                       aGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaa                                                                                                                >  2:358295/1‑150 (MQ=255)
                                        gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaaa                                                                                                               >  2:4691413/1‑150 (MQ=255)
                                        gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaaa                                                                                                               >  1:4202345/1‑150 (MQ=255)
                                        gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaaa                                                                                                               >  1:4841789/1‑150 (MQ=255)
                                         aTCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTaaaa                                                                                                              >  1:3214962/1‑150 (MQ=255)
                                           cTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAg                                                                                                            <  1:5292400/150‑1 (MQ=255)
                                           cTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAg                                                                                                            <  1:1247905/150‑1 (MQ=255)
                                             gTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAATCCGCGCCGATGAAGTTGCCGGTAAAAAGGc                                                                                                          >  2:3331989/1‑150 (MQ=255)
                                                      aaTTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACaa                                                                                                 >  2:1746448/1‑150 (MQ=255)
                                                      aaTTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACaa                                                                                                 >  1:2111345/1‑150 (MQ=255)
                                                       aTTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGGGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAg                                                                                                >  2:2135879/1‑150 (MQ=255)
                                                       aTTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAg                                                                                                >  1:3733095/1‑150 (MQ=255)
                                                        ttCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAgg                                                                                               <  2:1386688/150‑1 (MQ=255)
                                                         tCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGa                                                                                              <  1:316600/150‑1 (MQ=255)
                                                          cATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGAt                                                                                             <  1:3414341/150‑1 (MQ=255)
                                                          cATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGAt                                                                                             >  1:1702773/1‑150 (MQ=255)
                                                          cATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGAt                                                                                             >  2:3132366/1‑150 (MQ=255)
                                                             gCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCAc                                                                                          <  1:1074119/150‑1 (MQ=255)
                                                              ccGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCAcc                                                                                         <  1:3878196/150‑1 (MQ=255)
                                                                  tGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTAtt                                                                                     <  1:2754324/150‑1 (MQ=255)
                                                                  tGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTAtt                                                                                     <  2:1875080/150‑1 (MQ=255)
                                                                   ggCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTg                                                                                    <  1:134095/150‑1 (MQ=255)
                                                                   ggCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTg                                                                                    <  1:3097523/150‑1 (MQ=255)
                                                                    gCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGc                                                                                   >  2:2459101/1‑150 (MQ=255)
                                                                     cTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGca                                                                                  >  2:5220248/1‑150 (MQ=255)
                                                                          ccAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAg                                                                             <  1:261427/150‑1 (MQ=255)
                                                                           cAACACGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGt                                                                            >  1:2331004/1‑150 (MQ=255)
                                                                            aaCCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTg                                                                           <  1:3056460/150‑1 (MQ=255)
                                                                            aaCCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTg                                                                           >  2:740146/1‑150 (MQ=255)
                                                                               ccGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTg                                                                        <  1:4032580/150‑1 (MQ=255)
                                                                                cGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGt                                                                       >  2:2672943/1‑150 (MQ=255)
                                                                                  aaaGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGttt                                                                     <  2:1538479/150‑1 (MQ=255)
                                                                                   aaGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTg                                                                    >  1:5015541/1‑150 (MQ=255)
                                                                                   aaGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTg                                                                    <  2:302457/150‑1 (MQ=255)
                                                                                    aGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGa                                                                   >  2:5172675/1‑150 (MQ=255)
                                                                                    aGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGa                                                                   >  1:4338199/1‑150 (MQ=255)
                                                                                       tGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTCGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCt                                                                <  1:5078841/150‑1 (MQ=255)
                                                                                        gCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTg                                                               >  2:2760247/1‑150 (MQ=255)
                                                                                           ggTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTg                                                            >  1:1922827/1‑150 (MQ=255)
                                                                                            gTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGc                                                           <  2:3801361/150‑1 (MQ=255)
                                                                                            gTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGc                                                           <  2:610164/150‑1 (MQ=255)
                                                                                            gTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGc                                                           <  2:4271135/150‑1 (MQ=255)
                                                                                             tGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCa                                                          <  2:1818744/150‑1 (MQ=255)
                                                                                               cTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCaaa                                                        >  2:1888506/1‑150 (MQ=255)
                                                                                               cTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCaaa                                                        >  1:158277/1‑150 (MQ=255)
                                                                                                ttCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAAc                                                       <  2:3032186/150‑1 (MQ=255)
                                                                                                ttCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAAc                                                       <  1:3327529/150‑1 (MQ=255)
                                                                                                 tCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACa                                                      <  1:4783300/150‑1 (MQ=255)
                                                                                                 tCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACa                                                      <  2:4411221/150‑1 (MQ=255)
                                                                                                 tCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACa                                                      >  2:147017/1‑150 (MQ=255)
                                                                                                   ggTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGa                                                    >  2:2848125/1‑150 (MQ=255)
                                                                                                     tATTTGCCCTCTCTATCCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGTACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATa                                                  <  1:2543295/150‑1 (MQ=255)
                                                                                                      aTTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATAc                                                 <  2:397809/150‑1 (MQ=255)
                                                                                                        ttGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCa                                               <  2:3717567/150‑1 (MQ=255)
                                                                                                          gCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCAtt                                             <  1:3544410/150‑1 (MQ=255)
                                                                                                              tctctAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTaaaa                                         >  2:2253417/1‑150 (MQ=255)
                                                                                                               ctctAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTaaaaa                                        >  1:3826269/1‑150 (MQ=255)
                                                                                                                 ctAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATa                                      >  2:3653646/1‑150 (MQ=255)
                                                                                                                  tAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAg                                     >  2:4469714/1‑150 (MQ=255)
                                                                                                                      cGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTg                                 <  2:2266354/150‑1 (MQ=255)
                                                                                                                        tGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGAc                               >  2:4038623/1‑150 (MQ=255)
                                                                                                                        tGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGAc                               <  2:1839933/150‑1 (MQ=255)
                                                                                                                        tGGGGCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGAc                               <  2:4142173/150‑1 (MQ=255)
                                                                                                                         ggggCAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACg                              <  1:4734121/150‑1 (MQ=255)
                                                                                                                            gCAAGTAGCGAAATGGCGATATAACTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTTGTGCAAACAGATACCATTAAAAATAGTCTGACgctc                           <  1:2938858/150‑1 (MQ=255)
                                                                                                                             cAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcg                          >  1:4254916/1‑150 (MQ=255)
                                                                                                                             cAAGTAGCGCAATGGCGAAATACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcg                          <  2:4982683/150‑1 (MQ=255)
                                                                                                                                gTAGCGCACTGGCGAACTCCCGGTCTGAACGCGAAAACCGCGCCGATGAAGGGGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcggct                       <  1:705752/150‑1 (MQ=255)
                                                                                                                                    cgcAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCa                   <  2:2437402/150‑1 (MQ=255)
                                                                                                                                       aaTGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCatat                <  1:612298/150‑1 (MQ=255)
                                                                                                                                       aaTGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCatat                >  1:4668798/1‑150 (MQ=255)
                                                                                                                                        aTGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATAtt               <  1:2877643/150‑1 (MQ=255)
                                                                                                                                              aaaTACCGGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTgaag         >  2:5101168/1‑150 (MQ=255)
                                                                                                                                                aTACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTgaagaa       >  1:1460023/1‑150 (MQ=255)
                                                                                                                                                    cTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTGAAGAAGAtt   >  2:1166417/1‑150 (MQ=255)
                                                                                                                                                    cTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGATGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTGAAGAAGAtt   <  1:4944571/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCTCTCTAACCGTGGGGCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGGCTCGCATATTCTGAAGAAGATTA  >  CP009273/2546049‑2546347

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: