Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,692,994 T→A 91.1% 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,9940TA91.1% 479.0 / 31.3 158F474Y (TTC→TAC) bcsGDUF3260 family inner membrane protein associated with cellulose production
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base A (62/82);  total (71/87)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 165 total aligned reads displayed.

TGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCA  >  CP009273/3692851‑3693137
                                                                                                                                               |                                                                                                                                               
ttaaCTTGAGAAATCGGGTCGTATAGTTATGGTGGTCGTGGTGCCGGAACACGGCGGCGCTCTGAAGGGCGACAGAATGCATGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCAGTTGGGGTGAAAtacttc                                                                                                                                             <  1:2151297/146‑1 (MQ=255)
tGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtacttc                                                                                                                                             >  1:2172359/1‑148 (MQ=255)
tGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtacttc                                                                                                                                             <  1:211098/148‑1 (MQ=255)
tGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtacttc                                                                                                                                             >  1:2089433/1‑148 (MQ=255)
tGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtacttc                                                                                                                                             >  1:2172664/1‑148 (MQ=255)
  atattGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTTTCACCGACGTCCCCGTTGGGGTGAAATACTTCgg                                                                                                                                           <  2:1533971/145‑1 (MQ=255)
   aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGc                                                                                                                                          <  1:395559/148‑1 (MQ=255)
   aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGc                                                                                                                                          <  1:380674/148‑1 (MQ=255)
     ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCAt                                                                                                                                        >  1:2455103/1‑148 (MQ=255)
        agaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGaa                                                                                                                                     <  1:2417572/148‑1 (MQ=255)
          aaaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAgg                                                                                                                                   >  1:115189/1‑148 (MQ=255)
          aaaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAgg                                                                                                                                   >  1:2002323/1‑148 (MQ=255)
           aaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGc                                                                                                                                  <  2:2063926/148‑1 (MQ=255)
            aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                 <  1:1078715/148‑1 (MQ=255)
             tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                >  2:2636578/1‑148 (MQ=255)
                ggTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGc                                                                                                                             <  1:78898/148‑1 (MQ=255)
                     tAAAGTGATGGTGGTCGTGGTGACGGAACACGGCGGCGCGCTGAAGTGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGAACCGCATCAg                                                                                                                        >  2:2284717/1‑148 (MQ=255)
                       aaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGGTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggg                                                                                                                      >  2:2129899/1‑148 (MQ=255)
                          tGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGc                                                                                                                   >  1:1970871/1‑148 (MQ=255)
                             tggtggTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCAcc                                                                                                                >  1:2157778/1‑148 (MQ=255)
                                tggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGAt                                                                                                             <  2:2826797/148‑1 (MQ=255)
                                 ggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGAtt                                                                                                            >  2:551016/1‑148 (MQ=255)
                                 ggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGAtt                                                                                                            >  2:57947/1‑148 (MQ=255)
                                  gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTg                                                                                                           <  1:944265/148‑1 (MQ=255)
                                  gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTg                                                                                                           >  1:1543162/1‑148 (MQ=255)
                                   tCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGt                                                                                                          >  1:713499/1‑148 (MQ=255)
                                            cGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaaca                                                                                                 <  2:1798103/148‑1 (MQ=255)
                                              gAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaacaac                                                                                               >  1:3305536/1‑148 (MQ=255)
                                                 cacGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGa                                                                                            >  2:2778420/1‑148 (MQ=255)
                                                  acGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGag                                                                                           <  1:870585/148‑1 (MQ=255)
                                                  acGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGag                                                                                           <  1:3045048/148‑1 (MQ=255)
                                                     gcggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcag                                                                                        <  2:2809494/148‑1 (MQ=255)
                                                      cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                       <  1:3248591/148‑1 (MQ=255)
                                                           cgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCt                                                                                  <  2:649015/148‑1 (MQ=255)
                                                             cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                <  1:600387/148‑1 (MQ=255)
                                                             cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                <  2:797060/148‑1 (MQ=255)
                                                              tGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGc                                                                               <  1:329870/148‑1 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           >  2:2780971/1‑148 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           <  1:49767/148‑1 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           >  2:1265169/1‑148 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           >  1:1903786/1‑148 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           >  1:1631642/1‑148 (MQ=255)
                                                                  gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                           >  1:1408302/1‑148 (MQ=255)
                                                                   ggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtct                                                                          >  1:1443868/1‑148 (MQ=255)
                                                                    gCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtctc                                                                         <  1:797205/148‑1 (MQ=255)
                                                                    gCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtctc                                                                         <  2:377961/148‑1 (MQ=255)
                                                                        cAGAATGCAGGTATCTGGCCTACGTGATATCACTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGAt                                                                     <  1:2390038/148‑1 (MQ=255)
                                                                         aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATc                                                                    <  1:2616953/148‑1 (MQ=255)
                                                                         aGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATc                                                                    <  1:2616962/148‑1 (MQ=255)
                                                                               cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                              <  1:1725810/148‑1 (MQ=255)
                                                                               cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACAGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                              <  1:1095815/148‑1 (MQ=255)
                                                                               cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                              <  2:1146599/148‑1 (MQ=255)
                                                                                aGGTATCTGGCCTACGTGATATCCCTAGCCCCTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTggtggt                                                             <  2:1205692/148‑1 (MQ=255)
                                                                                  gTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTCTCTCCGATCTGGTGGTTc                                                           >  2:1717794/1‑148 (MQ=255)
                                                                                      cTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGt                                                       >  1:441637/1‑148 (MQ=255)
                                                                                          ccTACGTGATATCCCTAGCCCGTCTATCACTGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtctc                                                   <  2:1021134/148‑1 (MQ=255)
                                                                                          ccTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtctc                                                   >  1:2793276/1‑148 (MQ=255)
                                                                                           cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCg                                                  >  1:1636429/1‑148 (MQ=255)
                                                                                           cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCg                                                  <  2:577620/148‑1 (MQ=255)
                                                                                            tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                 <  1:804165/148‑1 (MQ=255)
                                                                                             aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGAt                                                >  1:3083349/1‑148 (MQ=255)
                                                                                             aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGAt                                                >  1:3085947/1‑148 (MQ=255)
                                                                                                tGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGc                                             >  1:404855/1‑148 (MQ=255)
                                                                                                 gATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCa                                            >  1:673311/1‑148 (MQ=255)
                                                                                                   tatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAg                                          >  2:337162/1‑148 (MQ=255)
                                                                                                   tatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAg                                          <  1:1114822/148‑1 (MQ=255)
                                                                                                   tatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAg                                          <  1:3139091/148‑1 (MQ=255)
                                                                                                    atCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGa                                         >  1:1765570/1‑148 (MQ=255)
                                                                                                    atCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGa                                         >  1:2892882/1‑148 (MQ=255)
                                                                                                     tCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAt                                        <  1:176687/148‑1 (MQ=255)
                                                                                                     tCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAt                                        <  2:225294/148‑1 (MQ=255)
                                                                                                      cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtt                                       >  2:170747/1‑148 (MQ=255)
                                                                                                      cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtt                                       <  1:2373044/148‑1 (MQ=255)
                                                                                                      cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtt                                       <  1:20292/148‑1 (MQ=255)
                                                                                                      cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtt                                       <  2:2609261/148‑1 (MQ=255)
                                                                                                        cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtttt                                     >  1:482354/1‑148 (MQ=255)
                                                                                                        cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtttt                                     <  1:1123378/148‑1 (MQ=255)
                                                                                                         tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTc                                    >  2:406316/1‑148 (MQ=255)
                                                                                                          aGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCa                                   <  1:1265169/148‑1 (MQ=255)
                                                                                                              cGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGa                               <  1:3074536/148‑1 (MQ=255)
                                                                                                                tCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAg                             >  2:2140436/1‑148 (MQ=255)
                                                                                                                tCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAg                             >  1:215005/1‑148 (MQ=255)
                                                                                                                 cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                            <  2:2576214/148‑1 (MQ=255)
                                                                                                                   aTCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACa                          <  2:2424015/148‑1 (MQ=255)
                                                                                                                    tCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACaa                         <  2:1970871/148‑1 (MQ=255)
                                                                                                                      aCCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATg                       <  2:2951869/148‑1 (MQ=255)
                                                                                                                       ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                      <  2:2455103/148‑1 (MQ=255)
                                                                                                                       ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                      <  1:3103110/148‑1 (MQ=255)
                                                                                                                          aCGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGa                   <  2:838297/148‑1 (MQ=255)
                                                                                                                            gTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACt                 <  2:2302874/148‑1 (MQ=255)
                                                                                                                              ccccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTgg               >  2:3075522/1‑148 (MQ=255)
                                                                                                                              ccccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTgg               >  1:133218/1‑148 (MQ=255)
                                                                                                                                      gggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACt       <  1:551016/148‑1 (MQ=255)
                                                                                                                                       ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc      <  2:3085947/148‑1 (MQ=255)
                                                                                                                                       ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc      <  2:2002323/148‑1 (MQ=255)
                                                                                                                                       ggTCAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGACAACGAAAAACCTATCAGCTACCTGGCTATCTCCGATCTGGTGGTTCTCTTTCTCGATTGCAAGATTTTCACCGAAGACAAAGTTGACTGGAAAAAACTc      >  2:697688/1‑148 (MQ=255)
                                                                                                                                        gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa     <  2:3057905/148‑1 (MQ=255)
                                                                                                                                         tGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCAc    <  1:1557354/148‑1 (MQ=255)
                                                                                                                                          gAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCAcc   >  2:2324591/1‑148 (MQ=255)
                                                                                                                                           aaaTACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCa  <  1:851323/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCA  >  CP009273/3692851‑3693137

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: