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,9940TA96.1% 530.2 / 6.7 165F474Y (TTC→TAC) bcsGDUF3260 family inner membrane protein associated with cellulose production
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base A (73/85);  total (75/90)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 172 total aligned reads displayed.

TACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAG  >  CP009273/3692848‑3693138
                                                                                                                                                  |                                                                                                                                                
tACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAttc                                                                                                                                                 <  1:1965450/148‑1 (MQ=255)
 aCTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGANAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtact                                                                                                                                                <  1:1803647/148‑1 (MQ=255)
 aCTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtact                                                                                                                                                <  1:1144828/148‑1 (MQ=255)
    gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                             >  2:84396/1‑148 (MQ=255)
     aaCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCgg                                                                                                                                            <  1:3257934/148‑1 (MQ=255)
     aaCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCgg                                                                                                                                            <  1:3257935/148‑1 (MQ=255)
     aaCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCgg                                                                                                                                            >  1:2508599/1‑148 (MQ=255)
      aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGc                                                                                                                                           <  2:1030332/148‑1 (MQ=255)
      aCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGc                                                                                                                                           <  2:1103236/148‑1 (MQ=255)
         tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                        <  2:2937215/148‑1 (MQ=255)
           agaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGaa                                                                                                                                      <  2:1600523/148‑1 (MQ=255)
               aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                  <  1:1642126/148‑1 (MQ=255)
               aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                  >  1:576744/1‑148 (MQ=255)
                tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                 <  1:1191404/148‑1 (MQ=255)
                tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                 >  1:685885/1‑148 (MQ=255)
                tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                 <  2:2154276/148‑1 (MQ=255)
                           aGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgggg                                                                                                                      >  1:250577/1‑148 (MQ=255)
                            gTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggggg                                                                                                                     >  2:2089649/1‑148 (MQ=255)
                            gTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGAAGGTATCTGGCCTACGTGATATCCCTAGACCGTCTATCACCGACGTCCCCGTTGGTGTGAAATACATCGGCATGAAGGCACCTCATCAggggg                                                                                                                     >  2:2089726/1‑148 (MQ=255)
                                        gtggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCACTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCa                                                                                                         <  1:2580970/148‑1 (MQ=255)
                                         tggtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCAt                                                                                                        >  2:2609410/1‑148 (MQ=255)
                                            tgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGa                                                                                                     <  2:1979238/148‑1 (MQ=255)
                                             gCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaa                                                                                                    <  1:1482902/148‑1 (MQ=255)
                                             gCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaa                                                                                                    <  1:921804/148‑1 (MQ=255)
                                               cGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaaca                                                                                                  <  2:2322162/148‑1 (MQ=255)
                                               cGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaaca                                                                                                  <  2:2322081/148‑1 (MQ=255)
                                                  aaCACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAAcc                                                                                               >  1:1666158/1‑148 (MQ=255)
                                                   acacGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCg                                                                                              <  2:2525114/148‑1 (MQ=255)
                                                      cggcggCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagc                                                                                           <  1:3450674/148‑1 (MQ=255)
                                                         cgtcgcgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                        >  1:2884642/4‑148 (MQ=255)
                                                         cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                        <  1:1885262/148‑1 (MQ=255)
                                                            cgcgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTTATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTc                                                                                     <  2:1390644/148‑1 (MQ=255)
                                                              cgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCt                                                                                   >  2:1034787/1‑148 (MQ=255)
                                                                cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                 <  1:663502/148‑1 (MQ=255)
                                                                cTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGACCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTgg                                                                                 >  1:1988832/1‑148 (MQ=255)
                                                                 tGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGc                                                                                >  2:116315/1‑148 (MQ=255)
                                                                   aaGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTa                                                                              >  2:3146370/1‑148 (MQ=255)
                                                                   aaGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTa                                                                              >  2:2520285/1‑148 (MQ=255)
                                                                     gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  1:898409/1‑148 (MQ=255)
                                                                     gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                            >  1:3140107/1‑148 (MQ=255)
                                                                      ggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtct                                                                           >  2:1867042/1‑148 (MQ=255)
                                                                        cGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTcc                                                                         <  2:889087/148‑1 (MQ=255)
                                                                          aCAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGa                                                                       >  1:2008624/1‑148 (MQ=255)
                                                                              aaTGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTCTCTCCGATCtg                                                                   >  2:3084914/1‑148 (MQ=255)
                                                                              aaTGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtg                                                                   <  1:84396/148‑1 (MQ=255)
                                                                                  cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                               >  1:2827519/1‑148 (MQ=255)
                                                                                    ggTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGtt                                                             >  1:1010673/1‑148 (MQ=255)
                                                                                       aacTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACACCCGAGCAGCTTCCTGGCTCTCTCCGATCTGGTGGTTcgc                                                          >  1:689660/3‑148 (MQ=255)
                                                                                        tCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                         >  1:1795439/1‑147 (MQ=255)
                                                                                         cTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgctt                                                        >  2:330372/1‑146 (MQ=255)
                                                                                         cTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGt                                                        >  1:3020360/1‑148 (MQ=255)
                                                                                             ccTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTtctc                                                    >  1:3166938/1‑148 (MQ=255)
                                                                                             ccTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTtctc                                                    >  1:3166896/1‑148 (MQ=255)
                                                                                              cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCg                                                   <  2:3254637/148‑1 (MQ=255)
                                                                                              cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCg                                                   <  1:2895973/148‑1 (MQ=255)
                                                                                              cTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCg                                                   <  1:2455688/148‑1 (MQ=255)
                                                                                               tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                  <  1:956756/148‑1 (MQ=255)
                                                                                               tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                  <  1:3037304/148‑1 (MQ=255)
                                                                                               tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGa                                                  <  2:528157/148‑1 (MQ=255)
                                                                                                aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGAt                                                 <  1:2133986/148‑1 (MQ=255)
                                                                                                aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGAt                                                 <  1:2385094/148‑1 (MQ=255)
                                                                                                 cGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATg                                                >  1:1291704/1‑148 (MQ=255)
                                                                                                    gATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCa                                             >  1:689576/1‑148 (MQ=255)
                                                                                                    gATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCa                                             <  2:1737494/148‑1 (MQ=255)
                                                                                                     atatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCaa                                            >  2:396937/1‑148 (MQ=255)
                                                                                                       atCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGa                                          >  1:2968469/1‑148 (MQ=255)
                                                                                                          ccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAttt                                       >  1:1593301/1‑148 (MQ=255)
                                                                                                            tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTc                                     <  2:3517250/148‑1 (MQ=255)
                                                                                                             aGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCa                                    <  2:626374/148‑1 (MQ=255)
                                                                                                                 cGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGa                                <  2:576744/148‑1 (MQ=255)
                                                                                                                  gTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGaa                               >  2:956392/1‑148 (MQ=255)
                                                                                                                   tCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAg                              <  1:883801/148‑1 (MQ=255)
                                                                                                                    cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                             <  1:1054500/148‑1 (MQ=255)
                                                                                                                    cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGa                             <  2:2918980/148‑1 (MQ=255)
                                                                                                                    cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGa                             >  2:483073/1‑148 (MQ=255)
                                                                                                                    cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGa                             >  1:699447/1‑148 (MQ=255)
                                                                                                                        cACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAAt                         >  2:1640984/1‑148 (MQ=255)
                                                                                                                         aCCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATg                        <  2:687189/148‑1 (MQ=255)
                                                                                                                         aCCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATg                        <  2:685871/148‑1 (MQ=255)
                                                                                                                          ccGACGTCCGCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  1:2274630/148‑1 (MQ=255)
                                                                                                                          ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACTGAAGACAATGt                       <  2:339951/148‑1 (MQ=255)
                                                                                                                          ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  1:2274594/148‑1 (MQ=255)
                                                                                                                          ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  2:1652231/148‑1 (MQ=255)
                                                                                                                          ccGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGt                       <  2:3144900/148‑1 (MQ=255)
                                                                                                                           cGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGtt                      >  2:180821/1‑148 (MQ=255)
                                                                                                                           cGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGtt                      >  2:2382498/1‑148 (MQ=255)
                                                                                                                               gTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACt                  <  2:1101906/148‑1 (MQ=255)
                                                                                                                                 ccccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTgg                >  1:2683528/1‑148 (MQ=255)
                                                                                                                                  cccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGa               >  1:3351888/1‑148 (MQ=255)
                                                                                                                                  cccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGa               <  2:2508599/148‑1 (MQ=255)
                                                                                                                                   ccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcg                                                           >  1:1104273/1‑103 (MQ=255)
                                                                                                                                     gTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaa            <  2:1195659/148‑1 (MQ=255)
                                                                                                                                     gTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaa            <  2:2176156/148‑1 (MQ=255)
                                                                                                                                      ttGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa           >  1:4617/1‑148 (MQ=255)
                                                                                                                                      ttGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa           >  1:4394/1‑148 (MQ=255)
                                                                                                                                         gggTGACATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACt        <  1:611717/148‑1 (MQ=255)
                                                                                                                                          ggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTc       <  2:3168066/148‑1 (MQ=255)
                                                                                                                                           gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa      <  2:784807/148‑1 (MQ=255)
                                                                                                                                           gTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCa      <  2:2956035/148‑1 (MQ=255)
                                                                                                                                               aaTACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAg  >  1:3181246/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAG  >  CP009273/3692848‑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: