Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 222,352 A→T 6.6% E69D (GAA→GAT)  PP_0170 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947222,3520AT 6.6% 180.3 / 2.3 60E69D (GAA→GAT) PP_0170ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (34/22);  new base T (2/2);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

TCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGCCGTGGC  >  NC_002947/222204‑222495
                                                                                                                                                    |                                                                                                                                               
tcgGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                                                                                                                 <  7:285606/149‑1 (MQ=255)
  ggCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGt                                                                                                                                               >  8:256217/1‑149 (MQ=255)
         tAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                                                                                                        >  7:89341/1‑149 (MQ=255)
               cccccGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAg                                                                                                                                  <  1:83846/149‑1 (MQ=255)
                        tgacccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGt                                                                                                                         >  5:318849/4‑149 (MQ=255)
                            ccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                                                                                                        >  5:322399/1‑130 (MQ=255)
                            ccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                                                                                                        <  6:322399/130‑1 (MQ=255)
                            ccAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATGGAGAAACACGGCAAGGAGCAGGGCATCGAGATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCt                                                                                                                     >  2:324911/1‑149 (MQ=255)
                              aGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCg                                                                                                                   >  8:302540/1‑149 (MQ=255)
                                  cGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCg                                                                                                               >  6:16915/1‑149 (MQ=255)
                                      ggCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCaa                                                                                                           <  1:272440/149‑1 (MQ=255)
                                         aaGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCt                                                                                                                                >  4:86330/1‑125 (MQ=255)
                                         aaGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCt                                                                                                                                <  3:86330/125‑1 (MQ=255)
                                               cGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAt                                                                                                   <  5:16915/148‑1 (MQ=255)
                                                      ccGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTcc                                                                                           >  2:143915/1‑149 (MQ=255)
                                                              ttCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGt                                                                                   <  8:89341/149‑1 (MQ=255)
                                                              ttCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGt                                                                                   <  7:302540/149‑1 (MQ=255)
                                                                  gCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTAATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGa                                                                                >  7:173849/1‑148 (MQ=255)
                                                                      cGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATcg                                                                           >  2:44157/1‑149 (MQ=255)
                                                                         cTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATcgccg                                                                        >  5:141285/1‑149 (MQ=255)
                                                                             ctgctgctgAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACTGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATcgccggcgc                                                                    >  2:287018/1‑149 (MQ=255)
                                                                              tgctgctgAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATcgccggcgcc                                                                   >  8:56599/1‑149 (MQ=255)
                                                                               gctgctgAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATcgccggcgccg                                                                  >  4:205787/1‑149 (MQ=255)
                                                                                  gctgAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGgtg                                                               >  3:122854/1‑149 (MQ=255)
                                                                                                     cAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGt                                            >  6:131099/1‑149 (MQ=255)
                                                                                                         aCCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTgg                                        >  1:21023/1‑149 (MQ=255)
                                                                                                          ccTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTggg                                       <  1:287018/149‑1 (MQ=255)
                                                                                                           cTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTc                                                                                            >  6:17158/1‑95 (MQ=255)
                                                                                                           cTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTc                                                                                            <  5:17158/95‑1 (MQ=255)
                                                                                                           cTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGa                                      <  3:205787/149‑1 (MQ=255)
                                                                                                           cTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCAGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGa                                      >  3:108351/1‑149 (MQ=255)
                                                                                                           cTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGACGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTggg                                       >  1:266652/1‑148 (MQ=255)
                                                                                                            tGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGAc                                     >  7:131401/1‑149 (MQ=255)
                                                                                                               ttGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGc                                  <  6:141285/149‑1 (MQ=255)
                                                                                                                 gagaAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCAc                                >  4:8176/1‑149 (MQ=255)
                                                                                                                   gaAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCAcc                               >  4:107388/1‑148 (MQ=255)
                                                                                                                      acacGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAgg                           <  1:143915/149‑1 (MQ=255)
                                                                                                                      acacGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAgg                           >  3:241464/1‑149 (MQ=255)
                                                                                                                      acacGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAgg                           >  7:166184/1‑149 (MQ=255)
                                                                                                                        acGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGc                         >  3:119273/1‑149 (MQ=255)
                                                                                                                          ggCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCg                       <  3:107388/149‑1 (MQ=255)
                                                                                                                            cAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCgg                      >  6:32072/1‑148 (MQ=255)
                                                                                                                            cAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGc                     <  2:88649/149‑1 (MQ=255)
                                                                                                                             aaGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGcc                        >  3:9483/1‑145 (MQ=255)
                                                                                                                             aaGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGcc                        <  4:9483/145‑1 (MQ=255)
                                                                                                                             aaGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCa                    <  4:119273/149‑1 (MQ=255)
                                                                                                                                 aGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAAc                <  4:122854/149‑1 (MQ=255)
                                                                                                                                     gggCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGAACGTCTGGGACCGCACCAAGGGCCGGCAGAACGTg             >  8:148515/1‑148 (MQ=255)
                                                                                                                                     gggCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGCACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTg             >  4:31195/1‑148 (MQ=255)
                                                                                                                                      ggCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCAGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGaa           <  4:108351/149‑1 (MQ=255)
                                                                                                                                      ggCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCGCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGaa           <  3:8176/149‑1 (MQ=255)
                                                                                                                                      ccaatcgacttcgaTGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGc                     >  1:61203/4‑139 (MQ=255)
                                                                                                                                      ccaatcgacttcgaTGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGc                     <  2:61203/136‑1 (MQ=255)
                                                                                                                                      ccaatcgacttcgaTGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTg             >  6:167455/4‑147 (MQ=255)
                                                                                                                                        catcgacatcgaAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAgg         <  5:131099/149‑1 (MQ=255)
                                                                                                                                         atcgacatcgaAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGGGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGc        <  6:318849/149‑1 (MQ=255)
                                                                                                                                           cgacatcgaAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTc                                                             >  1:11572/1‑94 (MQ=255)
                                                                                                                                           cgacatcgaAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTc                                                             <  2:11572/94‑1 (MQ=255)
                                                                                                                                            gacttcgaTGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGCCGt     <  5:167455/149‑1 (MQ=255)
                                                                                                                                             acatcgaAGACGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCCCCGGCCCCGGTGGCGGACCCCTGTTGCCCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGCCGTg    >  6:257825/1‑149 (MQ=255)
                                                                                                                                               atcgaAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGCCGTGGc  >  5:317010/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCCGTATTGCCGAACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGGTTCGGTGGACATCGCCGGCGCCGGTGTCGGACCGCTGTTGACCGTCTGGGACCGCACCAAGGGCCGGCAGAACGTGAAGGCCGTGGC  >  NC_002947/222204‑222495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑