Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,265 G→C 17.1% intergenic (+62/+82) 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,2650GC17.1% 210.7 / 36.5 82intergenic (+62/+82)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base C (7/7);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACA  >  NC_002947/1070118‑1070413
                                                                                                                                                    |                                                                                                                                                    
acATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCcgcgcg                                                                                                                                                    >  1:1412872/2‑150 (MQ=255)
  aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTTAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                                                                  <  2:1338504/150‑1 (MQ=255)
    ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                <  2:203299/150‑1 (MQ=255)
       ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                             >  2:976177/1‑150 (MQ=255)
          aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                          <  2:997302/150‑1 (MQ=255)
           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                                         >  1:322239/1‑150 (MQ=255)
            tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                        <  2:1732237/150‑1 (MQ=255)
            tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCACCGTTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                        <  1:460305/150‑1 (MQ=255)
              cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                      <  2:1209165/150‑1 (MQ=255)
              cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCAGATCGGGTTTGTTCGTAAGGGGGATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCTCCGAGCGCGGGGGGCTCgtcc                                                                                                                                      >  2:341144/1‑147 (MQ=255)
               tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                     <  1:824481/150‑1 (MQ=255)
               tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                      <  1:1076140/149‑1 (MQ=255)
               tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                      >  2:1076140/1‑149 (MQ=255)
                     cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                      >  2:1863969/1‑143 (MQ=255)
                     cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                      <  1:1863969/143‑1 (MQ=255)
                       gCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                   <  2:1871718/144‑1 (MQ=255)
                       gCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                   >  1:1871718/1‑144 (MQ=255)
                          tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGACGGCGGCGCACGGTCTCCCAGGCGagg                                                                                                                          >  2:1080326/1‑148 (MQ=255)
                             aTCGCTTTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGTCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       <  1:341144/150‑1 (MQ=255)
                               cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                     <  2:1801189/150‑1 (MQ=255)
                               cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                     <  2:1206773/150‑1 (MQ=255)
                                 ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                                                   >  2:819823/1‑150 (MQ=255)
                                  tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCGAGAGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                                                                  >  1:1416439/1‑150 (MQ=255)
                                       tgcgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGc                                                                                                             >  1:1940363/4‑150 (MQ=255)
                                         tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                               <  2:54201/114‑1 (MQ=255)
                                         tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                               >  1:54201/1‑114 (MQ=255)
                                          gggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                                                          >  2:1272618/1‑150 (MQ=255)
                                           ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTGTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATGGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                                         >  1:934907/1‑150 (MQ=255)
                                             tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGt                                                                                                       <  1:1284134/150‑1 (MQ=255)
                                                cTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                     >  2:954834/1‑133 (MQ=255)
                                                cTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                     <  1:954834/133‑1 (MQ=255)
                                                 tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                                                   >  1:1326562/1‑150 (MQ=255)
                                                 tttAAAGCTGAAGAAGGCTCGCAAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGTGGCGCTCGATATGACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                                                   >  1:793029/1‑150 (MQ=39)
                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                                                           >  1:751823/1‑150 (MQ=255)
                                                           aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCTATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGTGGCCCTGAc                                                                                         >  2:1547521/1‑150 (MQ=255)
                                                           aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         >  1:1597617/1‑150 (MQ=255)
                                                           aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         >  2:1239780/1‑150 (MQ=255)
                                                           aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         <  2:1242711/150‑1 (MQ=255)
                                                           aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         >  1:1347260/1‑150 (MQ=255)
                                                           aagaagGCGCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCCGTAACTGAAGTTCGTAGCGCCGCGCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         <  2:558697/150‑1 (MQ=255)
                                                             gaagGCTCGCTAGTCCAATGGCTTATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGAGTGACAATGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                                                       <  1:612040/150‑1 (MQ=37)
                                                             gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                                                       <  1:1080326/150‑1 (MQ=255)
                                                               agGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgg                                                                                                                                                   >  1:667735/1‑88 (MQ=255)
                                                               agGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgg                                                                                                                                                   <  2:667735/88‑1 (MQ=255)
                                                                 gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                                                   <  1:1828026/150‑1 (MQ=255)
                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       >  2:226439/1‑112 (MQ=255)
                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       <  1:226439/112‑1 (MQ=255)
                                                                     gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                                                               >  1:273999/1‑150 (MQ=255)
                                                                      ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                                              >  1:385125/1‑150 (MQ=255)
                                                                       cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                             >  1:321710/1‑150 (MQ=255)
                                                                         ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                                                           <  1:1272618/150‑1 (MQ=255)
                                                                         ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                                                           >  2:173732/1‑150 (MQ=255)
                                                                           ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                                                         <  2:273999/150‑1 (MQ=255)
                                                                            cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                                                        >  1:211810/1‑150 (MQ=255)
                                                                             aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                                                                       <  2:1347260/150‑1 (MQ=255)
                                                                              aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGTAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                                                      <  2:10879/150‑1 (MQ=255)
                                                                                 gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCGGACGAGTTCCCCTTCCATGCCCCTc                                                                   >  2:680858/1‑150 (MQ=255)
                                                                                    gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                                                <  2:385125/150‑1 (MQ=255)
                                                                                     aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAGGCCCCTCAAtg                                                               >  2:717300/1‑150 (MQ=255)
                                                                                       cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                             <  2:1412872/150‑1 (MQ=255)
                                                                                       cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                             <  1:1607091/150‑1 (MQ=255)
                                                                                         ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                                           <  2:751823/150‑1 (MQ=255)
                                                                                          gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGAGAGATCGGGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgca                                                          <  2:1326562/150‑2 (MQ=255)
                                                                                           tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCTGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCt                                                         <  1:819823/150‑1 (MQ=255)
                                                                                                     aGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgca                                               <  1:1239780/150‑1 (MQ=255)
                                                                                                      ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                                              >  2:501255/1‑150 (MQ=255)
                                                                                                            tggcGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGACCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAg                                        >  1:1937223/4‑150 (MQ=255)
                                                                                                               cGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg                                     >  2:812892/1‑150 (MQ=255)
                                                                                                               cGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                          >  1:1687785/1‑49 (MQ=11)
                                                                                                               cGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                          <  2:1687785/49‑1 (MQ=11)
                                                                                                                gTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAACGACTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                    <  1:122101/150‑1 (MQ=255)
                                                                                                                     gTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                               >  1:864784/1‑150 (MQ=255)
                                                                                                                          catcgcCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAac                          <  2:1016918/147‑1 (MQ=255)
                                                                                                                                cTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTa                    <  2:321710/150‑1 (MQ=255)
                                                                                                                                  gtgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAcccccgagc                  >  2:684320/1‑145 (MQ=255)
                                                                                                                                      gaTCGAGCGCCGCCCGCGCGGCGCTCAATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGc              >  2:1388814/1‑150 (MQ=255)
                                                                                                                                       aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgcg                                                                                      >  1:414471/1‑77 (MQ=255)
                                                                                                                                       aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgcg                                                                                      <  2:414471/77‑1 (MQ=255)
                                                                                                                                       aTCGAGCGCCGCCCGCGCGGCGCTCGCTCTCACAGACGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCa             >  2:711419/1‑150 (MQ=255)
                                                                                                                                         cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAcc           <  1:812892/150‑1 (MQ=255)
                                                                                                                                         cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAcc           >  2:351725/1‑150 (MQ=255)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCTTGCCCCTCaa                                                                 >  1:970062/1‑95 (MQ=37)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCTTGCCCCTCaa                                                                 <  2:970062/95‑1 (MQ=37)
                                                                                                                                            aGCGCCGCGAGGGCGCCGCTCGCACTGACAACCGTTGAAAGGCTAAAGGCATTCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCgg         <  2:211810/150‑1 (MQ=255)
                                                                                                                                            aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTa                    >  2:1537284/1‑139 (MQ=255)
                                                                                                                                            aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTa                    <  1:1537284/139‑1 (MQ=255)
                                                                                                                                              cgccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCa       <  2:864784/150‑1 (MQ=255)
                                                                                                                                               gccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                  >  1:667221/1‑106 (MQ=255)
                                                                                                                                               gccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                  <  2:667221/106‑1 (MQ=255)
                                                                                                                                                 cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTa    <  2:1597617/150‑1 (MQ=255)
                                                                                                                                                   cgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa  <  2:1940363/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACA  >  NC_002947/1070118‑1070413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: