Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,098 T→A 6.2% I50N (ATC→AAC)  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,0980TA 6.2% 195.7 / 2.8 65I50N (ATC→AAC) kdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base A (2/2);  total (27/38)
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.35e-01

ACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGT  >  NC_002947/1807953‑1808236
                                                                                                                                                 |                                                                                                                                            
aCCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGGGCGGGTTACCGAGAAGCTCGGTATcc                                                                                                                                            <  1:216317/148‑1 (MQ=255)
 cccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                          >  7:251374/1‑149 (MQ=255)
  ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                          >  5:165726/1‑148 (MQ=255)
  ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                          >  7:342764/1‑148 (MQ=255)
  ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                          >  2:282023/1‑148 (MQ=255)
  ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGt                                                                                                                                         <  3:72207/149‑1 (MQ=255)
   cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTa                                                                                                                                        <  6:165726/149‑1 (MQ=255)
    agaagaGCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACCTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCAc                                                                                                                                       <  6:1449/149‑1 (MQ=38)
     gaagaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACg                                                                                                                                      >  7:338152/1‑149 (MQ=255)
               cgcgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCTGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAgg                                                                                                                             <  6:62896/148‑1 (MQ=255)
                gcgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGc                                                                                                                            <  4:303140/148‑1 (MQ=255)
                 cgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCa                                                                                                                          <  4:19969/149‑1 (MQ=255)
                     ggTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGttt                                                                                                                      <  8:342764/149‑1 (MQ=255)
                      gTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                     <  7:192605/149‑1 (MQ=255)
                         aCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGa                                                                                                                   <  6:22479/148‑1 (MQ=255)
                           aTCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑c                                                                                                                >  2:143836/1‑147 (MQ=255)
                           aTCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑Caa                                                                                                              <  5:127299/149‑1 (MQ=255)
                            tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTc                                                                                                                                 >  6:150877/1‑131 (MQ=255)
                            tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTc                                                                                                                                 <  5:150877/131‑1 (MQ=255)
                            tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAg                                                                                                             >  5:243683/1‑149 (MQ=255)
                              gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAgg                                                                                                            <  2:166663/148‑1 (MQ=255)
                              gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGc                                                                                                           >  3:150152/1‑149 (MQ=255)
                                     ccAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑Caa                                                                                                              >  6:2634/1‑139 (MQ=255)
                                     ccAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑Caa                                                                                                              <  5:2634/139‑1 (MQ=255)
                                           aCAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTtcgtcg                                                                                              <  8:251374/149‑1 (MQ=255)
                                             aaGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCgg                                                                                             >  2:101169/1‑148 (MQ=255)
                                                ccGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAAc                                                                                         >  4:259006/1‑149 (MQ=255)
                                                 cGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACg                                                                                        <  1:282023/149‑1 (MQ=255)
                                                       tCCTGTGCGGCGGCATGACCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCGGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTAc                                                                                  <  8:175064/149‑1 (MQ=255)
                                                         cTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCg                                                                                >  7:274992/1‑149 (MQ=255)
                                                            ttCGGCGGCATGAACGTTCTGTATTGCCTTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAATCTCGGTATCCCGTACGTGTTCAAGGCGAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTg                                                                              <  7:167929/148‑1 (MQ=255)
                                                            ttCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGt                                                                               <  1:174297/147‑1 (MQ=255)
                                                                gcggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCt                                                                         <  1:143836/149‑1 (MQ=255)
                                                                    cATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCa                                                                     >  7:251587/1‑149 (MQ=255)
                                                                           gTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGAGGGTAGGCAAACCGTTCTTCGGTACCGACATACCGTGGCCCTGCCATGGa                                                                 >  1:147187/1‑148 (MQ=255)
                                                                           gTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                                >  4:321801/1‑147 (MQ=255)
                                                                            ttCTTGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAagag                                                             <  6:117991/149‑1 (MQ=255)
                                                                                ggAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAACAGGGGc                                                         >  3:347831/1‑149 (MQ=255)
                                                                                 gAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGc                                                         <  5:246743/148‑1 (MQ=255)
                                                                                   aTCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGa                                                      <  3:259006/149‑1 (MQ=255)
                                                                                      ccGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCACgtac                                                                                                                                   <  8:296334/71‑3 (MQ=14)
                                                                                      ccGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCACgtac                                                                                                                                   >  7:296334/1‑69 (MQ=14)
                                                                                           aCCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCtt                                               <  8:251587/148‑1 (MQ=255)
                                                                                           aCCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCtt                                               <  4:150152/148‑1 (MQ=255)
                                                                                                gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                           <  1:101169/147‑1 (MQ=255)
                                                                                                gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                           <  8:199896/147‑1 (MQ=255)
                                                                                                gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                           >  7:199896/1‑147 (MQ=255)
                                                                                                gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                           <  4:238917/147‑1 (MQ=255)
                                                                                                gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                           <  8:274992/147‑1 (MQ=255)
                                                                                                   aTGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAggg                                                           >  6:173459/1‑128 (MQ=255)
                                                                                                   aTGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAggg                                                           <  5:173459/128‑1 (MQ=255)
                                                                                                       aGGTCTTCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAg                                  <  6:75149/149‑1 (MQ=255)
                                                                                                               gAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGa                                                   <  7:263443/124‑1 (MQ=255)
                                                                                                               gAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGa                                                   >  8:263443/1‑124 (MQ=255)
                                                                                                               gAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCtt                          >  1:29910/1‑149 (MQ=255)
                                                                                                                        gTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGcc                 >  2:210195/1‑149 (MQ=255)
                                                                                                                           cGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGt              <  6:243683/149‑1 (MQ=255)
                                                                                                                            gggTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTg             <  8:338152/149‑1 (MQ=255)
                                                                                                                             ggatACCGAGCTTCTCGGTAACCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGa            >  6:395782/4‑149 (MQ=255)
                                                                                                                              gTTACCGTGATGCTGGGTATCCCGTACGTGTTCATGGCCAGTTTCGA‑‑CTAGGCCAACCGTTCGTCGGTATCGTCCTAGCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt           <  4:347831/149‑1 (MQ=255)
                                                                                                                                   cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCAcc       >  8:14185/1‑148 (MQ=255)
                                                                                                                                   cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCAcc       >  2:25954/1‑148 (MQ=255)
                                                                                                                                      gaAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAtg   >  5:390415/1‑149 (MQ=255)
                                                                                                                                       aaGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAATGCCAACCGTTCGTCGGTAACGTCATACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAtgt  <  2:368151/149‑1 (MQ=255)
                                                                                                                                        aGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAtgt  <  4:176396/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
ACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGA‑‑CAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGT  >  NC_002947/1807953‑1808236

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