Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,246:1 +GA 100% intergenic (+46/+96) 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.G93.4% 181.8 / 7.4 61intergenic (+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/3);  new base G (33/24);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029471,070,2462.A93.4% 185.3 / 7.2 61intergenic (+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/3);  new base A (33/23);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAC  >  NC_002947/1070098‑1070384
                                                                                                                                                          ||                                                                                                                                            
ccTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                <  7:268744/149‑1 (MQ=255)
 cTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                >  3:276840/1‑148 (MQ=255)
   gACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGTTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                >  5:145940/1‑146 (MQ=255)
      gccgccGATTCGCTTCTCGGTAATGCTGATCCCGGCATGGATCTCTGTGCTGTGGGTGGCTTTCAAGCTGATGAAGGGTCGCCAGTCCAAGGGCTGATGGGGTTTGTTCGTAAGGGGTATGCGGTTAGGTTTTCGGCGCG‑‑T‑‑GT‑‑Gagatc                                                                                                                                        <  1:155850/149‑5 (MQ=255)
       ccgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                >  7:257779/1‑142 (MQ=255)
             aTTCGCATCTCGGTAATTCTGATCCCGTCATGGATCGCTGTTCTGTGGGTGGCTTTCAAGCTGATGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGcgcc                                                                                                                                 <  4:296013/149‑1 (MQ=255)
                cGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCcgcg                                                                                                                              >  5:302035/1‑149 (MQ=255)
                 gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCcgcgc                                                                                                                             >  3:314265/1‑149 (MQ=255)
                 gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCcgcgc                                                                                                                             >  3:138649/1‑149 (MQ=255)
                      tcGTTAATGCTGATCGCGTCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCTCCAGGCGAATGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCgcgtgcg                                                                                                                        <  8:80369/149‑1 (MQ=255)
                        ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGggcggc                                                                                                                      >  5:131413/1‑149 (MQ=255)
                           aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGACAGACCCCCGGGCGGCGCt                                                                                                                   >  1:56699/1‑149 (MQ=255)
                              gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                >  7:352580/1‑149 (MQ=255)
                              gCTGATCCCGGCATGGATCGCTGTGCTGTGGCTGGCTTTCAACAAGAAGAAGGCTCGCCAGGCCAACGGCTGAACGGGTTTGTTCGTAAGGGGTATGCCGTTAGCTTTTCGGCGCCAGT‑‑GA‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                >  2:226921/1‑149 (MQ=37)
                               cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                >  2:293619/1‑148 (MQ=255)
                               cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                <  1:293619/148‑1 (MQ=255)
                                tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTGGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                              >  5:336150/1‑149 (MQ=255)
                                tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                  <  5:21450/145‑1 (MQ=38)
                                tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                  >  6:21450/1‑145 (MQ=38)
                                       ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                       <  7:76268/149‑1 (MQ=255)
                                         cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGc                                                                                                      >  4:300076/1‑148 (MQ=255)
                                         cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                                     <  3:300076/149‑1 (MQ=255)
                                          agggAGCGCTGTGCTGGGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGTGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                                    <  8:352580/147‑1 (MQ=38)
                                              aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                            <  4:391486/137‑1 (MQ=255)
                                              aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                            >  3:391486/1‑137 (MQ=255)
                                                     tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                                         <  3:47177/149‑1 (MQ=255)
                                                       ctgggGGGGGCTTTCAAGCTGATGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCGGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCTAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                        <  2:173464/148‑1 (MQ=37)
                                                          ggggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGCGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGc                                                                                     <  8:257779/147‑1 (MQ=255)
                                                             gtggCTTTCAAGCTGAATAAGTCTCGCCAGGCCATTGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                 <  2:56699/149‑1 (MQ=38)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                                >  7:356641/1‑149 (MQ=255)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                                >  5:116790/1‑149 (MQ=255)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCAAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCAGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTAACAGGCGATGAAAGCCTAAAGGCATGCg                                                                                >  6:163831/1‑149 (MQ=38)
                                                                  tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                            >  2:274748/1‑149 (MQ=255)
                                                                      aaGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcg                                                                         <  6:223276/148‑1 (MQ=255)
                                                                          tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGATTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                    >  1:111403/1‑149 (MQ=255)
                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGCTCTCCCCGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCGGAcgcggt                                                            >  7:353581/1‑147 (MQ=38)
                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                            <  6:302035/149‑1 (MQ=255)
                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGTCGGCCCTGACGCGtt                                                            >  7:261379/1‑149 (MQ=38)
                                                                                      gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                        >  4:28134/1‑149 (MQ=255)
                                                                                      gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCCCCGCGAGGGCGGCGCTCGATCTCACACGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                        >  4:257212/1‑149 (MQ=39)
                                                                                      gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GAGCTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                        <  8:334577/149‑1 (MQ=255)
                                                                                        cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                      <  1:274748/149‑1 (MQ=255)
                                                                                           gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATATCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                                   >  7:333642/1‑149 (MQ=255)
                                                                                            ccAAGTGCTGATCGGGTTTGTTCGTAAGGGGTATTCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGGCTGGCGGCGCTGACGCGTTCCCCTTCCAt                                                  <  6:131413/149‑1 (MQ=38)
                                                                                             cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                                 >  1:302981/1‑149 (MQ=255)
                                                                                                gggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                              <  4:276840/149‑1 (MQ=255)
                                                                                                      aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGGCTGGCGGCCCTGACGCGGTCCCCGTCCATGCCCCTCCAt                                        >  4:45837/1‑149 (MQ=38)
                                                                                                       tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                       >  7:132448/1‑149 (MQ=255)
                                                                                                       tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                       <  3:28134/149‑1 (MQ=255)
                                                                                                           gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   >  2:319302/1‑149 (MQ=255)
                                                                                                            tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   >  2:266503/1‑148 (MQ=255)
                                                                                                               gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                               >  4:260844/1‑149 (MQ=255)
                                                                                                               gTTCGTAAGGCGTATGCCGTTAGGTTTTCGGCGCC‑‑TGGGA‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                                <  2:305565/148‑1 (MQ=38)
                                                                                                                ttCGTAAGGGTTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                               <  6:145940/148‑1 (MQ=255)
                                                                                                                ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCTC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                               >  1:251137/1‑149 (MQ=38)
                                                                                                                   gTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATCATGCCTCcccc                            >  6:157435/1‑148 (MQ=38)
                                                                                                                       ggggTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATGTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgca                       >  3:200004/1‑149 (MQ=255)
                                                                                                                        gggTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                      <  8:132448/149‑1 (MQ=255)
                                                                                                                                cGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg              <  8:356641/149‑1 (MQ=255)
                                                                                                                                 tttAGGCTTTCATCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCCCGCGCGGCGCTCGCTCTCACAGGCGCTGCAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTGCCATGCCCCGCAATGCTGCCTCCGCCAAGCACGCACAgtgg             >  5:386267/2‑149 (MQ=21)
                                                                                                                                      gTTTTCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGAGCCGCGAGGGATGCGCTCGATCTCACATGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCAACTTACATGCACCTCAATGCTGCCACCGCCAAGCACGCACAGTGGTGTa         >  4:107516/1‑148 (MQ=255)
                                                                                                                                          tCGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCTCCGCGCGGGGGGCGCTCGATCTCAGAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCAGGCACAGTGGTGTACAGc     <  8:261379/148‑1 (MQ=255)
                                                                                                                                           cGGCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa   >  4:90842/1‑149 (MQ=255)
                                                                                                                                            ggCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAac  <  4:314265/149‑1 (MQ=255)
                                                                                                                                            ggCGCC‑‑T‑‑GTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAac  <  4:138649/149‑1 (MQ=255)
                                                                                                                                                          ||                                                                                                                                            
CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAC  >  NC_002947/1070098‑1070384

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: