Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,997,189:1 +G 11.5% 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.G11.5% 204.9 / 25.3 61intergenic (‑158/+42)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base . (22/32);  new base G (4/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

ACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT  >  NC_002947/5997043‑5997335
                                                                                                                                                     |                                                                                                                                                   
aCTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                   <  3:211608/149‑1 (MQ=255)
aCTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                   <  2:56162/149‑1 (MQ=255)
       gCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                            <  4:358988/149‑1 (MQ=255)
         cGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAgg                                                                                                                                          >  4:84263/1‑149 (MQ=255)
          ggCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGAAGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCAGCTACAACGTCTGTAg                                                                                                                                           >  1:111292/1‑149 (MQ=39)
          ggCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAg                                                                                                                                           >  2:62360/1‑147 (MQ=255)
               ggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                    <  8:316133/148‑1 (MQ=255)
                gggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                            <  1:199358/140‑1 (MQ=255)
                gggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                            >  2:199358/1‑140 (MQ=255)
                gggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGt                                                                                                                                  <  3:272954/149‑1 (MQ=255)
                gggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGt                                                                                                                                  >  8:22747/1‑149 (MQ=255)
                 ggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGtt                                                                                                                                 <  7:333945/149‑1 (MQ=255)
                    cTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                             >  8:237413/1‑135 (MQ=255)
                    cTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                             <  7:237413/135‑1 (MQ=255)
                       cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                    >  3:39459/1‑140 (MQ=255)
                       cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                    <  4:39459/140‑1 (MQ=255)
                       cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCACCTACAACGTCTGTAGGAGCGGGTTTACCcg                                                                                                                           >  3:364488/1‑149 (MQ=255)
                             aTGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                      >  2:264857/1‑132 (MQ=255)
                             aTGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                      <  1:264857/132‑1 (MQ=255)
                                 cGGTTTATGCATGGCCCGGGCTGTAACTGAGGTGCGGGGG‑‑TTTGATGAAGGTTCGCCCGGGTCAAGCTTTTTGCTGCGTGGGCGGGCGTCTTGGCTGGTAAACCCGCTCCTACAACGTCTTTAGGAGGGGGTTTACCCGCGAATGTTTc                                                                                                                 <  8:128232/149‑1 (MQ=255)
                                        tGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGCTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAt                                                                                                          <  7:19706/149‑1 (MQ=255)
                                             ggCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCg                                                                                                     <  3:84263/149‑1 (MQ=255)
                                               ccTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAc                                                                                                       >  3:158866/1‑145 (MQ=255)
                                               ccTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAc                                                                                                       <  4:158866/145‑1 (MQ=255)
                                                cTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCt                                                                                                   >  7:115354/1‑148 (MQ=255)
                                                   ggCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGaa                                                                                                <  2:111292/148‑1 (MQ=255)
                                                       tGAATGCCTGGGCGGGGG‑‑GTGGAAGACGGTAAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAAc                                                                                           <  6:279855/149‑1 (MQ=37)
                                                         aaTGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGaa                                                                                                <  7:319927/142‑1 (MQ=255)
                                                             gATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGGGCCTACAAGGGGGGTGGGGGCGGGGTTAGCCGCGAATGTTTCCCCCAATTACCGCGGAAGGAACGCCAGc                                                                                     >  2:4447/1‑149 (MQ=37)
                                                                     gTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                      <  4:51411/92‑1 (MQ=37)
                                                                     gTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                      >  3:51411/1‑92 (MQ=37)
                                                                      tGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCa                                                                            <  8:115354/149‑1 (MQ=255)
                                                                           gTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTc                                                                         <  7:22747/149‑1 (MQ=255)
                                                                            ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCAGCAAAAAAGTAGGTAGGAGCCGGGTTTACCGGCTAATTTTACCACAAAAAACGCTGAAGTAACCCCACCAAATccacctta                                                                         >  6:322320/1‑142 (MQ=255)
                                                                              gaagaaGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTGCGACAACGTGTGTTGGAGCGGGGTTTCCCGCGAATGTTTGCCACAATTACCGCTGGAGGAACGCCAGCAAGTCCTCATTAAgc                                                                      >  7:328257/1‑149 (MQ=38)
                                                                               aagaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTGCTACAACGTGTGTAGGAGAGGGGTTACAAGGGAATGTTTACCACAATTTAAGCTGAAGGAACGCCAGCAAGTCCTCATTAAgct                                                                     >  7:319815/1‑149 (MQ=37)
                                                                                      ttAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTgcgc                                                              <  1:331632/149‑1 (MQ=255)
                                                                                        aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCa                                                             <  4:364488/148‑1 (MQ=255)
                                                                                          cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                          <  3:229383/149‑1 (MQ=255)
                                                                                          cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                          <  1:62360/149‑1 (MQ=255)
                                                                                          cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGTCTGTAGGAGCGGGTTTACCCGCGATTGTTTCCGAGAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                          <  5:256576/149‑1 (MQ=255)
                                                                                           ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAg                                                         >  5:294258/1‑149 (MQ=255)
                                                                                               gTCAAGCTTTTTGCTGCCTGGGCGGGCCTCTTCGCGGGTAAACCCGATCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCa                                                      <  5:322320/148‑1 (MQ=255)
                                                                                                   aGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCGCTCCTACAACGTGTGTAGGAGCGGGTTTACCCGCGCATGTTTGCCACAATTTCCGCTGAAGGAACGCCCCCAAGTCCAAATTACCCAGCGGCGCATTAGTAACTgg                                                  >  6:122904/1‑147 (MQ=21)
                                                                                                    gCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc                                                                                                                 <  7:6422/84‑1 (MQ=255)
                                                                                                           gctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGt                                          <  7:128187/148‑1 (MQ=255)
                                                                                                           gctgcCTGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCGAATTCAGATGCTGCGCATGAGTCACTGCAAAACCGTgc                                        >  4:300040/1‑149 (MQ=38)
                                                                                                                       ggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc                             <  3:205135/149‑1 (MQ=255)
                                                                                                                       ggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc                             <  8:328257/149‑1 (MQ=255)
                                                                                                                          ctctTCGCGGGTAAACCCCCTCCTACAACGTCTGTAGGAGCGGGTTTTCCCGCGAATGTTTCACACAATTAACGCTGAAGGAAAGCCCGCCAGTCCTCATTCAGCTGCGGCGGATGAGTCAGTGGAAAACCGTGCGGCGCCCCGGGAt                           >  4:296334/1‑148 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg                 <  3:17031/144‑1 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg                 >  4:17031/1‑144 (MQ=255)
                                                                                                                                         cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc           >  8:63550/1‑149 (MQ=255)
                                                                                                                                         cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc           >  7:218499/1‑149 (MQ=255)
                                                                                                                                          ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          >  6:346327/1‑149 (MQ=255)
                                                                                                                                          ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          <  6:294258/149‑1 (MQ=255)
                                                                                                                                          ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg          >  4:318355/1‑149 (MQ=255)
                                                                                                                                                cTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCAc      <  7:63550/147‑1 (MQ=255)
                                                                                                                                                cTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGa    <  5:122904/149‑1 (MQ=255)
                                                                                                                                                 tACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGaa   >  8:266867/1‑149 (MQ=255)
                                                                                                                                                   cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt  <  1:4447/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGG‑‑GTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT  >  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: