Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,997,189:1 +G 7.0% intergenic (‑158/+42) PP_5252 ← / ← PP_5253 hypothetical protein/arylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,997,1891.G 7.0% 277.8 / 14.8 71intergenic (‑158/+42)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base . (33/33);  new base G (3/2);  total (36/35)
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.73e-01

ACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT  >  NC_002947/5997043‑5997335
                                                                                                                                                   |                                                                                                                                                  
acgaGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                  <  7:289001/146‑1 (MQ=255)
    ggggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGTCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAATCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                              <  5:288234/149‑1 (MQ=255)
    ggggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                              >  1:18473/1‑149 (MQ=255)
    ggggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                              <  3:66517/149‑1 (MQ=255)
     gggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTg                                                                                                                                             <  6:229174/149‑1 (MQ=255)
     gggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTg                                                                                                                                             <  8:264301/149‑1 (MQ=255)
      ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTg                                                                                                                                             >  4:118028/1‑148 (MQ=255)
      ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTAAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCACCAACAACCTCAGt                                                                                                                                            >  4:22464/1‑149 (MQ=255)
        gcGGCAGGGGGACTGCGATGGATGCCGGTTGATGCAAGTCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGTCCTCTTCGCGTGTAAACCCTCTCCTACAACGTCTGTAg                                                                                                                                          <  7:95232/148‑1 (MQ=255)
               ggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                   <  2:94911/149‑1 (MQ=255)
                  gACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGttt                                                                                                                                <  2:342650/149‑1 (MQ=255)
                  gACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGttt                                                                                                                                <  3:118028/149‑1 (MQ=255)
                     tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAcc                                                                                                                             >  2:130094/1‑149 (MQ=255)
                      gCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGa                                                                                                                                        >  7:312274/1‑137 (MQ=255)
                      gCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGa                                                                                                                                        <  8:312274/137‑1 (MQ=255)
                         aaGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCcgcg                                                                                                                         <  2:18473/149‑1 (MQ=255)
                                  ggTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTcc                                                                                                                <  2:34852/149‑1 (MQ=255)
                                   gTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCCCTACAACAACGTCTGTAGGAGCGGGTTTAACCGCCAATGtttac                                                                                                                >  2:116168/1‑146 (MQ=255)
                                          cAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTa                                                                                                        >  8:189662/1‑149 (MQ=255)
                                           aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCCCACCAAAAACGGCTGTTTGAGCGGGTTTTAACCCGAAAGTTTTCCACAAATa                                                                                                        >  5:63828/1‑148 (MQ=255)
                                            aGGCCTGGGCTTTAATGGATGTGGGGTGGGTTTAAGACGGTTAGCCGTGGTAAAGCTTTTTGCTGCCTGGGCCGGCCTATTCGCGGGTAACCCCGCTCCTACAACGTCTTTTGGTGGGGTTTTACCCGCGAATGTTTCCCACAATTAcc                                                                                                      <  6:63828/149‑1 (MQ=255)
                                            aGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTa                                                                                                        <  7:219671/147‑1 (MQ=255)
                                            aGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTa                                                                                                        >  8:219671/1‑147 (MQ=255)
                                              gCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGGTTTCCCGCGAATGTTTCCCACAATTACCGc                                                                                                    >  8:1182/1‑149 (MQ=255)
                                              gCCGGGGCTTGAATGGATGGGCGGTGGGTTGAAGACGGTACGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGc                                                                                                    <  2:265927/149‑1 (MQ=37)
                                                   ggCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGaa                                                                                                <  1:134424/148‑1 (MQ=255)
                                                   ggCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCAGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAg                                                                                               >  8:212164/1‑149 (MQ=255)
                                                     cTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                       >  8:188725/1‑123 (MQ=255)
                                                     cTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                       <  7:188725/123‑1 (MQ=255)
                                                           tGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTAAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCa                                                                                       <  3:22464/149‑1 (MQ=255)
                                                            ggATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCGCACCAAAAACGTGTGTATGGGCGGGGTTTCACGGGAAAGTTTGACACAAATTACGCTGAAGGAAACCCAg                                                                                      >  7:335088/1‑149 (MQ=255)
                                                               tgtgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                   <  2:33030/149‑1 (MQ=255)
                                                                   gggTGGGTTGATTAAGGTTTCCCCTTGGCATGATTTTTGCTTCGTGGGCCGTCCTGTTCGCGGGTAACCCCGCTACTAAAACGTCTGTTGGAGGGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTcc                                                                               <  8:148835/148‑1 (MQ=255)
                                                                     gTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCTCCTAGAACCTCAGTATGAGAGGGTTTACCCGCGCAATTTTCCCACAATTTACCCGGAAGGTACGCCAGCAAGtaagc                                                                             >  2:36282/1‑145 (MQ=255)
                                                                          ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCTGAAGGAACGCCAGCAAGTCCGCAtt                                                                          >  3:358142/1‑147 (MQ=255)
                                                                          ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAACGTATGTAGGAGCGGGTTTAACCGCGAATGTTTACCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTc                                                                         >  8:16533/1‑148 (MQ=255)
                                                                             aagaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCTGAAGGaa                                                                                            >  5:198931/1‑126 (MQ=255)
                                                                             aagaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCTGAAGGaa                                                                                            <  6:198931/126‑1 (MQ=255)
                                                                               gaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgc                                                                   >  4:311038/1‑149 (MQ=255)
                                                                                     tAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGACCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACCCCAGCAAGTCCTCATTCAGCTGCTGCGCa                                                             >  6:313900/1‑149 (MQ=255)
                                                                                         ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTACTACAACGTCTGTATGAGCCGGTTTAACCGCGAATGTTTCCCACAAATACCGCAGAAGGAACGCCAGCAAGTCCGCAATAACCAGCAGCGCAt                                                            >  6:131935/1‑146 (MQ=37)
                                                                                         ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGATACTACAACCTCACTTAGAACCCGTGTAACCCCGAATGTTTCCCACAAATAAACCAGAACGAACGCCAGCAAATACACAATCAGCAGCTGCGCaata                                                          >  6:224878/1‑145 (MQ=255)
                                                                                              tCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt                                                    <  3:304103/149‑1 (MQ=255)
                                                                                              tCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt                                                    <  3:311038/149‑1 (MQ=255)
                                                                                                      tttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaaaa                                              <  7:183294/147‑1 (MQ=255)
                                                                                                         gctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTGTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGTACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTg                                         >  7:73447/1‑149 (MQ=255)
                                                                                                          ctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAc                                                                                                                              >  6:392193/1‑63 (MQ=255)
                                                                                                          ctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAc                                                                                                                              <  5:392193/63‑1 (MQ=255)
                                                                                                                   ccggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCg                               <  5:313900/149‑1 (MQ=255)
                                                                                                                    cggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCgg                              <  7:189662/149‑1 (MQ=255)
                                                                                                                     ggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAATACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc                             <  1:130094/149‑1 (MQ=255)
                                                                                                                      gccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGTACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCa                            <  8:73447/149‑1 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACAACGTCTGTATTAGGGGGTTTACCCGCGAATGTTTCCCAAAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTTCAAAACCGTGCGGCGCCCCGGCATACACCtt                    <  7:15484/149‑1 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCt                     >  5:71342/1‑148 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCt                     <  6:71342/148‑1 (MQ=255)
                                                                                                                                 ggTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg                 >  2:212458/1‑149 (MQ=255)
                                                                                                                                    aaaCCCGAGCCTACACCGTCTTTCGCAGGGGTTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTc              <  5:67267/149‑1 (MQ=255)
                                                                                                                                      aCCCGATCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGACACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCgg            >  7:173964/1‑149 (MQ=255)
                                                                                                                                       cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCGACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAATCCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc           >  8:310601/1‑149 (MQ=38)
                                                                                                                                       cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTc                   <  1:209445/141‑1 (MQ=255)
                                                                                                                                       cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTc                   >  2:209445/1‑141 (MQ=255)
                                                                                                                                        ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          >  7:66961/1‑149 (MQ=255)
                                                                                                                                        ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          >  2:215686/1‑149 (MQ=255)
                                                                                                                                        ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          >  3:244910/1‑149 (MQ=255)
                                                                                                                                          gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGcc        >  7:226610/1‑149 (MQ=255)
                                                                                                                                             ccTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                          <  8:129775/96‑1 (MQ=255)
                                                                                                                                             ccTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                          >  7:129775/1‑96 (MQ=255)
                                                                                                                                              cTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCACTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACg     >  6:416118/1‑148 (MQ=255)
                                                                                                                                                aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAACACCGTGCGGCGCCCCGGCATAAACCTTCAGCTCGGCGCCACGAAt  >  4:261163/1‑149 (MQ=255)
                                                                                                                                                 cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt  <  6:136933/148‑1 (MQ=255)
                                                                                                                                                 cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt  >  5:136933/1‑148 (MQ=255)
                                                                                                                                                  aaCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc           <  5:156797/138‑1 (MQ=255)
                                                                                                                                                  aaCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCAACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc           >  6:156797/1‑138 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT  >  NC_002947/5997043‑5997335

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