Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,263 G→C 13.8% intergenic (+63/+79) 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,2630GC13.8% 143.0 / 16.1 58intergenic (+63/+79)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base C (4/4);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

TCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑T‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATA  >  NC_002947/1070118‑1070410
                                                                                                                                                       |                                                                                                                                                   
tctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCgcggg                                                                                                                                                    >  8:290400/1‑149 (MQ=255)
        aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                             <  4:126850/148‑1 (MQ=255)
        aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATACCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            <  4:138923/149‑1 (MQ=255)
        aTGCAGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCGAAGGGCTGATCGGGTTTGGTCGTAAGGGGTATGCCGTTAGGTTTGCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            <  1:249782/149‑1 (MQ=38)
           cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         <  7:240936/149‑1 (MQ=255)
           cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                          <  4:278848/148‑1 (MQ=38)
              aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                                      >  3:203589/1‑149 (MQ=255)
               tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     <  1:318166/149‑1 (MQ=255)
               tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     <  1:313083/149‑1 (MQ=255)
                cccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     >  5:293296/1‑148 (MQ=255)
                  cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                                                  >  7:332143/1‑149 (MQ=255)
                   gggaTGGATCGCTGTGCTGTGGGTGTCTTTCAAGCTTAATAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                 <  6:293296/146‑1 (MQ=38)
                    gCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGc                                                                                                                                <  5:133706/149‑1 (MQ=255)
                           tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                         <  8:106916/149‑1 (MQ=255)
                           tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                         <  7:122536/149‑1 (MQ=255)
                               tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  2:266753/1‑126 (MQ=38)
                               tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            <  1:266753/126‑1 (MQ=38)
                                      tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGGGATGAAAGCCTAAAGGc                                                                                                               <  8:332143/148‑1 (MQ=255)
                                              tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCGGCGCGGGCGGCGCTCGAt                                                                                                                                         >  1:283246/1‑114 (MQ=37)
                                              tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         <  2:283246/114‑1 (MQ=255)
                                              tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                                                      >  3:267793/1‑149 (MQ=255)
                                                         agaagGCTCGCCAGGCCAAGGGCTGATCGGGTGTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                           <  7:290400/149‑1 (MQ=255)
                                                           aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGCAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                                                         >  4:41170/1‑149 (MQ=255)
                                                             ggCTCGCCAGGCCTAGTGCTGTTCCGGTTTGTTCGTTAGGGGTTTGCCGTTAGGTTTTCGGCGCCTTTG‑‑T‑‑GATCGATCGCCGCGGGGGCGGCGCTCGATCTCACAGGCGATGACAGCCTAAAGGCATGCGTCTggc                                                                                                    <  7:222100/136‑1 (MQ=25)
                                                              gCTCGCCAGGCCAAGGGCTGATCGGGTTTCCTCCTAAGGGCTAACCCGTTAGGCTTTCGGCGCC‑‑TGACA‑‑GACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAGGAAAGc                                                                                                                       <  7:164967/116‑1 (MQ=18)
                                                               cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 >  7:255739/1‑89 (MQ=37)
                                                               cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 <  8:255739/89‑1 (MQ=37)
                                                                tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                                                    >  7:4252/1‑149 (MQ=255)
                                                                  gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                                                    <  8:4252/147‑1 (MQ=255)
                                                                      ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                                      >  5:95351/1‑93 (MQ=25)
                                                                      ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                                      <  6:95351/93‑1 (MQ=25)
                                                                        ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                                                            >  8:278122/1‑149 (MQ=255)
                                                                              gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCA‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                                                      >  6:16188/1‑149 (MQ=255)
                                                                                 gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          <  1:333812/78‑1 (MQ=21)
                                                                                 gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          >  2:333812/1‑78 (MQ=21)
                                                                                 gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCaa                                                                   >  7:70731/1‑149 (MQ=255)
                                                                                  agaggggTTGTTCGTAAGGGGTATGCCGGTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                                                  <  4:267793/146‑1 (MQ=38)
                                                                                     gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCA‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                                                                <  5:16188/148‑1 (MQ=255)
                                                                                            ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                                                         <  2:10392/148‑1 (MQ=255)
                                                                                               gTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                      >  3:133218/1‑148 (MQ=255)
                                                                                                tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    <  2:27290/149‑1 (MQ=255)
                                                                                                    gggTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                                                <  3:41170/149‑1 (MQ=255)
                                                                                                       tATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                                             >  5:66103/1‑149 (MQ=255)
                                                                                                          gCCGTTAGGTTTTCGGCGCCCGTG‑‑A‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCGTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAg                                          >  5:212274/1‑149 (MQ=38)
                                                                                                            cGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                        <  7:278122/149‑1 (MQ=255)
                                                                                                            cGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                        <  6:66103/149‑1 (MQ=255)
                                                                                                            cGTTAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                        >  7:60156/1‑149 (MQ=255)
                                                                                                              ttAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    >  7:152952/1‑135 (MQ=255)
                                                                                                              ttAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    <  8:152952/135‑1 (MQ=255)
                                                                                                              ttAGGTTTTCGGCGCC‑‑TG‑‑TGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                      <  4:133218/149‑1 (MQ=255)
                                                                                                                                  tgtgagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTACAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCGTCCATGCCCCTCAATGCTGCCGCCGCCAAGCACGCACAGGGGTGTACCGCCACCCCg                        >  7:178345/6‑147 (MQ=37)
                                                                                                                                   gtgagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAAAGCCACACCGGAg                     >  4:139050/5‑149 (MQ=255)
                                                                                                                                     gagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCg                   <  6:212274/147‑1 (MQ=255)
                                                                                                                                       gATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAacac                          <  7:234948/140‑1 (MQ=255)
                                                                                                                                       gATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAacac                          >  8:234948/1‑140 (MQ=255)
                                                                                                                                       gATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa                  <  8:60156/148‑1 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    >  7:256076/1‑113 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    <  8:256076/113‑1 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCATCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAg                 >  1:274187/1‑148 (MQ=255)
                                                                                                                                             gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                                 >  6:136237/1‑127 (MQ=255)
                                                                                                                                             gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                                 <  5:136237/127‑1 (MQ=255)
                                                                                                                                                     cgggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa    >  1:130015/1‑148 (MQ=255)
                                                                                                                                                     cgggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACAt   <  3:139050/149‑1 (MQ=255)
                                                                                                                                                      gggcggCGCTCGATCTCACAGGCGCTGAAAGCCGAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATa  <  8:178345/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TG‑‑T‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATA  >  NC_002947/1070118‑1070410

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: