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% 167.7 / 3.6 54intergenic (+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/2);  new base G (23/27);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,070,2462.A94.4% 169.7 / 3.7 54intergenic (+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/2);  new base A (24/27);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070104‑1070385
                                                                                                                                                 ||                                                                                                                                              
gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                                  >  2:313809/1‑143 (MQ=255)
gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatc                                                                                                                                          <  3:126983/149‑5 (MQ=255)
            cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCcgcgcg                                                                                                                              >  3:112402/1‑149 (MQ=255)
                  ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCgcggg                                                                                                                            <  5:148929/145‑1 (MQ=255)
                  ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGggcggc                                                                                                                        >  7:126831/1‑149 (MQ=255)
                   gaaaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGgcgc                                                                                                                      <  6:131356/147‑1 (MQ=39)
                        gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                  >  3:259182/1‑149 (MQ=255)
                        gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                  >  3:155972/1‑149 (MQ=255)
                             tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                              <  4:112402/148‑1 (MQ=255)
                             tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                             <  8:201894/149‑1 (MQ=255)
                             tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                             <  1:277531/149‑1 (MQ=255)
                                cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                          >  7:20545/1‑149 (MQ=255)
                                       gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATg                                                                                                    >  6:299399/1‑148 (MQ=255)
                                         tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                 <  3:274122/149‑1 (MQ=255)
                                           ggtgtgctgtgGGTGGCTGTCAAGCTGAAGAAGGGTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTATGGGGTTTGCCTTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                               <  4:213798/147‑1 (MQ=37)
                                            ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                              >  1:155765/1‑149 (MQ=255)
                                             tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                              >  8:271357/1‑148 (MQ=255)
                                               tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                                            >  5:160482/1‑148 (MQ=255)
                                                 ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                         <  4:259182/149‑1 (MQ=255)
                                                  tgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCTCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                                        >  1:64968/1‑149 (MQ=255)
                                                     gggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                    <  4:155972/149‑1 (MQ=255)
                                                       gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                  <  6:160482/149‑1 (MQ=255)
                                                        tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGc                                                                                                        <  3:191173/127‑1 (MQ=255)
                                                          gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                               >  6:223472/1‑149 (MQ=255)
                                                             ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAATCCTAACTGCATGCGTCTgg                                                                            >  8:326191/1‑149 (MQ=38)
                                                                aaGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATCCCGTTAGGTTTTCGGCGCCTGAGA‑‑GAT‑‑CGAGCGCCGCGCGGCCCGCGCTCGATCTCACAGGCGATGAAAGCCTAAATCCATGCGTCTggcgg                                                                         >  3:153085/1‑149 (MQ=255)
                                                                  gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCCCGCGGGCGGCGCGCGATCTCACAGGCGATTAAAGACGAAAAGCAATGCTGTGGCGGcc                                                                       >  2:172770/1‑149 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCTTGCGTCTGGCGGCCCTGAc                                                                  <  7:271357/149‑1 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                  <  3:39964/149‑1 (MQ=255)
                                                                        gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                                                 >  7:239549/1‑149 (MQ=255)
                                                                           ggCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGt                                                              >  5:296306/1‑149 (MQ=255)
                                                                             cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                            <  2:155765/149‑1 (MQ=255)
                                                                              tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                           >  2:85932/1‑149 (MQ=255)
                                                                              tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                           >  2:72154/1‑149 (MQ=255)
                                                                                    ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                     <  2:64968/149‑1 (MQ=255)
                                                                                       cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                                  <  1:72154/149‑1 (MQ=255)
                                                                                          gggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                               <  8:20545/149‑1 (MQ=255)
                                                                                            gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                             <  5:223472/149‑1 (MQ=255)
                                                                                              tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                                           <  8:126831/149‑1 (MQ=255)
                                                                                                aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                         <  1:172770/149‑1 (MQ=255)
                                                                                                    ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctg                                     >  7:117444/1‑149 (MQ=255)
                                                                                                            cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGc                              <  8:239549/148‑1 (MQ=255)
                                                                                                             gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTTAGAT‑‑AGAGCGCCGCGCGGGCGGCGCTCGATCTCACACGCGATGAAAGCCTAAAGGCATGAGTCTGGCGGCCCTGACGCGTTCCCCATACATGCCCCTCAATTCTGCCACCGcc                             >  4:10743/1‑148 (MQ=37)
                                                                                                               aaGGGGTATGCCGCTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                       <  7:164137/120‑1 (MQ=37)
                                                                                                               aaGGGGTATGCCGCTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                       >  8:164137/1‑120 (MQ=37)
                                                                                                                   ggTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg                      >  5:370498/1‑149 (MQ=255)
                                                                                                                       tGCCGTTAGGTTTTCGGGGCCT‑‑GTGAGAT‑‑CGAGCGCGGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTGCTCCAATCACGcaca                  <  4:153085/149‑1 (MQ=38)
                                                                                                                        gCCGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAg                 >  7:75880/1‑149 (MQ=255)
                                                                                                                         ccGTTATGATTTCGGCGCCT‑‑GT‑‑GATAGCGCGCGCCGCGCGGGCGGCGGTCGATCTATCAGGCGATGAAAGCCTAAATGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgt                <  8:290467/149‑1 (MQ=255)
                                                                                                                          cggtAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCTCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg               <  2:301194/146‑1 (MQ=38)
                                                                                                                          cGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                                       >  4:118376/1‑125 (MQ=255)
                                                                                                                          cGTTAGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                                       <  3:118376/125‑1 (MQ=255)
                                                                                                                              aGGTTTTCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGTCTCCGCCAAGCACGCACAgtggtg            <  1:85932/148‑1 (MQ=255)
                                                                                                                                    tCGGCGCCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCa     >  3:103802/1‑149 (MQ=255)
                                                                                                                                       gcgcCT‑‑GTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAaca  <  6:296306/149‑1 (MQ=255)
                                                                                                                                                 ||                                                                                                                                              
GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070104‑1070385

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: