Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 262,940 G→C 6.9% intergenic (+175/‑310) 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,9400GC 6.9% 223.7 / 6.4 73intergenic (+175/‑310)PP_0211/gabD‑Ihypothetical protein/succinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (41/27);  new base C (3/2);  total (44/29)
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.36e-01

AGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCCACAACA  >  NC_002947/262797‑263087
                                                                                                                                               |                                                                                                                                                    
aGAAGCAGGAAGACCTCCTAGCCAACGATCTGGCAGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGAACGGCAt                                                                                                                                                   >  7:184601/1‑147 (MQ=255)
  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                                                                                                                                       <  2:66032/135‑1 (MQ=255)
                       aaCGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgct                                                                                                                                       >  1:66032/1‑135 (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                                                                                                                                >  5:121160/1‑122 (MQ=255)
                                           ccTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTgctgagctga                                                                                                                                <  6:121160/122‑1 (MQ=255)
                                            cTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCAt                                                                                                    <  8:273231/149‑1 (MQ=255)
                                                  tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATccc                                                                                                 <  4:326161/146‑1 (MQ=255)
                                                  tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                              >  6:390227/1‑149 (MQ=255)
                                                  tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                              >  5:210152/1‑149 (MQ=255)
                                                  tttGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTc                                                                                              >  8:12788/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:221950/1‑149 (MQ=255)
                                                             cccGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCg                                                                                   >  8:216068/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                                                 <  2:333963/109‑1 (MQ=255)
                                                                                                                                     gctCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTc                                                 >  1:333963/3‑111 (MQ=255)
                                                                                                                                      ctCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTc          >  4:94911/2‑149 (MQ=255)
                                                                                                                                       tCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAAc                                          >  3:172291/1‑116 (MQ=255)
                                                                                                                                       tCAGCACGCCCTGATGCCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAAc                                          <  4:172291/116‑1 (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)
                                                                                                                                               |                                                                                                                                                    
AGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTACGATGCCGATGACAACAGCATCATCGTTTCAGCACGGCATCAGGGCGTGCTGAGCTGATCGACCTCGCCAAGGCGTTCTGATTCATCCCCTCCGGGAGGGCCGTCGGTCCCCCCTTCACATTTCTTTCAGTTCATTCGTCTAACCTGTACTGGCCCTATCGCTGGCAAGCCAGCTCCCACAACA  >  NC_002947/262797‑263087

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: