Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,089 A→T 6.5% K47M (AAG→ATG) ‡ kdsA‑I → 2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,808,0890AT 6.5% 182.8 / 2.7 62K47M (AAG→ATG) ‡kdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

CTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAT  >  NC_002947/1807943‑1808234
                                                                                                                                                  |                                                                                                                                                   
cTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGc                                                                                                                                                   >  3:32355/1‑149 (MQ=255)
 tGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCt                                                                                                                                                  >  5:199917/1‑149 (MQ=255)
   aTCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                >  4:67701/1‑149 (MQ=255)
   aTCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                <  6:199917/149‑1 (MQ=255)
           cccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                        <  5:301342/149‑1 (MQ=255)
            ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATccc                                                                                                                                         >  1:316292/1‑147 (MQ=255)
            ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATccc                                                                                                                                         <  2:316292/147‑1 (MQ=255)
            ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                        >  8:129178/1‑148 (MQ=255)
             cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTa                                                                                                                                      <  4:32355/149‑1 (MQ=255)
              agaagaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCAc                                                                                                                                     <  7:45531/149‑1 (MQ=39)
                      aTTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                <  1:134794/130‑1 (MQ=255)
                      aTTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                >  2:134794/1‑130 (MQ=255)
                      aTTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTc                                                                                                                              >  1:83919/1‑147 (MQ=255)
                       ttCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAg                                                                                                                           <  7:125096/149‑1 (MQ=255)
                          gcgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGc                                                                                                                         >  8:51746/1‑148 (MQ=255)
                           cgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCa                                                                                                                       <  7:194865/149‑1 (MQ=255)
                                 tAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCg                                                                                                                 <  2:46500/149‑1 (MQ=255)
                                  aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCaa                                                                                                                            >  6:351732/1‑137 (MQ=255)
                                  aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCaa                                                                                                                            <  5:351732/137‑1 (MQ=255)
                                   aCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAg                                                                                                                           <  4:166578/137‑1 (MQ=255)
                                   aCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAg                                                                                                                           >  3:166578/1‑137 (MQ=255)
                                      tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑g                                                                                                           <  2:320424/149‑1 (MQ=255)
                                        gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACgt‑gt                                                                                                                                >  4:152047/1‑127 (MQ=255)
                                        gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACgt‑gt                                                                                                                                <  3:152047/127‑1 (MQ=255)
                                         agaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGc                                                                                                         >  5:260435/1‑148 (MQ=255)
                                            tCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAAc                                                                                                     >  3:164868/1‑149 (MQ=255)
                                             cGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAAcc                                                                                                    >  3:257105/1‑149 (MQ=255)
                                              gCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCg                                                                                                   >  2:23032/1‑149 (MQ=255)
                                               ccAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGt                                                                                                  >  1:139490/1‑149 (MQ=255)
                                                 aaCGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTtc                                                                                                <  7:51746/149‑1 (MQ=255)
                                                    gACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTtcgtc                                                                                             <  1:23032/149‑1 (MQ=255)
                                                     aCAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCt                                                                                                                                                  >  4:30119/1‑97 (MQ=255)
                                                     aCAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCt                                                                                                                                                  <  3:30119/97‑1 (MQ=255)
                                                         gCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTaa                                                                                        >  7:107510/1‑149 (MQ=255)
                                                            gTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAgg                                                                                                                          <  5:298460/113‑1 (MQ=255)
                                                            gTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAgg                                                                                                                          >  6:298460/1‑113 (MQ=255)
                                                               cGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTcc                                                                                   <  3:67701/148‑1 (MQ=255)
                                                                   cTGTTCGGCGGCATGAACGTTCTGGAATCCGGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCg                                                                              >  1:30544/1‑149 (MQ=255)
                                                                      ttCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTgg                                                                           <  4:164868/149‑1 (MQ=255)
                                                                            ggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGt                                                                                                  <  3:242812/120‑1 (MQ=255)
                                                                            ggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGt                                                                                                  >  4:242812/1‑120 (MQ=255)
                                                                                   aCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                              <  2:83919/149‑1 (MQ=255)
                                                                                    cGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAACGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                              <  1:252760/149‑1 (MQ=255)
                                                                                      ttCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAagag                                                           <  8:107510/149‑1 (MQ=255)
                                                                                         tGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAgggg                                                        >  8:156683/1‑149 (MQ=255)
                                                                                           gAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAgggg                                                        <  5:61362/147‑1 (MQ=255)
                                                                                                ccGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTg                                                     >  5:307614/1‑145 (MQ=255)
                                                                                                 cGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGAt                                                >  7:2109/1‑149 (MQ=255)
                                                                                                  gTGACCTGGTGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATc                                               >  5:77091/1‑149 (MQ=255)
                                                                                                      ccTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                            <  6:307614/148‑1 (MQ=255)
                                                                                                          gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                         <  3:142626/147‑1 (MQ=255)
                                                                                                          gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                         <  8:2109/147‑1 (MQ=255)
                                                                                                              tGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATc                                   <  4:257105/149‑1 (MQ=255)
                                                                                                                 aGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAg                                <  6:260435/149‑1 (MQ=255)
                                                                                                                        cGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCt                         >  5:332019/1‑149 (MQ=255)
                                                                                                                        cGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCt                         >  8:73823/1‑149 (MQ=255)
                                                                                                                          ttgaacaCGTACGGGATACCGAGCTTCTCGGTAACCCGCACGTAGTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                            >  1:67409/7‑129 (MQ=25)
                                                                                                                          ttgaacaCGTACGGGATACCGAGCTTCTCGGTAACCCGCACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                            <  2:67409/122‑1 (MQ=25)
                                                                                                                                  gTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCCAGCGCACCTTCAACGTGcc               >  6:351497/1‑149 (MQ=255)
                                                                                                                                        gatACCGAGCTTCTCGGTAACCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt         >  7:47147/3‑149 (MQ=255)
                                                                                                                                              gagaAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGa   <  6:332019/149‑1 (MQ=255)
                                                                                                                                               agaAGCTCGGTATCCGGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAt  <  2:30544/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGT‑GTTCAAGGCCAGTTTCGACAA‑GGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAT  >  NC_002947/1807943‑1808234

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