Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,103 T→C 5.4% Y52H (TAC→CAC)  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,1030TC 5.4% 223.7 / 2.2 75Y52H (TAC→CAC) kdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/2);  total (37/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 = 9.62e-01

CAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGCCTTAC  >  NC_002947/1807956‑1808249
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                                            <  6:44970/142‑1 (MQ=255)
              cgTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTc                                                                                                                                            >  5:44970/1‑142 (MQ=255)
                  ggTAACATCGAGATCTCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                                  <  8:104944/148‑1 (MQ=255)
                    tAACATCGAGATCGCCAACGACAAGCCGTTGGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                                >  1:29677/1‑148 (MQ=255)
                     aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                                  >  6:110934/1‑145 (MQ=255)
                     aaCATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGtt                                                                                                                                  <  5:110934/145‑1 (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                                                                                                                          <  8:274105/149‑1 (MQ=255)
                         tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                                          <  2:42994/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                                                                                                                         <  3:228549/115‑1 (MQ=255)
                                                            ggcggcATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAgg                                                                                                                         >  4:228549/1‑115 (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                         >  2:96271/1‑149 (MQ=255)
                                                                                                                           ggTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGa                         >  4:152848/1‑149 (MQ=255)
                                                                                                                            gTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt                        <  5:96890/149‑1 (MQ=255)
                                                                                                                                 cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                                               <  8:209187/89‑1 (MQ=255)
                                                                                                                                 cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                                               >  7:209187/1‑89 (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)
                                                                                                                                             taacccGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGc       <  5:121766/146‑1 (MQ=255)
                                                                                                                                             tATCCCGTACGTGTTCAAGGCCAGTTTCGCCAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGCGGGGCTGAAGATCTTCGAAGAGATCACGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGc       >  2:242656/1‑149 (MQ=255)
                                                                                                                                             tATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAg        >  8:220607/1‑148 (MQ=255)
                                                                                                                                                 ccGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGCCtt    <  2:130539/148‑1 (MQ=255)
                                                                                                                                                 ccGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGCCTTa   >  6:54678/1‑149 (MQ=255)
                                                                                                                                                  cGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTCGCATCGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGCCtttg  <  7:220607/149‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGC‑GATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCACGAGCCTTAC  >  NC_002947/1807956‑1808249

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