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.4% 474.4 / 39.1 186intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (6/8);  new base G (98/73);  total (104/82)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01
*NC_0029471,070,2462.A92.5% 450.4 / 37.0 186intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (6/8);  new base A (98/74);  total (104/82)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 192 total aligned reads displayed.

GGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGC  >  NC_002947/1070122‑1070362
                                                                                                                             ||                                                                                                                      
ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                          <  1:859204/125‑1 (MQ=255)
ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                          >  2:859463/1‑125 (MQ=255)
  tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑a                                                                                                                    <  1:859205/125‑1 (MQ=255)
  tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑a                                                                                                                    >  2:859688/1‑125 (MQ=255)
   aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ag                                                                                                                   <  2:859207/125‑2 (MQ=255)
    aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑aga                                                                                                                  <  1:859208/125‑3 (MQ=255)
      gCTGATCCCGGCATGGATCGCTGTGCTGTGTGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agatc                                                                                                                <  1:859211/125‑5 (MQ=255)
      gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agatc                                                                                                                >  2:859505/1‑121 (MQ=255)
       cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑agatcg                                                                                                               >  2:859425/1‑120 (MQ=255)
           tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTGTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑aga                                                                                                                  >  2:859526/1‑116 (MQ=255)
                gCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcg                                                                                                        <  1:859217/125‑1 (MQ=255)
                gCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcg                                                                                                        >  2:859540/1‑125 (MQ=255)
                 cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcg                                                                                                        >  2:859449/1‑124 (MQ=255)
                    ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcggg                                                                                                    >  1:859430/1‑125 (MQ=255)
                    ggATCGCTGTGCTGTGGGTGGCTGTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcggg                                                                                                    <  2:859218/125‑1 (MQ=39)
                      aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                  >  1:859421/1‑125 (MQ=255)
                      aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggc                                                                                                   >  2:859529/1‑124 (MQ=255)
                          ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGgcgc                                                                                              >  1:859569/1‑125 (MQ=255)
                          ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTATTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                  >  1:859524/1‑121 (MQ=39)
                             tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCt                                                                                             <  2:859220/123‑1 (MQ=255)
                                tgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                         <  1:859223/124‑1 (MQ=255)
                                    ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                       <  2:859224/122‑1 (MQ=255)
                                     gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTgaca                                                                                   >  1:859577/1‑125 (MQ=38)
                                      tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                  >  1:859480/1‑125 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgggg                                                                               >  2:859468/1‑122 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                              <  1:859227/125‑1 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGAGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                 <  2:859228/122‑1 (MQ=38)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                         <  2:859232/125‑1 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                         <  2:859235/125‑1 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                         <  1:859231/125‑1 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                         <  1:859234/125‑1 (MQ=255)
                                                  tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                      >  2:859565/1‑125 (MQ=255)
                                                   gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                     <  2:859238/125‑1 (MQ=255)
                                                    aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                    >  1:859386/1‑125 (MQ=255)
                                                    aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTGGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCGGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                    >  2:859358/1‑125 (MQ=37)
                                                     agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                   <  1:859240/125‑1 (MQ=255)
                                                      gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                      >  2:859584/1‑121 (MQ=255)
                                                      gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaga                                                                  >  2:859337/1‑123 (MQ=255)
                                                         ggCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                  <  2:859242/122‑1 (MQ=255)
                                                          gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                              <  2:859244/125‑1 (MQ=255)
                                                           cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                             <  1:859246/125‑1 (MQ=255)
                                                           cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCAGCGCGGCCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                             <  1:859248/125‑1 (MQ=38)
                                                            tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                            >  1:859536/1‑125 (MQ=255)
                                                            tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                            >  2:859346/1‑125 (MQ=255)
                                                                cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                        <  2:859249/125‑1 (MQ=255)
                                                                cAGGCCAAGGGCTGAGCGGGTTTGTTCGTAAGGGGTATGGCGTTAGGTTTTCGGCGCTTGT‑‑GATATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                        <  2:859251/125‑1 (MQ=37)
                                                                 aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                       <  1:859253/125‑1 (MQ=255)
                                                                  ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                       <  1:859254/124‑1 (MQ=255)
                                                                  ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                      >  1:859451/1‑125 (MQ=255)
                                                                     cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcg                                                   >  1:859414/1‑125 (MQ=255)
                                                                     cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcg                                                   <  2:859256/125‑1 (MQ=255)
                                                                       aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcgg                                                  >  2:859423/1‑124 (MQ=255)
                                                                           cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                             >  2:859655/1‑125 (MQ=255)
                                                                            tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                             >  2:859556/1‑124 (MQ=255)
                                                                             gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                               >  1:859516/1‑121 (MQ=255)
                                                                             gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                           >  1:859341/1‑125 (MQ=255)
                                                                              aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                          >  2:859390/1‑125 (MQ=255)
                                                                              aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                          >  1:859440/1‑125 (MQ=255)
                                                                                 gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGggcggc                                                                                                >  2:859537/1‑68 (MQ=21)
                                                                                 gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGt                                       <  2:859259/125‑1 (MQ=255)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                        <  2:859261/91‑1 (MQ=25)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                 >  2:859622/1‑114 (MQ=255)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                      >  2:859512/1‑125 (MQ=255)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                      >  1:859493/1‑125 (MQ=255)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGt                                       >  1:859474/1‑124 (MQ=255)
                                                                                    tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                    >  2:859508/1‑125 (MQ=255)
                                                                                     ttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                   >  2:859416/1‑125 (MQ=255)
                                                                                       gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                 <  1:859263/125‑1 (MQ=255)
                                                                                         tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCTGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTAGCGGCCCTGACGCGTTCCCCTTc                               >  2:859702/1‑125 (MQ=37)
                                                                                           gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                             <  1:859265/125‑1 (MQ=255)
                                                                                           gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGAGGCCCTGACGCGTTCCCCTTCCa                             <  2:859266/125‑1 (MQ=37)
                                                                                            tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                            >  1:859499/1‑125 (MQ=255)
                                                                                             aaGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                           >  2:859595/1‑109 (MQ=255)
                                                                                             aaGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                           <  1:859268/125‑1 (MQ=255)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGCGCGCCGCGCGCGCGGCGCTCGATCTCATAGGCGATGAAAGCCTACAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                            >  2:859614/1‑122 (MQ=25)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGTCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                         <  1:859270/125‑1 (MQ=37)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                         >  1:859623/1‑125 (MQ=255)
                                                                                                gggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                        >  1:859527/1‑125 (MQ=255)
                                                                                                  gTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCTCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                      >  2:859495/1‑125 (MQ=38)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                  >  1:859546/1‑80 (MQ=18)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                     >  1:859376/1‑125 (MQ=255)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                     >  1:859352/1‑125 (MQ=255)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                     <  1:859273/125‑1 (MQ=255)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCATGACGCGTTCCCCTTCCATGCCCCTc                     >  1:859520/1‑125 (MQ=37)
                                                                                                    aTGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                       >  2:859488/1‑122 (MQ=255)
                                                                                                    aTGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                     <  2:859274/124‑1 (MQ=255)
                                                                                                     tGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCaa                   >  2:859477/1‑125 (MQ=255)
                                                                                                      gCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAATGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                     <  2:859277/122‑1 (MQ=37)
                                                                                                      gCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCaa                   <  2:859276/124‑1 (MQ=255)
                                                                                                        cGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                <  2:859280/125‑1 (MQ=255)
                                                                                                          ttAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctg              >  2:859406/1‑125 (MQ=255)
                                                                                                            aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                        >  1:859596/1‑113 (MQ=255)
                                                                                                            aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                <  2:859282/121‑1 (MQ=255)
                                                                                                              gTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc          <  1:859284/125‑1 (MQ=255)
                                                                                                               ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc         >  1:859503/1‑125 (MQ=255)
                                                                                                               ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc         <  2:859286/125‑1 (MQ=255)
                                                                                                                 ttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGc       <  2:859288/125‑1 (MQ=255)
                                                                                                                    ggCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc      >  1:859461/1‑123 (MQ=255)
                                                                                                                     gcgcCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                 >  2:859561/1‑111 (MQ=255)
                                                                                                                     gcgcCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg   <  2:859290/125‑1 (MQ=255)
                                                                                                                             ||                                                                                                                      
GGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGC  >  NC_002947/1070122‑1070362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: