Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,070,248 +GA intergenic (+48/+94) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G92.9% 221.5 / 9.8 71intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base G (36/29);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,070,2462.A93.0% 224.4 / 10.1 71intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base A (36/30);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  NC_002947/1070098‑1070380
                                                                                                                                                     ||                                                                                                                                         
ccTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCTAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                                             <  5:155890/149‑1 (MQ=255)
ccTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                                             >  6:185743/1‑149 (MQ=255)
 cTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt‑‑g                                                                                                                                          <  6:262072/149‑1 (MQ=255)
  tGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                                             >  1:322622/1‑147 (MQ=255)
     ctccgccgATTCTCATCTCGTTAATGCTGATCCCGGCATGTATCGCTGTGCTGTGGGTTGCTTTCACGCTGAAGACGGCTCGCCAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agat                                                                                                                                    <  6:362874/147‑4 (MQ=255)
     ctccgccgATTCGCATCTCGGTACTGCTGATCCCGGCATTGATCGCTGTGCTGTGGGTGGCTTTCACGCTGAAGACGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agat                                                                                                                                    <  6:368167/147‑4 (MQ=255)
     cgccgcCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agat                                                                                                                                    <  8:57165/149‑4 (MQ=255)
       ccgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agatc                                                                                                                                   >  2:101947/1‑144 (MQ=255)
            gATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAgcgc                                                                                                                               >  3:267728/1‑149 (MQ=255)
               tCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGcgccgc                                                                                                                            >  4:18567/1‑149 (MQ=255)
                cGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcg                                                                                                                           >  8:250824/1‑149 (MQ=255)
                 gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcgc                                                                                                                          >  3:24949/1‑149 (MQ=255)
                 gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcgc                                                                                                                          >  6:103487/1‑149 (MQ=255)
                   aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgg                                                                                                                        >  5:86965/1‑149 (MQ=255)
                        ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                     <  6:86965/147‑1 (MQ=255)
                        ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGggcggc                                                                                                                   <  1:56952/149‑1 (MQ=255)
                            aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                <  5:367303/148‑1 (MQ=255)
                            aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                               <  2:322622/149‑1 (MQ=255)
                               cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                            <  4:24949/148‑1 (MQ=255)
                               cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                           <  5:240080/149‑1 (MQ=38)
                                tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                          >  1:335529/1‑149 (MQ=255)
                                  aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                     <  6:208754/137‑1 (MQ=255)
                                  aTCCAGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                     >  5:208754/1‑137 (MQ=255)
                                      cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                    <  7:250824/149‑1 (MQ=255)
                                         cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                    >  1:104449/1‑146 (MQ=255)
                                         cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                    <  2:104449/146‑1 (MQ=255)
                                          aTGGATCGCAGCGCTGTGGCTCGCGTTCCAGCTGACGAAGGCTCCCCACGCCAAGCGCTGCTCGCCCGTGCCCGTCAGGCGTCTGCCGGAAGGCTTTCGGCGCCTGTCAC‑‑AGCCAGCGCCGCGCGGGCCGCGCTCGATCTCACAGGCGa                                                                                                <  4:211418/149‑1 (MQ=255)
                                            ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                               >  1:238707/1‑148 (MQ=255)
                                              aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGa                                                                                             >  2:228471/1‑148 (MQ=255)
                                               tCGCTGTGCTTTGGGTGGCTTTCAATCTGAAGAAGGCTCGCCAGGAGAAGGGCTGAGCGGGTTTGTTCGTAATGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGATGaaa                                                                                           <  6:87313/149‑1 (MQ=37)
                                               tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                           <  5:103487/149‑1 (MQ=255)
                                                 gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                         <  1:247074/149‑1 (MQ=255)
                                                          tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGTGGGCCGGGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGc                                                                                 <  2:112108/148‑1 (MQ=39)
                                                             gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                             >  2:327147/1‑149 (MQ=255)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                                                 <  3:72543/113‑1 (MQ=255)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                                                 >  4:72543/1‑113 (MQ=255)
                                                               ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                                  <  6:188448/111‑1 (MQ=255)
                                                               ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                                  >  5:188448/1‑111 (MQ=255)
                                                                gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                          >  6:297477/1‑149 (MQ=255)
                                                                  tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        >  2:313074/1‑149 (MQ=255)
                                                                  tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCATAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        >  3:245313/1‑149 (MQ=255)
                                                                   ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                       >  7:276414/1‑149 (MQ=255)
                                                                    tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        <  5:185743/147‑1 (MQ=255)
                                                                           gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                                                 >  4:363021/1‑100 (MQ=37)
                                                                           gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                                                 <  3:363021/100‑1 (MQ=37)
                                                                              gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAGCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCcg                                                            >  5:201269/1‑149 (MQ=38)
                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                    >  2:24787/1‑121 (MQ=255)
                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                    <  1:24787/121‑1 (MQ=255)
                                                                                   cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                       <  1:327147/149‑1 (MQ=255)
                                                                                      gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCTGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                     <  1:101947/149‑1 (MQ=38)
                                                                                       ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                                  <  7:14245/134‑1 (MQ=255)
                                                                                       ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                                  >  8:14245/1‑134 (MQ=255)
                                                                                       ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                                   >  6:325246/1‑149 (MQ=255)
                                                                                           gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                               >  3:320489/1‑149 (MQ=255)
                                                                                           gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                               >  5:351256/1‑149 (MQ=255)
                                                                                           gCCAAGGGCTGATCGGGTTTGTTCCTAAGGGGTATGCCGTTAGGTTTTCAGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                               >  1:88670/1‑149 (MQ=38)
                                                                                           gCCAAGGGCTGATCGGGTTGGGTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GACATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCTGCCCTGACGCGTTCCCCTTc                                                 <  7:328057/147‑1 (MQ=38)
                                                                                               aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                                           <  2:335529/149‑1 (MQ=255)
                                                                                                gggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                          >  6:80437/1‑149 (MQ=255)
                                                                                                 ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTTGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                          >  4:15254/1‑148 (MQ=255)
                                                                                                    tggtCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATAGAGCGCCGCGCGGGCGGCGCTCGATCTAACAGGCGATGAAAGCCAAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCATCCATTccccac                                       >  2:131335/4‑146 (MQ=37)
                                                                                                      aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                    <  4:267728/149‑1 (MQ=255)
                                                                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATGTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTGCCCGTGCATGCCCCTACAtgct                                 >  6:200988/1‑149 (MQ=37)
                                                                                                                ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGGTCCCCTTCCATGCGCCGCAAAGCTGCCTCCg                          >  6:8353/1‑149 (MQ=37)
                                                                                                                   gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                        <  5:325246/148‑1 (MQ=255)
                                                                                                                    tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                      <  5:297477/149‑1 (MQ=255)
                                                                                                                      aGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc                    <  5:80437/149‑1 (MQ=255)
                                                                                                                              gCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAg            >  7:271587/1‑149 (MQ=255)
                                                                                                                                cGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg          <  6:351256/149‑1 (MQ=255)
                                                                                                                                cGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg          >  4:232756/1‑149 (MQ=255)
                                                                                                                                 gTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg         >  1:132248/1‑149 (MQ=255)
                                                                                                                                    aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggtg       <  1:313074/148‑1 (MQ=255)
                                                                                                                                      gTTTTTGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc    <  3:15254/149‑1 (MQ=255)
                                                                                                                                      gTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCCAGCCCGCACAGTGGGGTa     >  4:241865/1‑148 (MQ=38)
                                                                                                                                        tttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAg  >  8:332948/1‑149 (MQ=255)
                                                                                                                                                     ||                                                                                                                                         
CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  NC_002947/1070098‑1070380

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: