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,9940TA95.3% 677.9 / 17.0 213F474Y (TTC→TAC) bcsGDUF3260 family inner membrane protein associated with cellulose production
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base A (114/89);  total (116/97)
Fisher's exact test for biased strand distribution p-value = 4.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 216 total aligned reads displayed.

ACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAGT  >  CP009273/3692849‑3693139
                                                                                                                                                 |                                                                                                                                                 
aCTGTACTTGAGAAATCGGGTCGTAATGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTATATCACCTTCGTCCCCGTTGGGGTGAAAttct                                                                                                                                                 <  1:1243264/148‑1 (MQ=255)
   gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCg                                                                                                                                              >  1:59514/1‑148 (MQ=255)
    aaCTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCgg                                                                                                                                             <  1:409120/148‑1 (MQ=255)
        tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                         <  2:3559759/148‑1 (MQ=255)
        tGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATg                                                                                                                                         >  2:1246783/1‑148 (MQ=255)
          agaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGaa                                                                                                                                       <  2:824750/148‑1 (MQ=255)
           gaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGTCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAg                                                                                                                                      <  2:2799681/148‑1 (MQ=255)
           gaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAg                                                                                                                                      <  1:2384559/148‑1 (MQ=255)
           gaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAg                                                                                                                                      <  2:1374778/148‑1 (MQ=255)
            aaaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAgg                                                                                                                                     <  1:2427905/148‑1 (MQ=255)
             aaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGc                                                                                                                                    >  2:2909640/1‑148 (MQ=255)
              aTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCa                                                                                                                                   <  2:2612250/148‑1 (MQ=255)
               tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                  >  2:656876/1‑148 (MQ=255)
               tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAc                                                                                                                                  >  2:656834/1‑148 (MQ=255)
                cGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAcc                                                                                                                                 >  2:1779437/1‑148 (MQ=255)
                cGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCAcc                                                                                                                                 >  2:2526571/1‑148 (MQ=255)
                 gggTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCg                                                                                                                                <  2:2042352/148‑1 (MQ=255)
                  ggTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGc                                                                                                                               >  1:1714746/1‑148 (MQ=255)
                      gTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCa                                                                                                                           <  2:3256500/148‑1 (MQ=255)
                        aaaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgg                                                                                                                         <  1:3438630/148‑1 (MQ=255)
                        aaaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgg                                                                                                                         >  1:3439966/1‑148 (MQ=255)
                         aaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggg                                                                                                                        >  1:2548808/1‑148 (MQ=255)
                         aaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggg                                                                                                                        >  1:3363211/1‑148 (MQ=255)
                          aGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgggg                                                                                                                       >  2:1157516/1‑148 (MQ=255)
                          aGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgggg                                                                                                                       >  2:1043270/1‑148 (MQ=255)
                          aGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAgggg                                                                                                                       >  2:1157688/1‑148 (MQ=255)
                           gTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAggggg                                                                                                                      <  2:3509588/148‑1 (MQ=255)
                           gTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAggggg                                                                                                                      >  1:3069618/1‑148 (MQ=255)
                            tGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGc                                                                                                                     >  2:3331581/1‑148 (MQ=255)
                            tGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGc                                                                                                                     <  2:1574057/148‑1 (MQ=255)
                                 gtggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGa                                                                                                                >  2:1955899/1‑148 (MQ=255)
                                 gtggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGa                                                                                                                <  2:860921/148‑1 (MQ=255)
                                   ggtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGAtt                                                                                                              <  2:2610745/148‑1 (MQ=255)
                                    gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTg                                                                                                             >  1:2180132/1‑148 (MQ=255)
                                     tCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGt                                                                                                            >  1:2942850/1‑148 (MQ=255)
                                          gtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCg                                                                                                       >  2:2967070/1‑148 (MQ=255)
                                          gtgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATAGTCATCg                                                                                                       >  1:2294781/1‑148 (MQ=255)
                                           tgCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGa                                                                                                      <  2:1352049/148‑1 (MQ=255)
                                            gCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaa                                                                                                     <  1:1649111/148‑1 (MQ=255)
                                              cGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaaca                                                                                                   <  1:276263/148‑1 (MQ=255)
                                                gAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGaacaac                                                                                                 >  1:206669/1‑148 (MQ=255)
                                                 aaCACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAAGAAcc                                                                                                >  2:3025560/1‑148 (MQ=255)
                                                    acGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATAGAACAACCGag                                                                                             <  1:1804188/148‑1 (MQ=255)
                                                      ggcggcGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagca                                                                                           >  1:3306094/1‑148 (MQ=255)
                                                      ggcggcGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagca                                                                                           >  1:1520825/1‑148 (MQ=255)
                                                        cggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGagcagc                                                                                         <  2:1490854/148‑1 (MQ=255)
                                                         ggcgCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCt                                                                                        <  1:3026555/148‑1 (MQ=255)
                                                           cgcgcTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTc                                                                                      >  2:2338052/1‑148 (MQ=255)
                                                                 gAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCt                                                                                <  1:885354/148‑1 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                             >  1:755071/1‑148 (MQ=255)
                                                                    gggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtc                                                                             >  1:3015065/1‑148 (MQ=255)
                                                                     ggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtct                                                                            <  2:3427398/148‑1 (MQ=255)
                                                                     ggCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTAtct                                                                            <  1:1043270/148‑1 (MQ=255)
                                                                       cGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTcc                                                                          >  2:184117/1‑148 (MQ=255)
                                                                       cGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAAGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTcc                                                                          <  2:2113553/148‑1 (MQ=255)
                                                                                gCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtg                                                                 >  1:2691389/1‑148 (MQ=255)
                                                                                gCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtg                                                                 >  2:1624371/1‑148 (MQ=255)
                                                                                 cAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCtggtgg                                                                >  2:1542317/1‑148 (MQ=255)
                                                                                     tATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcg                                                            >  2:2160591/1‑148 (MQ=255)
                                                                                      aTCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgc                                                           <  1:3118815/148‑1 (MQ=255)
                                                                                       tCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgct                                                          <  1:1165459/148‑2 (MQ=255)
                                                                                          ggCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTcgctatc                                                       >  2:2444793/1‑144 (MQ=255)
                                                                                          ggCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtc                                                       >  2:2564896/1‑148 (MQ=255)
                                                                                           gCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTtct                                                      <  1:956054/148‑1 (MQ=255)
                                                                                           gCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTtct                                                      >  1:1759461/1‑148 (MQ=255)
                                                                                              tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                   <  1:3419952/148‑1 (MQ=255)
                                                                                              tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGa                                                   <  2:484163/148‑1 (MQ=255)
                                                                                              tACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGa                                                   >  2:1023746/1‑148 (MQ=255)
                                                                                               aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGAt                                                  >  2:3319228/1‑148 (MQ=255)
                                                                                               aCGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGAt                                                  >  2:337827/1‑148 (MQ=255)
                                                                                                    atatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCaa                                             >  2:1191269/1‑148 (MQ=255)
                                                                                                    atatCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCaa                                             >  2:1968483/1‑148 (MQ=255)
                                                                                                      atCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGa                                           >  1:1567957/1‑148 (MQ=255)
                                                                                                       tCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAt                                          >  1:2536957/1‑148 (MQ=255)
                                                                                                       tCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAt                                          <  2:655218/148‑1 (MQ=255)
                                                                                                        cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtt                                         >  1:1574357/1‑148 (MQ=255)
                                                                                                        cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtt                                         <  1:985961/148‑1 (MQ=255)
                                                                                                        cccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAtt                                         >  2:317013/1‑148 (MQ=255)
                                                                                                         ccTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGAttt                                        >  1:3295349/1‑148 (MQ=255)
                                                                                                          cTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGAtttt                                       >  1:2220227/1‑148 (MQ=255)
                                                                                                           tAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTc                                      >  2:1285309/1‑148 (MQ=255)
                                                                                                             gCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAc                                    <  1:502758/148‑1 (MQ=255)
                                                                                                             gCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAc                                    >  1:3180220/1‑148 (MQ=255)
                                                                                                             gCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAc                                    >  2:42465/1‑148 (MQ=255)
                                                                                                              cccGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCAcc                                   >  1:774988/1‑148 (MQ=255)
                                                                                                               ccGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCg                                  >  1:475756/1‑148 (MQ=255)
                                                                                                               ccGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCg                                  >  2:2710875/1‑148 (MQ=255)
                                                                                                                cGTCTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGa                                 <  1:1285309/148‑1 (MQ=255)
                                                                                                                   cTATCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGa                              <  2:59514/148‑1 (MQ=255)
                                                                                                                      tCACCGACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACaa                           <  1:1922414/148‑1 (MQ=255)
                                                                                                                           gACGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTg                      <  1:184605/148‑1 (MQ=255)
                                                                                                                            aCGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGa                     <  1:1976208/148‑1 (MQ=255)
                                                                                                                            aCGTCCCCGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGa                     <  2:1981521/148‑1 (MQ=255)
                                                                                                                                ccccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTgg                 >  2:380766/1‑148 (MQ=255)
                                                                                                                                 cccGTTGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGa                >  1:923984/1‑148 (MQ=255)
                                                                                                                                     ttGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa            <  2:2929955/148‑1 (MQ=255)
                                                                                                                                     ttGGGGTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGaaaaa            <  2:2666665/148‑1 (MQ=255)
                                                                                                                                        gggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACt         <  1:1483218/148‑1 (MQ=255)
                                                                                                                                        gggTGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACt         <  2:351988/148‑1 (MQ=255)
                                                                                                                                           tGAAATACTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCTTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCAc      <  1:1432264/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTGTCATCGAACAACCGAGCAGCTTCCTGGCTATCTCCGATCTGGTGGTTCGCGTTCTCGATGGCAAGATTTTCACCGAAGACAATGTTGACTGGAAAAAACTCACCAGT  >  CP009273/3692849‑3693139

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