Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 231,378 C→A 89.7% R26S (CGC→AGC)  yafS → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273231,3780CA89.7% 834.0 / 88.1 316R26S (CGC→AGC) yafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (9/24);  new base A (148/135);  total (157/159)
Fisher's exact test for biased strand distribution p-value = 9.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

TAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCC  >  CP009273/231232‑231524
                                                                                                                                                  |                                                                                                                                                  
tAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGc                                                                                                                                                  >  2:37308/1‑149 (MQ=255)
tAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGc                                                                                                                                                  >  1:2261939/1‑149 (MQ=255)
  aTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCaa                                                                                                                                                >  2:1516434/1‑149 (MQ=255)
   ttGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAg                                                                                                                                               <  2:3000418/149‑1 (MQ=255)
      tcatcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcg                                                                                                                                            <  1:764932/149‑1 (MQ=255)
      acatcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcg                                                                                                                                            >  2:735116/2‑149 (MQ=255)
       catcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGgcgc                                                                                                                                           <  1:4363607/149‑1 (MQ=255)
       catcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGgcgc                                                                                                                                           <  1:3537327/149‑1 (MQ=255)
           aaaGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGa                                                                                                                                       <  1:3689931/149‑1 (MQ=255)
             aGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAgc                                                                                                                                     <  2:1246826/149‑1 (MQ=255)
                    aaTACTGTTAAGATTAATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCt                                                                                                                              >  2:885915/1‑149 (MQ=255)
                     aTACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTc                                                                                                                             >  2:20027/1‑149 (MQ=255)
                        cTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAc                                                                                                                          <  1:1098084/149‑1 (MQ=255)
                          gTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAccc                                                                                                                        <  1:105277/149‑1 (MQ=255)
                          gTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAccc                                                                                                                        >  2:2559400/1‑149 (MQ=255)
                            tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGt                                                                                                                      <  1:1196698/149‑1 (MQ=255)
                            tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGt                                                                                                                      <  1:1941053/149‑1 (MQ=255)
                             aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                     <  2:2334920/149‑1 (MQ=255)
                               gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                   <  2:3032335/149‑1 (MQ=255)
                                aTTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGtt                                                                                                                  >  2:4154754/1‑149 (MQ=255)
                                 ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTc                                                                                                                 <  2:411280/149‑1 (MQ=255)
                                 ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTc                                                                                                                 <  2:3835165/149‑1 (MQ=255)
                                 ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTc                                                                                                                 <  2:2805296/149‑1 (MQ=255)
                                   cATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCAc                                                                                                               <  2:4716236/149‑1 (MQ=255)
                                    aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACt                                                                                                              >  2:4617018/1‑149 (MQ=255)
                                    aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACt                                                                                                              >  1:3623146/1‑149 (MQ=255)
                                       gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                           >  2:2712948/1‑149 (MQ=255)
                                       gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                           >  1:337924/1‑149 (MQ=255)
                                        aTTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaaa                                                                                                          >  2:2447702/1‑149 (MQ=255)
                                         ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAAt                                                                                                         <  2:777053/149‑1 (MQ=255)
                                         ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAAt                                                                                                         <  2:1370887/149‑1 (MQ=255)
                                         ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAAt                                                                                                         <  1:682241/149‑1 (MQ=255)
                                           aCCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGt                                                                                                       <  2:2060956/149‑1 (MQ=255)
                                            ccTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTa                                                                                                      >  2:878470/1‑149 (MQ=255)
                                            ccTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTa                                                                                                      <  2:4339549/149‑1 (MQ=255)
                                             ctctCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTAt                                                                                                     <  1:3566198/149‑1 (MQ=255)
                                             ctctCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTAt                                                                                                     <  1:837656/149‑1 (MQ=255)
                                              tctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATg                                                                                                    >  1:3115076/1‑149 (MQ=255)
                                                tcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                  >  1:165441/1‑149 (MQ=255)
                                                tcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                  <  2:710558/149‑1 (MQ=255)
                                                tcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTGCACTAAACTGTATGGt                                                                                                  >  1:2521132/1‑149 (MQ=255)
                                                   gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGtttt                                                                                               >  2:980362/1‑149 (MQ=255)
                                                      tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAACATGTATGGTTTTCAt                                                                                            >  2:4029861/1‑149 (MQ=255)
                                                       gAAACGGAATGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATc                                                                                           >  1:26184/1‑149 (MQ=255)
                                                          aCGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGc                                                                                        <  1:3205152/149‑1 (MQ=255)
                                                           cGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt                                                                                       <  1:4480015/149‑1 (MQ=255)
                                                              aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa                                                                                    <  1:1354565/149‑1 (MQ=255)
                                                              aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa                                                                                    <  1:46055/149‑1 (MQ=255)
                                                               aGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAg                                                                                   >  2:3777101/1‑149 (MQ=255)
                                                                 gTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAt                                                                                 <  1:2002316/149‑1 (MQ=255)
                                                                  ttGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                                >  1:1585985/1‑149 (MQ=255)
                                                                  ttGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                                >  2:4527314/1‑149 (MQ=255)
                                                                     tgatgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGc                                                                             >  1:3978812/1‑149 (MQ=255)
                                                                         gaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAtt                                                                         <  1:964721/149‑1 (MQ=255)
                                                                           aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTa                                                                       >  1:2715794/1‑149 (MQ=255)
                                                                           aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTa                                                                       >  1:3864905/1‑149 (MQ=255)
                                                                             ccGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAg                                                                     >  2:3552528/1‑149 (MQ=255)
                                                                              cGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgc                                                                    <  1:499627/149‑1 (MQ=255)
                                                                              cGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgc                                                                    <  2:251567/149‑1 (MQ=255)
                                                                               ggCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                                   >  1:4126237/1‑149 (MQ=255)
                                                                                   agagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                               <  2:1802542/149‑1 (MQ=255)
                                                                                   agagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTCAGATTGGCAATTTAAGCGCAGa                                                               >  2:3071619/1‑149 (MQ=255)
                                                                                       taccTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGTCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATc                                                           >  1:206416/3‑149 (MQ=255)
                                                                                       tCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATc                                                           >  1:2482336/1‑149 (MQ=255)
                                                                                        cccTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCa                                                          >  1:610536/1‑149 (MQ=255)
                                                                                         ccTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCaa                                                         <  1:2031692/149‑1 (MQ=255)
                                                                                           tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt                                                       <  2:2662135/149‑1 (MQ=255)
                                                                                             aagcTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTGTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGc                                                     >  2:1402126/4‑149 (MQ=255)
                                                                                             aaaCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGc                                                     >  2:2956169/1‑149 (MQ=255)
                                                                                                 tGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAg                                                 <  1:681047/149‑1 (MQ=255)
                                                                                                 tGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAg                                                 >  1:2392529/1‑149 (MQ=255)
                                                                                                   tCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCg                                               <  1:1958917/149‑1 (MQ=255)
                                                                                                      tGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                                            <  1:2179281/149‑1 (MQ=255)
                                                                                                       ggCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcg                                           >  1:131096/1‑149 (MQ=255)
                                                                                                          tCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                        >  1:947244/1‑149 (MQ=255)
                                                                                                            cTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGtt                                      >  1:4666492/1‑149 (MQ=255)
                                                                                                            cTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGtt                                      >  1:3249527/1‑149 (MQ=255)
                                                                                                             tGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                                     <  1:2877229/149‑1 (MQ=255)
                                                                                                                ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                                  <  2:566044/149‑1 (MQ=255)
                                                                                                                ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                                  <  2:1943691/149‑1 (MQ=255)
                                                                                                                ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGAGAAGCGTGCGCGGTTtctc                                  >  1:4074151/1‑149 (MQ=255)
                                                                                                                 tGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtca                                 <  1:4430310/149‑1 (MQ=255)
                                                                                                                   cTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatc                               <  2:337924/149‑1 (MQ=255)
                                                                                                                         cGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGa                         <  1:1369294/149‑1 (MQ=255)
                                                                                                                          gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGaa                        <  2:2302604/149‑1 (MQ=255)
                                                                                                                            tttGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATg                      >  2:4202443/1‑149 (MQ=255)
                                                                                                                              tGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGtt                    <  2:3194704/149‑1 (MQ=255)
                                                                                                                               gCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGttt                   >  2:3252470/1‑149 (MQ=255)
                                                                                                                                  cTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGACTGTTTCTg                >  2:3193463/1‑149 (MQ=255)
                                                                                                                                   tGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgc               <  2:2906163/149‑1 (MQ=255)
                                                                                                                                    ggggAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgcg              >  2:3316700/1‑149 (MQ=255)
                                                                                                                                       gAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCaa           >  1:4512705/1‑149 (MQ=255)
                                                                                                                                         aaGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAgg         <  1:3593673/149‑1 (MQ=255)
                                                                                                                                           gCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa       >  1:4249008/1‑149 (MQ=255)
                                                                                                                                           gCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa       <  2:174315/149‑1 (MQ=255)
                                                                                                                                           gCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa       <  2:835598/149‑1 (MQ=255)
                                                                                                                                            cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAt      <  1:37939/149‑1 (MQ=255)
                                                                                                                                            cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAt      <  2:1025445/149‑1 (MQ=255)
                                                                                                                                            cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAt      >  1:1464022/1‑149 (MQ=255)
                                                                                                                                                aTCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGccc  >  1:4727239/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCC  >  CP009273/231232‑231524

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: