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.G94.3% 215.9 / 6.6 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 . (4/0);  new base G (31/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01
*NC_0029471,070,2462.A94.4% 219.4 / 6.4 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 . (4/0);  new base A (31/36);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 5.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  NC_002947/1070105‑1070380
                                                                                                                                                  ||                                                                                                                                         
ccgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑G‑‑agatc                                                                                                                                   >  2:21087/1‑144 (MQ=255)
  gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑G‑‑agatcgag                                                                                                                                >  4:254657/1‑142 (MQ=255)
        tCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGcgccgc                                                                                                                            >  2:2713/1‑149 (MQ=255)
           cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCcgcgcg                                                                                                                         >  5:28895/1‑149 (MQ=255)
           cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCcgcgcg                                                                                                                         <  3:254657/149‑1 (MQ=255)
                 ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGTGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCgcgggcgg                                                                                                                    <  4:30191/148‑1 (MQ=255)
                 ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                     <  6:339941/147‑1 (MQ=255)
                 ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCgcggg                                                                                                                       <  3:369820/145‑1 (MQ=255)
                 ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCccgcgcgg                                                                                                                    >  1:191809/1‑148 (MQ=39)
                 ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCccgcgcgg                                                                                                                    <  2:191809/148‑1 (MQ=39)
                   tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGCTGGCTTTCAAGCTGAAGACGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                                                 <  8:297130/149‑1 (MQ=39)
                     aTGCTGATCCCGGCATGGATCGCTGTGCTGTGCGTGGCTTAAAACCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGATAGGTTTTCGGCGCCAGTGA‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                               >  4:367662/1‑149 (MQ=38)
                        cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                            >  8:289918/1‑149 (MQ=255)
                           aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑TGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                          <  5:283875/148‑1 (MQ=255)
                              ccGGCTTGGATCGCTGTGCTGTGGGTGGCTTTCATGCTGAAGAAGGCTCCCCCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                        <  1:344232/149‑1 (MQ=37)
                              ccGGCATGTATCGCTGTGCTGTGGGTGGCTTTCAAGCTTAAGAAGGCT‑‑CTCAAGGCCAAGGGCTGATCGGGTTTTTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                      <  6:283315/149‑1 (MQ=38)
                               cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                <  8:108185/138‑1 (MQ=38)
                               cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                >  7:108185/1‑138 (MQ=38)
                               cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGATAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGAGCGGGCGGCGCTCTATCACACAg                                                                                                     >  4:250178/1‑149 (MQ=38)
                                      gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCCCCGCGCGGGCGGCGCTCGATCTCACAGGCGATg                                                                                              >  7:342972/1‑148 (MQ=255)
                                           ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                >  1:389883/1‑126 (MQ=38)
                                           ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                <  2:389883/126‑1 (MQ=38)
                                                   tgggtggCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGATTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                                >  1:104775/1‑149 (MQ=255)
                                                   cgggtcgCTTTCAAGCTGAAGAAGGCT‑‑CTCCAGGCCAAGGGCTGATCGTGTTTGGTCGTAAGGGGTATGCCGGTAGGTTTTCGGCGCC‑‑TGTTAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                                <  6:354443/148‑1 (MQ=37)
                                                       tggCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                         <  5:272994/136‑1 (MQ=255)
                                                       tggCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                            >  3:383386/1‑149 (MQ=255)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                                     >  7:320755/1‑136 (MQ=255)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                                     <  8:320755/136‑1 (MQ=255)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        >  2:92598/1‑149 (MQ=255)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        >  2:149037/1‑149 (MQ=255)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                              <  2:159061/112‑1 (MQ=37)
                                                           tttCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                              >  1:159061/1‑112 (MQ=37)
                                                             taaaGATGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTAAGGTTTTCGGCGCC‑‑TGT‑‑GACATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                               >  7:263434/3‑109 (MQ=37)
                                                             tCAAGCTGACGAAGCCT‑‑CGCCAGGCCAAGGGCGGATCGGGTTTGTTCGTACGGGGTAGGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                        <  6:256143/147‑1 (MQ=37)
                                                             tCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                               >  7:263433/1‑109 (MQ=255)
                                                             tCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                               <  8:263433/109‑1 (MQ=255)
                                                             tCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATGTCACAGGCCATGAAAGCCTAAAGGCATGCGTCTggc                                                                      >  3:74890/1‑149 (MQ=38)
                                                                     aagaagGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGa                                                              >  3:135999/1‑149 (MQ=255)
                                                                            cT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                   <  3:74867/121‑1 (MQ=255)
                                                                            cT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                   >  4:74867/1‑121 (MQ=255)
                                                                            cT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                       <  1:149037/149‑1 (MQ=255)
                                                                             t‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                      >  4:375694/1‑149 (MQ=255)
                                                                                 gCCAGGCCAAGGTCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGTGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                      <  1:92598/148‑1 (MQ=38)
                                                                                     ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                <  4:383386/149‑1 (MQ=255)
                                                                                        cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                             <  8:342972/149‑1 (MQ=255)
                                                                                        cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                             >  6:207667/1‑149 (MQ=255)
                                                                                         aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                         <  5:211959/104‑1 (MQ=255)
                                                                                         aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                         >  6:211959/1‑104 (MQ=255)
                                                                                                 aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                                      <  1:2713/147‑1 (MQ=255)
                                                                                                  tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                   <  7:289918/149‑1 (MQ=255)
                                                                                                   cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                                  >  8:220176/1‑149 (MQ=255)
                                                                                                    gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                 >  1:49190/1‑149 (MQ=255)
                                                                                                      gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                               >  4:375321/1‑149 (MQ=255)
                                                                                                      gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                               >  7:392782/1‑149 (MQ=255)
                                                                                                      gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAagctgc                               >  1:212616/1‑149 (MQ=255)
                                                                                                           ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                           <  6:28895/148‑1 (MQ=255)
                                                                                                               tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                        >  2:268064/1‑147 (MQ=255)
                                                                                                                  gggtTATGCCGTTAGGGTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAGTGCTGCCTCCGCCAAgca                   <  3:367662/149‑1 (MQ=38)
                                                                                                                   gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                  <  4:135999/149‑1 (MQ=255)
                                                                                                                   gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                  <  2:104775/149‑1 (MQ=255)
                                                                                                                        tGCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgcac              <  3:375694/148‑1 (MQ=255)
                                                                                                                         gCCGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAg            >  1:60103/1‑149 (MQ=255)
                                                                                                                           cGTTAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg          >  1:63788/1‑149 (MQ=255)
                                                                                                                            tttAGGCTTTCATCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg         >  4:104925/2‑149 (MQ=21)
                                                                                                                             ttAGGTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTACCCGTACCAGCCCCGAAATGCTGCCACCGCCAAGCACGCACAgtggg        >  7:138376/1‑148 (MQ=25)
                                                                                                                             ttAGGTTTTCGGCGCC‑‑TGTGAG‑‑AGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt        <  3:250178/149‑1 (MQ=255)
                                                                                                                              tAGGCTTTCATCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt        <  3:104925/148‑1 (MQ=21)
                                                                                                                                 gTTTTCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc    <  2:49190/149‑1 (MQ=255)
                                                                                                                                  ttttCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa   <  7:220176/149‑1 (MQ=255)
                                                                                                                                    ttCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAg  <  2:63788/148‑1 (MQ=255)
                                                                                                                                    ttCGGCGCC‑‑TGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAg  <  5:207667/148‑1 (MQ=255)
                                                                                                                                                  ||                                                                                                                                         
CCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCT‑‑CGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  NC_002947/1070105‑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: