Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 355,222 A→C 6.9% intergenic (‑48/+37) cbcW ← / ← cbcX choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947355,2220AC 6.9% 474.4 / 19.0 145intergenic (‑48/+37)cbcW/cbcXcholine/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (68/67);  new base C (3/7);  total (71/74)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 159 total aligned reads displayed.

GCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGG  >  NC_002947/355073‑355371
                                                                                                                                                     |                                                                                                                                                     
gCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGggcaggca                                                                                                                                                       <  2:4582500/150‑1 (MQ=255)
   ttGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACgg                                                                                                                                                    >  1:4408232/1‑150 (MQ=255)
    tGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGt                                                                                                                                                   <  2:587378/150‑1 (MQ=255)
    tGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGt                                                                                                                                                   <  1:1770873/150‑1 (MQ=255)
     gcgcGTTGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTc                                                                                                                                                  <  1:4828941/150‑1 (MQ=255)
       gcgATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCt                                                                                                                                                <  2:341897/150‑1 (MQ=255)
         gATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGc                                                                                                                                              >  1:2885922/1‑150 (MQ=255)
         gATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGc                                                                                                                                              >  2:2308852/1‑150 (MQ=255)
            ggCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGcccc                                                                                                                                           <  1:4627589/150‑1 (MQ=255)
             gCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCg                                                                                                                                          <  1:1394908/150‑1 (MQ=255)
              cGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCgg                                                                                                                                         <  2:2931713/150‑2 (MQ=255)
              cGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCgg                                                                                                                                         <  2:181189/150‑2 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCggg                                                                                                                                        <  1:2100760/150‑3 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                        >  2:5217423/1‑150 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                        <  1:3750966/150‑1 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                        >  2:4004216/1‑150 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                        <  2:3901896/150‑1 (MQ=255)
               gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCggg                                                                                                                                        <  1:2692110/150‑3 (MQ=255)
                  gAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTa                                                                                                                                     >  1:5077042/1‑150 (MQ=255)
                         ttCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTAc                                                                                                                              <  1:3381574/150‑1 (MQ=255)
                         ttCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTAc                                                                                                                              <  1:3381225/150‑1 (MQ=255)
                          tCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACt                                                                                                                             <  2:3048198/150‑1 (MQ=255)
                             cgcCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCg                                                                                                                          <  1:769069/150‑1 (MQ=255)
                              gcCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGt                                                                                                                         <  1:5450728/150‑1 (MQ=255)
                              gcCGTTCTGGGTCAGCCATTCGACGAACGACGCGAGGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGt                                                                                                                         <  1:3693216/150‑1 (MQ=255)
                                       ggTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTa                                                                                                                <  1:4655031/150‑1 (MQ=255)
                                        gTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGCCACGGTCCTGCCCCg                                                                                                                                          >  1:5024635/1‑123 (MQ=255)
                                         tCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTAtt                                                                                                              <  2:2953722/150‑1 (MQ=255)
                                          cAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                             <  2:3041860/150‑1 (MQ=255)
                                          cAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                             <  2:3491449/150‑1 (MQ=255)
                                                ttCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAg                                                                                                       >  1:1238384/1‑150 (MQ=255)
                                                ttCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAg                                                                                                       >  2:2190988/1‑150 (MQ=255)
                                                   gacgaACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCtt                                                                                                    >  2:3309826/1‑150 (MQ=255)
                                                   gacgaACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCtt                                                                                                    <  1:2537593/150‑1 (MQ=255)
                                                   gacgaACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCtt                                                                                                    >  1:2191510/1‑150 (MQ=255)
                                                   gacgaACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCtt                                                                                                    >  1:4527671/1‑150 (MQ=255)
                                                     cgaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTgg                                                                                                  <  1:3739010/150‑1 (MQ=255)
                                                      gaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGc                                                                                                 >  2:3603679/1‑150 (MQ=255)
                                                        acgacgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCt                                                                                               >  1:3904491/1‑150 (MQ=255)
                                                          gacgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTg                                                                                             >  1:3084865/1‑150 (MQ=255)
                                                          gacgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTg                                                                                             <  1:3267133/150‑1 (MQ=255)
                                                           acgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGa                                                                                            <  1:551801/150‑1 (MQ=255)
                                                            cgcgATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGAc                                                                                           <  1:4955336/150‑1 (MQ=255)
                                                            cgcgATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGAc                                                                                           <  1:3309826/150‑1 (MQ=255)
                                                             gcgATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGAcc                                                                                          >  2:3286744/1‑150 (MQ=255)
                                                               gATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGc                                                                                        <  2:5077042/150‑1 (MQ=255)
                                                                aTGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGcc                                                                                       <  1:4151790/150‑1 (MQ=255)
                                                                    aCTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTcca                                                                                   >  1:1804924/1‑150 (MQ=255)
                                                                      tGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccagg                                                                                 >  1:1857337/1‑150 (MQ=255)
                                                                      tGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccagg                                                                                 >  2:834835/1‑150 (MQ=255)
                                                                           cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggccagg                                                                            >  2:3642248/1‑150 (MQ=255)
                                                                            aGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGt                                                                                                          <  1:1057487/119‑1 (MQ=255)
                                                                              gggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCtt                                                                         <  1:1100724/150‑1 (MQ=255)
                                                                              gggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCtt                                                                         <  2:4408232/150‑1 (MQ=255)
                                                                              gggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCtt                                                                         <  1:2580594/150‑1 (MQ=255)
                                                                                 tATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGcc                                                                      <  2:4793034/150‑1 (MQ=255)
                                                                                        tGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACg                                                               >  1:3950032/1‑150 (MQ=255)
                                                                                           tCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACggtg                                                            <  2:5055474/150‑1 (MQ=255)
                                                                                             gATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAGATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACggtggt                                                          >  1:5357608/1‑150 (MQ=255)
                                                                                             gATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACggtggt                                                          <  2:1718703/150‑1 (MQ=255)
                                                                                                aaGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAAc                                                       <  2:1383751/150‑1 (MQ=255)
                                                                                                  gCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGc                                                     >  1:5160320/1‑150 (MQ=255)
                                                                                                  gCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGCAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGc                                                     >  1:2432324/1‑150 (MQ=255)
                                                                                                   cATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACggtggt                                                          <  2:706612/144‑1 (MQ=255)
                                                                                                   cATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGcc                                                    >  1:2644029/1‑150 (MQ=255)
                                                                                                     tGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCgg                                                  >  1:652181/1‑150 (MQ=255)
                                                                                                       aTAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGgcc                                                >  2:45713/1‑150 (MQ=255)
                                                                                                               cTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGc                                        >  2:2473464/1‑150 (MQ=255)
                                                                                                                    aCCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCa                                   >  1:1446123/1‑150 (MQ=255)
                                                                                                                    aCCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCa                                   >  1:1268418/1‑150 (MQ=255)
                                                                                                                    aCCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCa                                   >  1:3016000/1‑150 (MQ=255)
                                                                                                                        gCAGAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCtt                               <  2:2985769/150‑1 (MQ=255)
                                                                                                                          aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGc                             >  2:5383492/1‑150 (MQ=255)
                                                                                                                            aGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCt                           >  1:420974/1‑150 (MQ=255)
                                                                                                                            aGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCt                           <  2:1855429/150‑1 (MQ=255)
                                                                                                                              aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                         >  1:3570762/1‑150 (MQ=255)
                                                                                                                              aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                         >  2:1898235/1‑150 (MQ=255)
                                                                                                                               ttGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                        >  1:3823462/1‑150 (MQ=255)
                                                                                                                               ttGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                        <  1:1856143/150‑1 (MQ=255)
                                                                                                                                   aaTCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGt                    <  2:4239660/150‑1 (MQ=255)
                                                                                                                                     tCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTCGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTcc                  >  2:3545068/1‑150 (MQ=255)
                                                                                                                                      cAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCt                 <  1:2308852/150‑1 (MQ=255)
                                                                                                                                       aaCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCa                                       >  1:4934040/1‑127 (MQ=255)
                                                                                                                                       aaCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCa                                       <  2:4934040/127‑1 (MQ=255)
                                                                                                                                        aCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTg               <  2:204165/150‑1 (MQ=255)
                                                                                                                                        aCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTg               >  1:3451226/1‑150 (MQ=255)
                                                                                                                                         cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa              >  1:849197/1‑150 (MQ=255)
                                                                                                                                         cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa              >  2:1487314/1‑150 (MQ=255)
                                                                                                                                         cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa              <  1:2500/150‑1 (MQ=255)
                                                                                                                                         cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa              <  2:3814835/150‑1 (MQ=255)
                                                                                                                                          tcGGGGCAGGACCGTGCCTGCCCCGGGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGac             >  1:1380741/2‑149 (MQ=255)
                                                                                                                                          ccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGAt             >  2:1830183/1‑150 (MQ=255)
                                                                                                                                          ccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGAt             >  1:4351228/1‑150 (MQ=255)
                                                                                                                                          ccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGAt             >  1:445505/1‑150 (MQ=255)
                                                                                                                                                caggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCAttt       <  1:2367048/150‑1 (MQ=255)
                                                                                                                                                 aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGccagccag                                          <  2:645006/114‑1 (MQ=255)
                                                                                                                                                   gcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCtt    <  1:5092327/150‑1 (MQ=255)
                                                                                                                                                    caCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCGACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTg   >  1:4739190/1‑150 (MQ=255)
                                                                                                                                                    caCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTg   >  2:1804615/1‑150 (MQ=255)
                                                                                                                                                     aCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTgg  <  2:492133/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGG  >  NC_002947/355073‑355371

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: