Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,098 T→A 6.7% 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.7% 221.6 / 4.7 74I50N (ATC→AAC) kdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (3/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGC  >  NC_002947/1807951‑1808244
                                                                                                                                                    |                                                                                                                                                  
tGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATc                                                                                                                                                   <  1:106028/149‑1 (MQ=255)
  aCCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACCTGCGGGTTACCGAGAAGCTCGGTATccc                                                                                                                                                 <  1:195750/149‑1 (MQ=255)
    gcAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                                                >  6:130791/2‑148 (MQ=255)
     cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTa                                                                                                                                              <  5:130791/149‑1 (MQ=255)
      agaagaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                                             >  3:308852/1‑149 (MQ=255)
       gaagaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGGTTGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                                             <  1:263391/149‑1 (MQ=255)
           aTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACgtgt                                                                                                                                         >  3:296541/1‑148 (MQ=255)
               ttCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAg                                                                                                                                    <  5:66921/149‑1 (MQ=255)
                   cgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTc                                                                                                                                       >  5:44970/1‑142 (MQ=255)
                   cgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTc                                                                                                                                       <  6:44970/142‑1 (MQ=255)
                       ggTAACATCGAGATCTCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                             <  8:104944/148‑1 (MQ=255)
                         tAACATCGAGATCGCCAACGACAAGCCGTTGGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                           >  1:29677/1‑148 (MQ=255)
                          aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                             <  5:110934/145‑1 (MQ=255)
                          aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                             >  6:110934/1‑145 (MQ=255)
                           aCATCGTGATCGCCTACGACAAGCCGTTAGTTCTGTTCTTCTTCATGAACGTTATTGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGATCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGa                                                                                                                         <  3:32869/148‑1 (MQ=255)
                           aCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGa                                                                                                                         <  2:29677/148‑1 (MQ=255)
                           aCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGa                                                                                                                         <  4:296541/148‑1 (MQ=255)
                             aTCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGAc                                                                                                                        >  2:327632/1‑147 (MQ=255)
                              tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                                     <  2:42994/149‑1 (MQ=255)
                              tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                                     <  8:274105/149‑1 (MQ=255)
                                       ccAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGt                                                                                                            >  6:96890/1‑149 (MQ=255)
                                         aaCGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                          >  7:322214/1‑149 (MQ=255)
                                          aCGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtcg                                                                                                         >  1:134399/1‑149 (MQ=255)
                                                 gCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTaa                                                                                                  >  8:76582/1‑149 (MQ=255)
                                                 gCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCACgta                                                                                                                                          >  6:151459/1‑108 (MQ=25)
                                                 gCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCACgta                                                                                                                                          <  5:151459/109‑2 (MQ=25)
                                                     ttCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTc                                                                                              <  1:327632/149‑1 (MQ=255)
                                                        gTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTa                                                                                           >  3:125991/1‑149 (MQ=255)
                                                              ttCGGCGTCATGAACGTTCTGGAATCCCGTGACCTGGC‑TATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCATACCGt                                                                                       <  3:135868/147‑1 (MQ=255)
                                                              ttCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTgg                                                                                     >  1:130539/1‑149 (MQ=255)
                                                                 ggcggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAgg                                                                                                                    >  4:228549/1‑115 (MQ=255)
                                                                 ggcggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAgg                                                                                                                    <  3:228549/115‑1 (MQ=255)
                                                                  gcggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                          <  6:258291/124‑1 (MQ=255)
                                                                  gcggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                          >  5:258291/1‑124 (MQ=255)
                                                                   cggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTg                                                                                <  4:308852/149‑1 (MQ=255)
                                                                        tGAACGTTCTGGAATCGCTTGAGCTGGC‑GATGAAGGTCTGCTAAGAGTACGTGCGGGTTATCGAGAAGCTCTGTATGCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCcttgcttg                                                                           <  4:125991/149‑1 (MQ=255)
                                                                          aaCGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                                          >  3:278669/1‑148 (MQ=255)
                                                                           aCGTTCTGTAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                                        <  2:140490/149‑1 (MQ=255)
                                                                            cGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCgg                                                                                                     >  6:285240/1‑119 (MQ=255)
                                                                            cGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCgg                                                                                                     <  5:285240/119‑1 (MQ=255)
                                                                              ttCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAagag                                                                     >  2:324501/1‑149 (MQ=255)
                                                                                cTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAgg                                                                    <  2:334812/148‑1 (MQ=255)
                                                                                   gAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCt                                                                >  8:147699/1‑149 (MQ=255)
                                                                                   gAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCt                                                                <  2:134399/149‑1 (MQ=255)
                                                                                      tCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTg                                                               >  4:121420/1‑147 (MQ=255)
                                                                                                tGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCg                                                     >  5:326721/1‑147 (MQ=255)
                                                                                                  gC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAaga                                                 <  4:278669/149‑1 (MQ=255)
                                                                                                     gATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATATAGGAagag                                                >  3:31761/1‑148 (MQ=255)
                                                                                                       tGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGAt                                              >  6:193581/1‑148 (MQ=255)
                                                                                                       tGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATc                                             <  5:193581/149‑1 (MQ=255)
                                                                                                             tCTGCGCAGAGTTAGTGCGGGGTACCGAGTATCTCGGTATTCCGGACGTGTTGAAGGGCATTTTGGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTTGCCCTGGCATGGTAGAGGGGCTGAAGATCTTCGAAGAGATCAAgccc                                       <  4:31761/149‑3 (MQ=255)
                                                                                                               tGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCa                                      >  7:174275/1‑148 (MQ=255)
                                                                                                                     gagTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCaa                               <  3:121420/149‑1 (MQ=255)
                                                                                                                      agTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCTGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAAc                              <  7:76582/149‑1 (MQ=255)
                                                                                                                        tACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGt                            >  1:329709/1‑149 (MQ=255)
                                                                                                                        caCGTACGGGATACCGAGCTTCTCGGTAACCCGCACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGt                            >  1:334495/2‑149 (MQ=38)
                                                                                                                          cGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGc                          <  1:324501/149‑1 (MQ=255)
                                                                                                                           gTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGcc                         >  6:160481/1‑149 (MQ=255)
                                                                                                                               gggATACCGAGCTTCTCGGTAACCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGt                      >  6:121766/1‑148 (MQ=255)
                                                                                                                                ggTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGa                    >  4:152848/1‑149 (MQ=255)
                                                                                                                                ggTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGa                    >  2:96271/1‑149 (MQ=255)
                                                                                                                                 gTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt                   <  5:96890/149‑1 (MQ=255)
                                                                                                                                      cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                                          >  7:209187/1‑89 (MQ=255)
                                                                                                                                      cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                                          <  8:209187/89‑1 (MQ=255)
                                                                                                                                      cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCAcc               >  2:25787/1‑148 (MQ=255)
                                                                                                                                      cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCg              <  6:326721/149‑1 (MQ=255)
                                                                                                                                      cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTCCCGGTGATCAcc               >  3:248956/1‑148 (MQ=255)
                                                                                                                                        agaAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAt            <  7:147699/149‑1 (MQ=255)
                                                                                                                                            gCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCa                 >  7:110298/1‑140 (MQ=255)
                                                                                                                                            gCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCa                 <  8:110298/140‑1 (MQ=255)
                                                                                                                                                ggTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGa    <  8:174275/149‑1 (MQ=255)
                                                                                                                                                ggTAACCCGCACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGa    <  2:334495/149‑1 (MQ=255)
                                                                                                                                                  tATCCCGTACGTGTTCAAGGCCAGTTTCGCCAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGCGGGGCTGAAGATCTTCGAAGAGATCACGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGc  >  2:242656/1‑149 (MQ=255)
                                                                                                                                                  tATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAg   >  8:220607/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGC  >  NC_002947/1807951‑1808244

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