Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,887,963 T→A 10.0% S241S (TCT→TCA yieE → phosphopantetheinyl transferase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,887,9630TA10.0% 517.1 / 42.8 180S241S (TCT→TCAyieEphosphopantetheinyl transferase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (65/95);  new base A (10/8);  total (75/105)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 185 total aligned reads displayed.

CTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGC  >  CP009273/3887817‑3888109
                                                                                                                                                  |                                                                                                                                                  
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACGATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  2:1827042/148‑1 (MQ=255)
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  1:291100/148‑1 (MQ=255)
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  2:1190872/148‑1 (MQ=255)
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  2:2934407/148‑1 (MQ=255)
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  2:3158646/148‑1 (MQ=255)
cTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTa                                                                                                                                                   <  2:3158688/148‑1 (MQ=255)
    tgtgGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGTATGATGCGTTTTGCCCAACTCTCTACCGt                                                                                                                                               <  2:904339/148‑1 (MQ=255)
      tgGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGa                                                                                                                                             <  1:762566/148‑1 (MQ=255)
      tgGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGa                                                                                                                                             >  1:236950/1‑148 (MQ=255)
       ggTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGaa                                                                                                                                            <  1:3309571/148‑1 (MQ=255)
         tCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAgg                                                                                                                                          <  2:126587/148‑1 (MQ=255)
          ccGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGc                                                                                                                                         <  2:722673/148‑1 (MQ=255)
              ggCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCtttt                                                                                                                                     <  1:1060088/148‑1 (MQ=255)
               gCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCttttt                                                                                                                                    <  2:1548035/148‑1 (MQ=255)
                cGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTc                                                                                                                                   <  2:1872860/148‑1 (MQ=255)
                 ggTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCg                                                                                                                                  <  2:2884187/148‑1 (MQ=255)
                 ggTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCg                                                                                                                                  <  2:2885421/148‑1 (MQ=255)
                   tCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGcc                                                                                                                                >  2:2446195/1‑148 (MQ=255)
                    cACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCa                                                                                                                               >  1:1721053/1‑148 (MQ=255)
                              aTTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGa                                                                                                                     <  1:214410/148‑1 (MQ=255)
                              aTTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGa                                                                                                                     <  1:214449/148‑1 (MQ=255)
                               ttGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGAc                                                                                                                    >  2:333619/1‑148 (MQ=255)
                                tGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACa                                                                                                                   <  1:1139716/148‑1 (MQ=255)
                                  aaaaGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAAc                                                                                                                 <  1:1113702/148‑1 (MQ=255)
                                    aaGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACca                                                                                                               >  1:1649576/1‑148 (MQ=255)
                                      gCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACcaca                                                                                                             >  2:1683172/1‑148 (MQ=255)
                                      gCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACcaca                                                                                                             >  2:420729/1‑148 (MQ=255)
                                          aGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAg                                                                                                         <  2:2942159/148‑1 (MQ=255)
                                            tgtCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGtc                                                                                                       >  1:967979/1‑148 (MQ=255)
                                              tCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGtcat                                                                                                     <  1:1312905/148‑1 (MQ=255)
                                               cTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGtcatc                                                                                                    <  2:1158309/148‑1 (MQ=255)
                                                 gggAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTcatcat                                                                                                  <  1:95942/148‑1 (MQ=255)
                                                 gggAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTcatcat                                                                                                  >  1:101300/1‑148 (MQ=255)
                                                   gAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGt                                                                                                >  1:1865913/1‑148 (MQ=255)
                                                   gAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGt                                                                                                <  1:518872/148‑1 (MQ=255)
                                                     gTTAGATGTCAAACACGTCTGGAAAAGCCTGCCGTATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCTAATTGATAGACAACCACAGGAGTCATCATGTCt                                                                                              <  1:1278658/148‑1 (MQ=255)
                                                     gTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCt                                                                                              <  2:2527561/148‑1 (MQ=255)
                                                     gTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCt                                                                                              <  1:1281842/148‑1 (MQ=255)
                                                      ttAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTg                                                                                             >  2:2672514/1‑148 (MQ=255)
                                                          aTGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGaaaa                                                                                         >  2:308456/1‑148 (MQ=255)
                                                            ggCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAAt                                                                                       >  1:55389/1‑148 (MQ=255)
                                                                aaCACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCa                                                                                   >  1:2959389/1‑148 (MQ=255)
                                                                  cacGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAgg                                                                                 <  1:2544606/148‑1 (MQ=255)
                                                                      gCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAggtggt                                                                             <  2:675892/148‑1 (MQ=255)
                                                                       cTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgtt                                                                            <  2:1343272/148‑1 (MQ=255)
                                                                        tGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgtta                                                                           >  1:2784993/1‑148 (MQ=255)
                                                                        tGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgtta                                                                           >  1:2013441/1‑148 (MQ=255)
                                                                         ggAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgttac                                                                          >  1:428433/1‑148 (MQ=255)
                                                                           aaaaGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgttacgt                                                                        <  1:103034/148‑1 (MQ=255)
                                                                              aGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGgttacgttac                                                                     >  2:3270812/1‑148 (MQ=255)
                                                                               gCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACt                                                                    >  1:532442/1‑148 (MQ=255)
                                                                                ccTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTg                                                                   >  1:162794/1‑148 (MQ=255)
                                                                                ccTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTg                                                                   >  2:1986470/1‑148 (MQ=255)
                                                                                  tGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTggg                                                                 <  1:2712448/148‑1 (MQ=255)
                                                                                   gCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTgggg                                                                <  1:155484/148‑1 (MQ=255)
                                                                                    ccGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGa                                                               <  1:2216628/148‑1 (MQ=255)
                                                                                      ggATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGc                                                             <  2:2720166/148‑1 (MQ=255)
                                                                                      ggATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGc                                                             <  2:2836683/148‑1 (MQ=255)
                                                                                      ggATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGc                                                             <  1:3294424/148‑1 (MQ=255)
                                                                                       gATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGcc                                                            <  1:2044891/148‑1 (MQ=255)
                                                                                       gATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGcc                                                            <  1:334549/148‑1 (MQ=255)
                                                                                       gATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGGTACGTTACTGGGGAGcc                                                            >  2:334218/1‑148 (MQ=255)
                                                                                          atTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTgc                                                         >  2:543688/1‑148 (MQ=255)
                                                                                            tCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTgcgc                                                       >  1:234397/1‑148 (MQ=255)
                                                                                            tCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTgcgc                                                       >  2:1415775/1‑148 (MQ=255)
                                                                                              acaGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTACAGGTGGTTACGTTACTGGGGAGCCTGCGCaa                                                     <  2:101300/148‑1 (MQ=255)
                                                                                              acaGCCGCACCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCaa                                                     <  1:2077517/148‑1 (MQ=255)
                                                                                                  ccGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTCCTGGGGAGCCTGCGCAAAGGc                                                 >  2:2152670/1‑148 (MQ=255)
                                                                                                   cgcgCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCt                                                <  1:420729/148‑1 (MQ=255)
                                                                                                      gcCAACAATCCTACCAGCCAGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCAt                                             <  1:1640907/148‑1 (MQ=255)
                                                                                                       ccAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCAtt                                            >  2:2342045/1‑148 (MQ=255)
                                                                                                        caacaaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCAttt                                           >  2:1203028/1‑148 (MQ=255)
                                                                                                           caaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAAt                                        <  1:1405526/148‑1 (MQ=255)
                                                                                                           caaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAAt                                        >  1:2036668/1‑148 (MQ=255)
                                                                                                            aaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATg                                       <  1:46008/148‑1 (MQ=255)
                                                                                                            aaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATg                                       >  1:1495015/1‑148 (MQ=255)
                                                                                                              tCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGc                                     <  2:950008/148‑1 (MQ=255)
                                                                                                                cTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCAt                                   >  2:824554/1‑148 (MQ=255)
                                                                                                                  aCCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATgg                                 >  1:625407/1‑148 (MQ=255)
                                                                                                                  aCCAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATgg                                 >  2:2503411/1‑148 (MQ=255)
                                                                                                                    cAGCCGGATGATGCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGtt                               >  2:2534466/1‑148 (MQ=255)
                                                                                                                          gatgatGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGt                         <  1:3064024/148‑1 (MQ=255)
                                                                                                                             gatgCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTAcc                      >  2:2888977/1‑148 (MQ=255)
                                                                                                                             gatgCGTTTTGCCCAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTAcc                      <  1:1912944/148‑1 (MQ=255)
                                                                                                                              atgCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTAccc                     >  2:1505687/1‑148 (MQ=255)
                                                                                                                               tgCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCt                    >  1:3065791/1‑148 (MQ=255)
                                                                                                                                    tttGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGa               >  1:2841566/1‑148 (MQ=255)
                                                                                                                                    tttGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGa               >  2:474700/1‑148 (MQ=255)
                                                                                                                                    attGCCCAACTCTCGACCGTGAAGGCTTGTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGa               <  2:2039208/147‑1 (MQ=255)
                                                                                                                                     ttGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGaa              >  1:3095650/1‑148 (MQ=255)
                                                                                                                                      tGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGaaa             <  1:2687388/148‑1 (MQ=255)
                                                                                                                                        cccAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAAt           <  2:1602706/148‑1 (MQ=255)
                                                                                                                                        cccAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAAt           <  2:1602908/148‑1 (MQ=255)
                                                                                                                                        cccAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAAt           <  2:1603209/148‑1 (MQ=255)
                                                                                                                                          cAACTCTCAACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTg         <  2:236950/148‑1 (MQ=255)
                                                                                                                                            aCTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCt       >  1:3216042/1‑148 (MQ=255)
                                                                                                                                              tctctACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTcc     >  2:769482/1‑148 (MQ=255)
                                                                                                                                              tctctACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTcc     >  1:100899/1‑148 (MQ=255)
                                                                                                                                                tctACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCgg   <  2:2064900/148‑1 (MQ=255)
                                                                                                                                                 ctACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGc  >  2:2272869/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGGTGTGGTCCGTGGCGGTCACGCCCACGATTGAAAAGCTCAGTGTCTGGGAGTTAGATGGCAAACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGC  >  CP009273/3887817‑3888109

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: