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.9% 959.7 / 6.0 290F474Y (TTC→TAC) bcsGDUF3260 family inner membrane protein associated with cellulose production
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base A (144/140);  total (148/142)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 297 total aligned reads displayed.

TACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAGTG  >  CP009273/3692848‑3693140
                                                                                                                                                  |                                                                                                                                                  
tACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtactt                                                                                                                                                 >  1:1527295/1‑150 (MQ=255)
       cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                          >  2:2342130/1‑150 (MQ=255)
       cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                          <  1:1114551/150‑1 (MQ=255)
       cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                          >  1:823765/1‑150 (MQ=255)
       cTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                          <  2:1497858/150‑1 (MQ=255)
        ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCCGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGa                                                                                                                                         <  1:2039273/150‑1 (MQ=255)
        ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGa                                                                                                                                         >  2:4314956/1‑150 (MQ=255)
         tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGaa                                                                                                                                        >  1:4278538/1‑150 (MQ=255)
         tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGaa                                                                                                                                        >  2:2179199/1‑150 (MQ=255)
          gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAg                                                                                                                                       <  1:41094/150‑1 (MQ=255)
          gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAg                                                                                                                                       <  1:643434/150‑1 (MQ=255)
           agaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTGGCCGTTGGGGTGAAATACTTCGGCATGAAgg                                                                                                                                      <  2:3338984/150‑1 (MQ=255)
            gaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGc                                                                                                                                     >  2:4871656/1‑150 (MQ=255)
            gaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGc                                                                                                                                     >  1:4772106/1‑150 (MQ=255)
             aaaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                    <  1:3251760/150‑1 (MQ=255)
              aaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                   <  2:4049150/150‑1 (MQ=255)
               aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAcc                                                                                                                                  <  2:903391/150‑1 (MQ=255)
                tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCg                                                                                                                                 >  1:3528428/1‑150 (MQ=255)
                tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCg                                                                                                                                 <  1:334956/150‑1 (MQ=255)
                   ggTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCAt                                                                                                                              <  1:539137/150‑1 (MQ=255)
                     tCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCa                                                                                                                            <  2:1018840/150‑1 (MQ=255)
                      cGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAg                                                                                                                           <  1:4445565/150‑1 (MQ=255)
                        tAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggg                                                                                                                         <  1:2645495/150‑1 (MQ=255)
                              gATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCAca                                                                                                                   <  2:1668082/150‑2 (MQ=255)
                                tggtggTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGa                                                                                                                 <  1:3667455/150‑1 (MQ=255)
                                  gtggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGAtt                                                                                                               <  1:4965092/150‑1 (MQ=255)
                                    ggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGt                                                                                                             >  1:5074452/1‑150 (MQ=255)
                                     gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTc                                                                                                            <  1:2758758/150‑1 (MQ=255)
                                      tCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCa                                                                                                           >  1:3044874/1‑150 (MQ=255)
                                        gtggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATc                                                                                                         >  1:1000440/1‑150 (MQ=255)
                                         tggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCg                                                                                                        >  1:2357356/1‑150 (MQ=255)
                                          ggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGa                                                                                                       <  2:1664978/150‑1 (MQ=255)
                                          ggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGa                                                                                                       >  1:1468778/1‑150 (MQ=255)
                                               cGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaacaac                                                                                                  >  1:582611/1‑150 (MQ=255)
                                                  aaCACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGa                                                                                               >  1:522714/1‑150 (MQ=255)
                                                       ggcggcGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                          >  2:1889035/1‑150 (MQ=255)
                                                       ggcggcGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                          >  2:3856453/1‑150 (MQ=255)
                                                         cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCtt                                                                                        <  1:2258809/150‑1 (MQ=255)
                                                         cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCtt                                                                                        <  1:4822725/150‑1 (MQ=255)
                                                         cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCtt                                                                                        >  2:4288291/1‑150 (MQ=255)
                                                              cgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                   >  2:3415709/1‑150 (MQ=255)
                                                               gcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGc                                                                                  >  1:756462/1‑150 (MQ=255)
                                                                cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCt                                                                                 <  2:1337887/150‑1 (MQ=255)
                                                                cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCt                                                                                 <  2:2399649/150‑1 (MQ=255)
                                                                  gAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAt                                                                               >  1:4496234/1‑150 (MQ=255)
                                                                  gAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAt                                                                               >  1:3124347/1‑150 (MQ=255)
                                                                   aaGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                              >  2:743481/1‑150 (MQ=255)
                                                                    aGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtct                                                                             >  1:2707553/1‑150 (MQ=255)
                                                                      ggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTcc                                                                           <  2:3183987/150‑1 (MQ=255)
                                                                       gCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCg                                                                          <  2:140383/150‑1 (MQ=255)
                                                                        cGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGa                                                                         <  2:4103906/150‑1 (MQ=255)
                                                                         gACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGAt                                                                        >  2:3818814/1‑150 (MQ=255)
                                                                            aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtg                                                                     >  1:489911/1‑150 (MQ=255)
                                                                             gAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTACCTCCGATCtgg                                                                    >  2:4579795/1‑150 (MQ=255)
                                                                              aaTGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggt                                                                   >  2:4686458/1‑150 (MQ=255)
                                                                               aTGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtg                                                                  >  1:1192422/1‑150 (MQ=255)
                                                                                tGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                                 >  2:1170118/1‑150 (MQ=255)
                                                                                  cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGtt                                                               <  2:2874682/150‑1 (MQ=255)
                                                                                   aGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTc                                                              <  1:4314703/150‑1 (MQ=255)
                                                                                   aGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTc                                                              <  2:3075385/150‑1 (MQ=255)
                                                                                      tATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                           >  2:5056406/1‑149 (MQ=255)
                                                                                      tATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                           <  1:677220/150‑2 (MQ=255)
                                                                                       aTCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgctt                                                          <  1:2577652/150‑3 (MQ=255)
                                                                                          tGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGTAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtct                                                       >  2:2713727/1‑150 (MQ=255)
                                                                                               tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATg                                                  <  2:68289/150‑1 (MQ=255)
                                                                                                 cGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGc                                                <  2:353047/150‑1 (MQ=255)
                                                                                                  gTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCa                                               >  2:3905813/1‑150 (MQ=255)
                                                                                                   tGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCaa                                              <  1:1289324/150‑1 (MQ=255)
                                                                                                    gATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAg                                             <  1:4594421/150‑1 (MQ=255)
                                                                                                       atCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtt                                          <  1:437642/150‑1 (MQ=255)
                                                                                                         cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtttt                                        >  2:2635177/1‑150 (MQ=255)
                                                                                                         cccTAGCCCGTCTATCACCGACGTCCCCGGTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtttt                                        >  1:1942815/1‑150 (MQ=255)
                                                                                                           cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCa                                      <  2:2749868/150‑1 (MQ=255)
                                                                                                           cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCa                                      <  1:2179199/150‑1 (MQ=255)
                                                                                                           cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCa                                      <  2:5106228/150‑1 (MQ=255)
                                                                                                            tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCTGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAc                                     >  1:2810750/1‑150 (MQ=255)
                                                                                                             aGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAcc                                    >  1:5147688/1‑150 (MQ=255)
                                                                                                              gCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCg                                   <  1:3756304/150‑1 (MQ=255)
                                                                                                                ccttCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGaa                                 >  1:1060739/4‑150 (MQ=255)
                                                                                                                 cGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAg                                >  1:3397277/1‑150 (MQ=255)
                                                                                                                  gTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                               >  1:2976259/1‑150 (MQ=255)
                                                                                                                    cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACa                             >  2:3871646/1‑150 (MQ=255)
                                                                                                                     tATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACaa                            >  1:3198398/1‑150 (MQ=255)
                                                                                                                       tCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATg                          <  2:4215321/150‑1 (MQ=255)
                                                                                                                        cACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                         >  2:1788669/1‑150 (MQ=255)
                                                                                                                         aCCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACACTGtt                        >  2:1949535/1‑150 (MQ=255)
                                                                                                                             aCGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACt                    >  2:4895146/1‑150 (MQ=255)
                                                                                                                             aCGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACt                    >  1:2443322/1‑150 (MQ=255)
                                                                                                                               gTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTgg                  >  1:2439955/1‑150 (MQ=255)
                                                                                                                                  cccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaa               <  2:232795/150‑1 (MQ=255)
                                                                                                                                    cGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa             >  1:1877953/1‑150 (MQ=255)
                                                                                                                                    cGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa             >  2:4983681/1‑150 (MQ=255)
                                                                                                                                       tGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACt          <  2:22409/150‑1 (MQ=255)
                                                                                                                                        ggggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc         >  1:1680999/1‑150 (MQ=255)
                                                                                                                                        ggggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc         <  2:582611/150‑1 (MQ=255)
                                                                                                                                          ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCAc       >  2:4381286/1‑150 (MQ=255)
                                                                                                                                          ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCAc       >  2:3282605/1‑150 (MQ=255)
                                                                                                                                            tGAAATACTTCGGCATGAAGGCACCGCATCTGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGTCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGTAAAAACTCACCa     <  1:4298314/150‑1 (MQ=255)
                                                                                                                                             gAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAg    <  2:4772106/150‑1 (MQ=255)
                                                                                                                                              aaaTACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAGt   >  1:3570496/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAGTG  >  CP009273/3692848‑3693140

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: ‑