Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,246,185 C→G 8.1% M30I (ATG→ATC)  alr ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,246,1850CG 8.1% 202.6 / 4.5 69M30I (ATG→ATC) alralanine racemase
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base G (3/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAG  >  NC_002947/4246038‑4246332
                                                                                                                                                   |                                                                                                                                                   
aTAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCa                                                                                                                                                    >  7:101585/1‑149 (MQ=255)
    gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                  <  4:25356/147‑1 (MQ=255)
    gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                  >  3:25356/1‑147 (MQ=255)
         ggCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACg                                                                                                                                           >  5:167309/1‑149 (MQ=255)
         ggCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACg                                                                                                                                           <  8:101585/149‑1 (MQ=255)
            cTTGAGCACGGCGCACAGCTTGGACTTTCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtg                                                                                                                                        <  6:199441/149‑1 (MQ=255)
                            aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAgggg                                                                                                                        <  4:323552/149‑1 (MQ=255)
                            aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAgggg                                                                                                                        >  2:34493/1‑149 (MQ=255)
                               ttGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                     >  1:102469/1‑149 (MQ=255)
                                     ttGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                         >  6:348110/1‑139 (MQ=255)
                                     ttGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                         <  5:348110/139‑1 (MQ=255)
                                       gCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCg                                                                                                              <  4:224709/148‑1 (MQ=255)
                                         cGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTg                                                                                                           <  5:368338/149‑1 (MQ=255)
                                               gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAg                                                                                                      >  8:212395/1‑148 (MQ=255)
                                               gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAg                                                                                                      >  8:218158/1‑148 (MQ=255)
                                                    gCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGc                                                                                                <  7:212395/149‑1 (MQ=255)
                                                    gCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGc                                                                                                <  7:218158/149‑1 (MQ=255)
                                                     ccTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTGGGTGCCGTTGTCCATCGACAACTGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGc                                                                                                >  5:35990/1‑148 (MQ=255)
                                                      cTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGc                                                                                                >  8:24134/1‑147 (MQ=255)
                                                      cTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGc                                                                                                >  4:349549/1‑147 (MQ=255)
                                                       tGCAGCGTGCGGATGTTGTGCTGCAGGGCACTGGCGCTGACTGCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGGGTTGGTGCCGTTGTCGATGGACAACggcac                                                                                                                                       <  8:262471/107‑4 (MQ=255)
                                                         cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAg                                                                                             <  5:284145/147‑1 (MQ=255)
                                                         cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAg                                                                                             <  3:233584/147‑1 (MQ=255)
                                                           gCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGAt                                                                                         <  6:322754/149‑1 (MQ=255)
                                                                cGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTCTCCGGTGATCAGAAGTGCCAGGGATGCAGc                                                                                    <  1:34493/149‑1 (MQ=255)
                                                                   aTGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCTACAACGGTGGTGCAGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCAAGGGATGCAGc                                                                                    >  8:373454/1‑146 (MQ=255)
                                                                      ttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAg                                                                              >  3:404650/1‑149 (MQ=255)
                                                                      ttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAg                                                                              <  3:354832/149‑1 (MQ=255)
                                                                      ttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAg                                                                              <  3:100181/149‑1 (MQ=255)
                                                                        gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                            >  8:244523/1‑149 (MQ=255)
                                                                         tgctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                            <  7:24134/148‑1 (MQ=255)
                                                                            tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                               <  4:143918/142‑1 (MQ=255)
                                                                            tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                               >  3:143918/1‑142 (MQ=255)
                                                                            tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACg                                                                        >  5:184156/1‑149 (MQ=255)
                                                                              cAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGc                                                                      <  1:265976/149‑1 (MQ=255)
                                                                                  gcgctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGtaa                                                                  <  6:167309/149‑1 (MQ=255)
                                                                                  gcgctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGaaa                                                                  <  1:283971/149‑1 (MQ=255)
                                                                                      tggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCa                                                                                                       >  2:216481/1‑108 (MQ=255)
                                                                                      tggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCa                                                                                                       <  1:216481/108‑1 (MQ=255)
                                                                                                ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTa                                                    >  2:217281/1‑149 (MQ=255)
                                                                                                   gACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCTACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCTGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                                 <  7:373454/149‑1 (MQ=255)
                                                                                                     cccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgt                                               >  8:51417/1‑149 (MQ=255)
                                                                                                        aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgtgt                                             >  4:146551/1‑148 (MQ=255)
                                                                                                          gCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCt                                          <  2:102469/149‑1 (MQ=255)
                                                                                                           cATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCtt                                         <  7:51417/149‑1 (MQ=255)
                                                                                                            aTTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTg                                        >  2:375850/1‑149 (MQ=255)
                                                                                                               gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCa                                                                                                       <  8:72323/83‑1 (MQ=255)
                                                                                                               gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCa                                                                                                       >  7:72323/1‑83 (MQ=255)
                                                                                                                 tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                   <  6:184156/149‑1 (MQ=255)
                                                                                                                  gTTTTGCACGGTCAGGGTGTTGGTGACGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                  >  6:109190/1‑149 (MQ=255)
                                                                                                                         aCGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTgcgcgc                           >  4:250904/1‑149 (MQ=255)
                                                                                                                                 ggTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                      >  4:193113/1‑146 (MQ=255)
                                                                                                                                  gtgtTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGa                  >  8:222145/1‑149 (MQ=255)
                                                                                                                                     ttGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACg                >  1:381143/1‑148 (MQ=255)
                                                                                                                                     ttGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACg                >  1:387190/1‑148 (MQ=255)
                                                                                                                                     ttGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCAAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACg                >  5:176313/1‑148 (MQ=255)
                                                                                                                                      tGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                  >  1:47098/1‑129 (MQ=255)
                                                                                                                                      tGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                  <  2:47098/129‑1 (MQ=255)
                                                                                                                                        gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAgg            >  2:125326/1‑149 (MQ=255)
                                                                                                                                        gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg             <  4:404650/148‑1 (MQ=255)
                                                                                                                                        gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg             <  3:146551/148‑1 (MQ=255)
                                                                                                                                         tGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt           <  7:244523/149‑1 (MQ=255)
                                                                                                                                         caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCtt                        >  3:43617/3‑136 (MQ=255)
                                                                                                                                         caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCtt                        <  4:43617/134‑1 (MQ=255)
                                                                                                                                         caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt           >  5:116220/3‑149 (MQ=255)
                                                                                                                                          gCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg          >  2:44315/1‑149 (MQ=255)
                                                                                                                                          accGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTT‑ACAGACGCAGgtgt         >  7:147971/2‑149 (MQ=255)
                                                                                                                                              ttGTCCATCGACAAGGTTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGtt      <  5:109190/149‑1 (MQ=255)
                                                                                                                                                gTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGttt     <  3:250904/148‑1 (MQ=255)
                                                                                                                                                  ccATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAg  >  7:74510/1‑149 (MQ=255)
                                                                                                                                                   cATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAg  <  7:222145/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAG  >  NC_002947/4246038‑4246332

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