Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 262,944 C→G 6.6% intergenic (+179/‑306) PP_0211 → / → gabD‑I hypothetical protein/succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947262,9440CG 6.6% 232.5 / 5.9 76intergenic (+179/‑306)PP_0211/gabD‑Ihypothetical protein/succinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (42/29);  new base G (3/2);  total (45/31)
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.95e-01

AAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCCACAACACAA  >  NC_002947/262799‑263090
                                                                                                                                                 |                                                                                                                                                   
aaGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAgg                                                                                                                                                  >  1:222484/1‑149 (MQ=255)
  gCAGGACGAAAGCCTAGACAACGAACTGGCCGACGCTAACCTGATCATTGGCGACACCAAGAGTGTGGGTACAGGGCGCAACTGCGCAAAGTTCCACCGATTACGATGCCGATGAAAACAGCATCATCATTTCAGCACCGCAACAgtcc                                                                                                                                                >  6:393897/1‑146 (MQ=255)
    aGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGt                                                                                                                                              >  6:136401/1‑149 (MQ=255)
      gacgacCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgc                                                                                                                                           <  1:155380/149‑1 (MQ=255)
                     aaCGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgct                                                                                                                                          >  1:66032/1‑135 (MQ=255)
                     aaCGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgct                                                                                                                                          <  2:66032/135‑1 (MQ=255)
                        gAACTGGCCGAGGGTAACCTGCTCTTTGGAGACACACCGAGTGTGGGTACAGGGCGCAACTGCGCAACGAGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCACCACGGCATCAGGGCGTGATGAGATGATCGACCTCGc                                                                                                                         >  8:189149/1‑149 (MQ=255)
                          aCTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGcc                                                                                                                        >  1:117187/1‑148 (MQ=255)
                           cTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGTTGAGCTGATCGACCTCGCCaa                                                                                                                      <  6:53909/149‑1 (MQ=255)
                                cGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGt                                                                                                                 >  7:146006/1‑149 (MQ=255)
                                 gAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGtt                                                                                                                <  5:136401/149‑1 (MQ=255)
                                   gggTTACCTGCTCTTTGGCGACACCCCTAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCt                                                                                                              <  7:189149/149‑1 (MQ=255)
                                     gTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGa                                                                                                            >  3:326161/1‑149 (MQ=255)
                                      tAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGGATCATCGTTTCAGCACGGCATCAGGGCGTAGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGa                                                                                                            <  1:254797/149‑1 (MQ=255)
                                        aCCTGCTCTTTGGAGACACCCCGATTGTGGGTACTCTGCGCAACTTGGCAACTTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCATCACGGCATCAGGGCTTGCTTAGCTGATCGACCTCGCCAAGGCGTTCTGAtt                                                                                                          <  5:393897/148‑1 (MQ=255)
                                         ccTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgctgagctga                                                                                                                                   <  6:121160/122‑1 (MQ=255)
                                         ccTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgctgagctga                                                                                                                                   >  5:121160/1‑122 (MQ=255)
                                          cTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCAt                                                                                                       <  8:273231/149‑1 (MQ=255)
                                                tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATccc                                                                                                    <  4:326161/146‑1 (MQ=255)
                                                tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                                 >  5:210152/1‑149 (MQ=255)
                                                tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                                 >  8:12788/1‑149 (MQ=255)
                                                tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                                 >  6:390227/1‑149 (MQ=255)
                                                 ttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTcc                                                                                                >  6:116340/1‑149 (MQ=255)
                                                  tGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCg                                                                                               >  8:273009/1‑149 (MQ=255)
                                                   ggCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGTATCATCCTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCg                                                                                               <  7:21639/148‑1 (MQ=255)
                                                   ggCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCgg                                                                                              >  5:294089/1‑149 (MQ=255)
                                                        cacCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAggg                                                                                         <  7:1250/149‑1 (MQ=255)
                                                          ccccGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGcc                                                                                       >  8:31370/1‑149 (MQ=255)
                                                           cccGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCg                                                                                      >  8:216068/1‑149 (MQ=255)
                                                           cccGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCg                                                                                      >  8:221950/1‑149 (MQ=255)
                                                            ccGAGTGTGGGTACAGGGCGCAACTGGTCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGCGAGTGCTGAGCTGATCGACATCGCCAAGGCTATCACATTCATCCCCTCCGGGAGGGACGt                                                                                     >  3:13650/1‑149 (MQ=255)
                                                             cGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTc                                                                                    >  2:32682/1‑149 (MQ=255)
                                                                     ggTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTccccc                                                                            >  1:350733/1‑149 (MQ=255)
                                                                      gTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTcccccc                                                                           >  6:98753/1‑149 (MQ=255)
                                                                          aGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTca                                                                       >  3:61904/1‑149 (MQ=255)
                                                                           gggCCCACCCGGGCAACGTGCCACCGAATCCGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCAGCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTcac                                                                      <  6:387305/149‑1 (MQ=255)
                                                                            ggCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTACGGGAGGGCCGTCGGTCCCCCCTTcaca                                                                     >  4:402837/1‑149 (MQ=255)
                                                                                 aaCTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACAtttct                                                                >  2:261508/1‑149 (MQ=255)
                                                                                   cTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCCCCGGGAGGGCCGTCGGTCCCCCCTTCACAtttcttt                                                              <  8:146006/149‑1 (MQ=255)
                                                                                    tGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTg                                                                                                             <  7:320862/101‑1 (MQ=255)
                                                                                    tGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTg                                                                                                             >  8:320862/1‑101 (MQ=255)
                                                                                    tGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGCCGTCCTCAGCTGATCTACCTCGCCAAGGCGTTCTCATTCATCACCTCCGGGACGGCCGTCCGTCCCACCATAACAtttcttt                                                              >  7:317771/1‑148 (MQ=255)
                                                                                      ggCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCa                                                            <  5:98753/148‑1 (MQ=255)
                                                                                      ggCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCa                                                            <  4:61904/148‑1 (MQ=255)
                                                                                       gCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGt                                                          <  2:198995/149‑1 (MQ=255)
                                                                                        cAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCATCGGTCCCCCCTTCACATTTCTTTCAGtt                                                         >  8:299076/1‑149 (MQ=255)
                                                                                             tGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTc                                                    <  7:273009/149‑1 (MQ=255)
                                                                                              gCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCg                                                   <  2:350733/149‑1 (MQ=255)
                                                                                                cACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTc                                                 >  6:352205/1‑149 (MQ=255)
                                                                                                 aCCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTc                                                 >  1:353767/1‑148 (MQ=255)
                                                                                                    gATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCGTTCACATTTCTTTCAGTTCATTCGTCTAAc                                             <  1:32682/149‑1 (MQ=255)
                                                                                                     acgatgATGCTGTTGTCATCGGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAAcc                                            >  4:245317/6‑149 (MQ=255)
                                                                                                         cGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTa                                        <  7:12788/149‑1 (MQ=255)
                                                                                                          gATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTAc                                       >  8:18581/1‑149 (MQ=255)
                                                                                                               cGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGcc                                  <  2:222484/149‑1 (MQ=255)
                                                                                                                    acaacaGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATc                             >  4:136702/1‑149 (MQ=255)
                                                                                                                     caacaGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCg                            >  5:250461/1‑149 (MQ=255)
                                                                                                                       tcggCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCt                          <  3:245317/146‑1 (MQ=255)
                                                                                                                             tcatcGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAg                    >  3:193466/1‑149 (MQ=255)
                                                                                                                                  gTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGc               >  8:53566/1‑149 (MQ=255)
                                                                                                                                   tttCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGGCAGCt              >  8:252720/1‑149 (MQ=255)
                                                                                                                                   tttCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCt              >  4:251014/1‑149 (MQ=255)
                                                                                                                                   tttCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCATCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCt              <  7:299076/149‑1 (MQ=255)
                                                                                                                                   gctCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTc                                                    >  1:333963/3‑111 (MQ=255)
                                                                                                                                   gctCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTc                                                    <  2:333963/109‑1 (MQ=255)
                                                                                                                                    ctCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTc             >  4:94911/2‑149 (MQ=255)
                                                                                                                                     tCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAAc                                             <  4:172291/116‑1 (MQ=255)
                                                                                                                                     tCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAAc                                             >  3:172291/1‑116 (MQ=255)
                                                                                                                                         cACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCGATGGTGTTCTGATTCATCCCCTGCGGGAGGGGCGTCGGTCAGGCCTTCTCATTTCTTTGAGTTCATTCGTCTTACCTGTTCTTGGCCTATTTCTGGCAAGCCAGCttcgac         <  4:13650/148‑6 (MQ=255)
                                                                                                                                         cACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCcac         >  5:355100/1‑148 (MQ=255)
                                                                                                                                          aCGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTTGCCCTATCGCTGGCAAGCCAGCTCCaaca        >  3:256403/1‑148 (MQ=255)
                                                                                                                                          aCGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCcacaa       <  1:261508/149‑1 (MQ=255)
                                                                                                                                            ggCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCcacaaca     <  6:210152/149‑1 (MQ=255)
                                                                                                                                               aTCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCcacaacaca   >  8:151539/1‑148 (MQ=255)
                                                                                                                                                tCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTg                                     <  4:200295/113‑1 (MQ=255)
                                                                                                                                                tCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCcacaacacaa  >  1:388468/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCCACAACACAA  >  NC_002947/262799‑263090

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