Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,018,014 A→T 9.9% S632C (AGC→TGC)  bcsB → cellulose synthase and translocator subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,018,0140AT 9.9% 201.7 / 12.8 71S632C (AGC→TGC) bcsBcellulose synthase and translocator subunit
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (3/4);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

ACTTGCTGCTGATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTGGCGATGTAGGC  >  NC_002947/3017865‑3018163
                                                                                                                                                     |                                                                                                                                                     
aCTTGCTGCTGATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCaa                                                                                                                                                       <  1:268462/150‑1 (MQ=255)
   tgctgcTGATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGTCCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCCGGGCCGCAGCTCATCCTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAagca                                                                                                                                                    <  2:108931/150‑1 (MQ=255)
     ctgctgATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCTTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGTAagcagc                                                                                                                                                  <  2:707589/150‑1 (MQ=255)
      tgctgATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGcc                                                                                                                                                 >  1:1200773/1‑150 (MQ=255)
      tgctgATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGcc                                                                                                                                                 >  2:496589/1‑150 (MQ=255)
      tgctgATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGcc                                                                                                                                                 >  1:231760/1‑150 (MQ=255)
       gctgATCGGCAGCCTGCCCGAGGCCCTGCGTGACGCGCCCGACCTTGGCCTGTTGCTCATCGCCCAGCGTGACTGGTTGCTGCAGGGACGCAGCGCATCGTTGCCCGGTGGCCAGGGTTTCGATACTGAGCCAGTGGCTGTAAGCAGCcg                                                                                                                                                >  2:627963/1‑150 (MQ=255)
        ctgATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCcgc                                                                                                                                               >  2:257741/1‑150 (MQ=255)
          gATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccaccg                                                                                                                                             <  2:1184106/150‑5 (MQ=255)
           aTCGGCAGCCTGCCCGAGGCCCTGAACGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTg                                                                                                                                            >  1:1261403/1‑150 (MQ=255)
              ggCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCg                                                                                                                                         >  1:154996/1‑150 (MQ=255)
                   ccTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAg                                                                                                                                    <  1:951522/150‑1 (MQ=255)
                        ccGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               >  2:706961/1‑132 (MQ=255)
                        ccGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               <  1:706961/134‑3 (MQ=255)
                        ccGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               <  1:706568/134‑3 (MQ=255)
                        ccGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               >  2:706568/1‑132 (MQ=255)
                        ccGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGccc                                                                                                                               >  1:327504/1‑150 (MQ=255)
                             cgaacccgcgACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGgcg                                                                                                                          >  2:23328/7‑150 (MQ=255)
                                 tGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGa                                                                                                                      <  2:614148/150‑1 (MQ=255)
                                      gACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGCCTGGTTGCTGAAGGGCCGCAGCGCATCGTTGCCCGGTGGACAGCGTTTCGATACTGTGCAGGGTGCTTGTAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCg                                                                                                                 <  1:739463/150‑1 (MQ=255)
                                      gACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCg                                                                                                                 >  1:243888/1‑150 (MQ=255)
                                      gACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCg                                                                                                                 <  1:1215998/150‑1 (MQ=255)
                                       aCGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCgg                                                                                                                <  1:1235133/150‑1 (MQ=255)
                                           gcCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCAt                                                                                                            <  2:1122283/150‑1 (MQ=255)
                                           gcCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCAt                                                                                                            <  1:329794/150‑1 (MQ=255)
                                            cccGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATc                                                                                                           >  2:1331287/1‑150 (MQ=255)
                                              cGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGAGGCCATCAc                                                                                                         >  2:189242/1‑150 (MQ=255)
                                                aCCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCg                                                                                                       <  2:1300449/150‑1 (MQ=255)
                                                 ccTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCgg                                                                                                      <  1:1340468/150‑1 (MQ=255)
                                                       ccTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGaa                                                                                                <  2:327504/150‑1 (MQ=255)
                                                             gCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGcc                                                                                          <  2:1200773/150‑1 (MQ=255)
                                                             gCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGcc                                                                                          >  2:905353/1‑150 (MQ=255)
                                                              cTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCg                                                                                         <  2:231760/150‑1 (MQ=255)
                                                              cTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGAGTGCGGCCATCACCGGGCTGAAGTCGCCg                                                                                         >  1:1026648/1‑150 (MQ=255)
                                                               tCAGCGCCCAGCGTGACTGTTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGt                                                                                        >  2:59789/1‑150 (MQ=255)
                                                               tCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGt                                                                                        >  1:294234/1‑150 (MQ=255)
                                                               tCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGAGGAACGGCGGAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGCCCATCACCGGGCTGAAGTCGCCGt                                                                                        >  1:798552/1‑150 (MQ=255)
                                                                 aGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTc                                                                                      >  2:1195355/1‑150 (MQ=255)
                                                                 aGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTc                                                                                      <  2:264916/150‑1 (MQ=255)
                                                                  gcgcCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTcc                                                                                     <  1:1209365/150‑1 (MQ=255)
                                                                      ccAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGa                                                                                 <  1:289019/150‑1 (MQ=255)
                                                                      ccAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGa                                                                                 >  2:793537/1‑150 (MQ=255)
                                                                           gTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGNGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGagcagc                                                                            <  1:627963/150‑1 (MQ=255)
                                                                           gTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGagcagc                                                                            <  1:496589/150‑1 (MQ=255)
                                                                                tGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCa                                                                                                          <  2:355822/115‑1 (MQ=255)
                                                                                tGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCa                                                                                                          >  1:355822/1‑115 (MQ=255)
                                                                                tGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAgt                                                                       <  2:392720/150‑1 (MQ=255)
                                                                                tGGTTGCGGCAGGGCCGCAGCGCATCGTTACCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAgt                                                                       <  2:174676/150‑1 (MQ=255)
                                                                                  gTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCgg                                                                                                                <  1:781898/107‑1 (MQ=255)
                                                                                  gTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCgg                                                                                                                >  2:781898/1‑107 (MQ=255)
                                                                                       tgcAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               <  1:803448/71‑3 (MQ=255)
                                                                                       tgcAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCTTGCAGccac                                                                                                                                               >  2:803448/1‑69 (MQ=255)
                                                                                       tgcAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTgg                                                                >  2:1132277/1‑150 (MQ=255)
                                                                                          aGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGa                                                                                 <  2:1262669/130‑1 (MQ=255)
                                                                                          aGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGa                                                                                 >  1:1262669/1‑130 (MQ=255)
                                                                                           gggCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCtt                                                            <  1:793537/150‑1 (MQ=255)
                                                                                              ccGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCt                                                         <  1:1298268/150‑1 (MQ=255)
                                                                                              ccGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAACCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGACGAACAACGCGCCGCGCCAGGATGTGGGCTTGCt                                                         <  2:798417/150‑1 (MQ=255)
                                                                                               cGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCttagtg                                                        >  2:366723/1‑146 (MQ=255)
                                                                                                   gcgcATCGTTGCCCGGTGGCCAGCGTTTCGATCCTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCCCCCAGGCGCATGATGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCa                                                    >  1:1033180/1‑150 (MQ=255)
                                                                                                     gcATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAac                                                  <  1:59789/150‑1 (MQ=255)
                                                                                                     gcATCGTTGCCCGGTGGCCAGCGTTTCGATACTCAGCCGGCGTCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATCGCGGCCATCACCAGGCCGAAGTCGCCGTTCCACGAGCAGCCCAGTGTTCTGGCCTTGCTGGCCacc                                                  >  2:1007390/1‑148 (MQ=255)
                                                                                                      cATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAaca                                                 <  1:189242/150‑1 (MQ=255)
                                                                                                       aTCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGGAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAacag                                                <  2:1026648/150‑1 (MQ=255)
                                                                                                        tCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAacagc                                               >  1:1264275/1‑150 (MQ=255)
                                                                                                                cGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCAcc                                                                                                        <  2:86160/85‑1 (MQ=255)
                                                                                                                cGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCAcc                                                                                                        >  1:86160/1‑85 (MQ=255)
                                                                                                                       cAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCa                                >  2:958790/1‑150 (MQ=255)
                                                                                                                       cAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCa                                >  1:1031959/1‑150 (MQ=255)
                                                                                                                                 atcgtGATCCGGTAGCTGAAAGCAGTCGCGTGGCGGTCAGCGCGCGGGGGCCGATTGCGCCCGTCACCGCTCTGAAGTCGCCGTTCCACGAGCAGCGCAATGTGATAGCCTTGCTGCCCCACAGCGACCGCGGTTACGCATTGTTgcgct                      >  2:296275/5‑149 (MQ=255)
                                                                                                                                         ccGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTCTCTTTTTTTTGTTCAACAGTGACGGCGATTACGCATTGCTGCGCGACATGCtt              <  1:958790/150‑1 (MQ=255)
                                                                                                                                               gCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGa                     >  2:191984/1‑137 (MQ=255)
                                                                                                                                               gCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGa                     <  1:191984/137‑1 (MQ=255)
                                                                                                                                                 tGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTGGCGATGt      >  1:19095/1‑150 (MQ=255)
                                                                                                                                                    aaGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTGGCGATGTAgg   >  2:361782/1‑150 (MQ=255)
                                                                                                                                                     agcagcCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTGGCGATGTAGGc  <  2:1031959/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTTGCTGCTGATCGGCAGCCTGCCCGAGGCCCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTGGCGATGTAGGC  >  NC_002947/3017865‑3018163

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: