Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,692,994 T→A 100% F474Y (TTC→TAC)  bcsG → DUF3260 family inner membrane protein associated with cellulose production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,692,9940TA97.0% 573.4 / 4.8 180F474Y (TTC→TAC) bcsGDUF3260 family inner membrane protein associated with cellulose production
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base A (87/87);  total (87/93)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 188 total aligned reads displayed.

ACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAG  >  CP009273/3692849‑3693138
                                                                                                                                                 |                                                                                                                                                
aCTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtact                                                                                                                                                <  1:2949445/148‑1 (MQ=255)
aCTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtact                                                                                                                                                <  1:2161568/148‑1 (MQ=255)
aCTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtact                                                                                                                                                <  1:2948651/148‑1 (MQ=255)
 cTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtactt                                                                                                                                               <  2:2308017/148‑1 (MQ=255)
 cTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtactt                                                                                                                                               >  1:790213/1‑148 (MQ=255)
   gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                             >  1:3427238/1‑148 (MQ=255)
   gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                             >  1:3426806/1‑148 (MQ=255)
   gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                             >  2:3074861/1‑148 (MQ=255)
   gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                             >  2:381868/1‑148 (MQ=255)
     aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGc                                                                                                                                           <  2:2837355/148‑1 (MQ=255)
     aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGc                                                                                                                                           <  2:2837450/148‑1 (MQ=255)
     aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGc                                                                                                                                           <  2:2837418/148‑1 (MQ=255)
      cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCa                                                                                                                                          >  2:3131872/1‑148 (MQ=255)
      cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCa                                                                                                                                          >  2:2129172/1‑148 (MQ=255)
       ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCa                                                                                                                                          >  1:3242343/1‑147 (MQ=255)
       ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCAt                                                                                                                                         >  2:1996847/1‑148 (MQ=255)
        tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                        >  2:3189615/1‑148 (MQ=255)
        tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                        >  2:2743130/1‑148 (MQ=255)
         gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGa                                                                                                                                       <  1:106162/148‑1 (MQ=255)
         gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCAGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGa                                                                                                                                       >  2:1806205/1‑148 (MQ=255)
            aaaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAgg                                                                                                                                    >  2:2748426/1‑148 (MQ=255)
              aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                  >  1:357181/1‑148 (MQ=255)
              aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                  >  1:476215/1‑148 (MQ=255)
                   gTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCa                                                                                                                             >  2:2615547/1‑148 (MQ=255)
                        aaaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgg                                                                                                                        <  1:2908193/148‑1 (MQ=255)
                         aaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggg                                                                                                                       >  1:1859626/1‑148 (MQ=255)
                          aGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgggg                                                                                                                      >  1:1013146/1‑148 (MQ=255)
                           gTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggggg                                                                                                                     <  1:892372/148‑1 (MQ=255)
                             gATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCa                                                                                                                   >  1:1016561/1‑148 (MQ=255)
                             gATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCa                                                                                                                   >  1:1017258/1‑148 (MQ=255)
                                 gtggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGa                                                                                                               >  1:1190011/1‑148 (MQ=255)
                                     tCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGt                                                                                                           <  2:2611168/148‑1 (MQ=255)
                                        tggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCAt                                                                                                        <  2:2891971/148‑1 (MQ=255)
                                           tgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGa                                                                                                     >  1:281251/1‑148 (MQ=255)
                                            gCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaa                                                                                                    <  1:1459952/148‑1 (MQ=255)
                                            gCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaa                                                                                                    <  2:2392795/148‑1 (MQ=255)
                                                gAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaacaac                                                                                                >  1:2551225/1‑148 (MQ=255)
                                                gAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaacaac                                                                                                >  2:2670525/1‑148 (MQ=255)
                                                  acacGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCg                                                                                              <  2:1933359/148‑1 (MQ=255)
                                                   cacGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGa                                                                                             <  2:454986/148‑1 (MQ=255)
                                                       gcggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCAGCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcag                                                                                         >  2:1384828/1‑148 (MQ=255)
                                                       gcggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcag                                                                                         <  2:1762575/148‑1 (MQ=255)
                                                       gcggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcag                                                                                         <  1:1560070/148‑1 (MQ=255)
                                                           cgcgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTc                                                                                     <  1:34293/148‑1 (MQ=255)
                                                           cgcgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTc                                                                                     <  1:2891396/148‑1 (MQ=255)
                                                               cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCCCCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                 <  1:3276562/148‑1 (MQ=255)
                                                                 gAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCt                                                                               <  1:2743282/148‑1 (MQ=255)
                                                                   aGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAt                                                                             <  1:2378954/148‑1 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            <  2:2211893/148‑1 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            <  2:1160876/148‑1 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  1:521822/1‑148 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  2:2811636/1‑148 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  2:1854800/1‑148 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATAATTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  2:321569/1‑148 (MQ=255)
                                                                       cGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTcc                                                                         <  2:156030/148‑1 (MQ=255)
                                                                        gACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCg                                                                        <  1:3171219/148‑1 (MQ=255)
                                                                         aCAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGa                                                                       <  1:118394/148‑1 (MQ=255)
                                                                           aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATc                                                                     >  2:2184438/1‑148 (MQ=255)
                                                                           aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATc                                                                     <  2:1179707/148‑1 (MQ=255)
                                                                           aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATc                                                                     <  2:1179739/148‑1 (MQ=255)
                                                                            gAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCt                                                                    >  2:1843684/1‑148 (MQ=255)
                                                                              aTGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtgg                                                                  <  2:3277471/148‑1 (MQ=255)
                                                                               tGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggt                                                                 <  1:3500014/148‑1 (MQ=255)
                                                                                gCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtg                                                                <  2:1921967/148‑1 (MQ=255)
                                                                                 cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                               <  2:1264223/148‑1 (MQ=255)
                                                                                  aGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTggtggt                                                              <  1:3278851/148‑1 (MQ=255)
                                                                                   ggTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGtt                                                             <  2:2639134/148‑1 (MQ=255)
                                                                                    gTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTc                                                            >  2:1012314/1‑148 (MQ=255)
                                                                                      aTCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgc                                                          >  1:1101548/1‑148 (MQ=255)
                                                                                       tCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                         <  2:1611546/148‑2 (MQ=255)
                                                                                       tCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                         <  2:1611569/148‑2 (MQ=255)
                                                                                        cTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgctt                                                        <  2:661027/148‑3 (MQ=255)
                                                                                        cTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgctt                                                        >  2:1048958/1‑146 (MQ=255)
                                                                                          ggCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtc                                                      >  2:1495968/1‑148 (MQ=255)
                                                                                           gCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtct                                                     >  1:1235318/1‑148 (MQ=255)
                                                                                             cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCg                                                   <  1:2843032/148‑1 (MQ=255)
                                                                                             cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCg                                                   <  1:633242/148‑1 (MQ=255)
                                                                                              tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                  <  1:1802975/148‑1 (MQ=255)
                                                                                                           tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTc                                     >  2:921220/1‑148 (MQ=255)
                                                                                                           tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTc                                     >  2:1247780/1‑148 (MQ=255)
                                                                                                              cccGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAcc                                  >  1:1350901/1‑148 (MQ=255)
                                                                                                                   cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                             <  1:2717515/148‑1 (MQ=255)
                                                                                                                   cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                             <  1:2717916/148‑1 (MQ=255)
                                                                                                                     aTCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACa                           <  1:554837/148‑1 (MQ=255)
                                                                                                                        aCCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATg                        <  1:2670525/148‑1 (MQ=255)
                                                                                                                         ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  2:2098149/148‑1 (MQ=255)
                                                                                                                         ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  1:1320392/148‑1 (MQ=255)
                                                                                                                         ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  2:326785/148‑1 (MQ=255)
                                                                                                                             cGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGAc                   >  2:3239288/1‑148 (MQ=255)
                                                                                                                                 cccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGa               >  1:132044/1‑148 (MQ=255)
                                                                                                                                  ccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGACTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaa              >  2:146973/1‑148 (MQ=255)
                                                                                                                                    gTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaa            >  1:742050/1‑148 (MQ=255)
                                                                                                                                     ttGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa           <  1:2090686/148‑1 (MQ=255)
                                                                                                                                      tGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaaa          >  1:283654/1‑148 (MQ=255)
                                                                                                                                      tGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaaa          >  2:1855573/1‑148 (MQ=255)
                                                                                                                                         ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc       >  2:1755434/1‑148 (MQ=255)
                                                                                                                                          gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa      >  2:2030766/1‑148 (MQ=255)
                                                                                                                                          gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa      >  1:2243122/1‑148 (MQ=255)
                                                                                                                                          gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa      >  2:258583/1‑148 (MQ=255)
                                                                                                                                             aaaTACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCa   >  2:756009/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAG  >  CP009273/3692849‑3693138

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: