Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,747,455 C→A 6.1% A173A (GCG→GCT)  PP_5043 ← PhoPQ‑activated pathogenicity‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,747,4550CA 6.1% 212.6 / 3.3 66A173A (GCG→GCT) PP_5043PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base A (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGGCAATG  >  NC_002947/5747307‑5747602
                                                                                                                                                    |                                                                                                                                                   
ccATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGgc                                                                                                                                                     >  2:248857/1‑149 (MQ=255)
       cggcAGCCATGCACTCCACCCGCGTTTCGAGGCAGCGGTGCTGATGAACTGCTCGAGGCGGAGCTCGGGCAACGCGGCAGCCAGGAGGTCCATGCCCCTGCTGGCGCCAGCGGCGATGGGGACAGGCCGCGGCAGCGCCGGCGccgccg                                                                                                                                              >  8:267850/4‑149 (MQ=255)
        gTCAGCCATGCACTACAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGc                                                                                                                                             >  3:333010/1‑149 (MQ=255)
                   aCTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACg                                                                                                                                    >  1:23178/1‑147 (MQ=255)
                   aCTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCa                                                                                                                                  <  6:363160/149‑1 (MQ=255)
                     tCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGc                                                                                                                                >  7:148512/1‑149 (MQ=255)
                     tCCAACAGCGTTTCGAGGCACAGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGATGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGc                                                                                                                                >  7:315321/1‑149 (MQ=255)
                        aaCCGCGTTTCGGGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCaa                                                                                                                             >  8:69280/1‑149 (MQ=255)
                        aaCCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCaa                                                                                                                             >  7:349864/1‑149 (MQ=255)
                         aCCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAAt                                                                                                                            >  2:322215/1‑149 (MQ=255)
                               tttCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGcc                                                                                                                       >  6:266211/1‑148 (MQ=255)
                                     ggCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTAt                                                                                                                <  1:192048/149‑1 (MQ=255)
                                       cACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATgg                                                                                                              <  2:23178/149‑1 (MQ=255)
                                       cACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATgg                                                                                                              <  1:248857/149‑1 (MQ=255)
                                        aCCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATggg                                                                                                             >  3:143759/1‑149 (MQ=255)
                                               atgatgaACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACgg                                                                                                       <  7:228349/148‑1 (MQ=255)
                                                tgatgaACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGTGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGca                                                                                                     <  7:87741/149‑1 (MQ=255)
                                                 gatgaACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGcat                                                                                                    >  3:264641/1‑149 (MQ=255)
                                                   tgaACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGcatca                                                                                                  >  1:313849/1‑149 (MQ=255)
                                                          cTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTTGCGACATCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTcgc                                                                                           <  4:143759/149‑1 (MQ=255)
                                                          cTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGTCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTcgc                                                                                           <  8:315321/149‑1 (MQ=255)
                                                          cTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTc                                                                                             <  8:287447/147‑1 (MQ=255)
                                                           tCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTcgcg                                                                                          >  4:142575/1‑149 (MQ=255)
                                                           tCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTcgcg                                                                                          >  5:148845/1‑149 (MQ=255)
                                                               ggCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATc                                                                                      <  2:44091/149‑1 (MQ=255)
                                                                            cAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCAt                                                                         >  1:167483/1‑149 (MQ=255)
                                                                             aaCGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCAt                                                                         >  7:215674/1‑148 (MQ=255)
                                                                             aaCGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCAt                                                                         >  7:209969/1‑148 (MQ=255)
                                                                                   tCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACgg                                                                    >  4:83093/1‑147 (MQ=255)
                                                                                   tCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGaa                                                                  <  4:333010/149‑1 (MQ=255)
                                                                                     aGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAgg                                                                >  1:231211/1‑149 (MQ=255)
                                                                                       ccATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGt                                                               >  6:210359/1‑148 (MQ=255)
                                                                                          tGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTACATCACGGAAGGTcaca                                                           <  8:349864/149‑1 (MQ=255)
                                                                                             ggTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATAc                                                        <  3:83093/149‑1 (MQ=255)
                                                                                              gTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACt                                                       >  5:52165/1‑149 (MQ=255)
                                                                                              gTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACt                                                       >  3:324713/1‑149 (MQ=255)
                                                                                                 cATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGt                                                    >  5:174325/1‑149 (MQ=255)
                                                                                                   tGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCTCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTg                                                  <  5:196082/149‑1 (MQ=255)
                                                                                                   tGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTg                                                  <  2:313849/149‑1 (MQ=255)
                                                                                                     ccccTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTgg                                                 <  5:266211/148‑1 (MQ=255)
                                                                                                     ccccTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTgg                                                 <  3:142575/148‑1 (MQ=255)
                                                                                                      cccTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTggg                                                >  6:60471/1‑148 (MQ=255)
                                                                                                      ccatgctgGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTggg                                                >  4:262394/4‑148 (MQ=255)
                                                                                                              gCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGc                                        <  6:52165/148‑1 (MQ=255)
                                                                                                              gCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGc                                        <  5:210359/148‑1 (MQ=255)
                                                                                                              gCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCt                                       <  7:267850/149‑1 (MQ=255)
                                                                                                                gACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTg                                      >  2:159623/1‑148 (MQ=255)
                                                                                                                     cGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGAt                                <  6:148845/149‑1 (MQ=255)
                                                                                                                     cGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGCGCACATCGCTGATccct                                >  7:361444/1‑146 (MQ=255)
                                                                                                                       gCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGTACATCGCTGATCGATGc                              <  4:264641/149‑1 (MQ=255)
                                                                                                                          aTTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCACAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCca                            >  5:7991/1‑148 (MQ=255)
                                                                                                                           ttGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCcac                           >  7:282991/1‑148 (MQ=255)
                                                                                                                              gCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCgcgc                                                                                         >  3:168326/1‑83 (MQ=255)
                                                                                                                              gCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCgcgc                                                                                         <  4:168326/83‑1 (MQ=255)
                                                                                                                              gCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCacaaca                       <  1:322215/149‑1 (MQ=255)
                                                                                                                                 caTGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTc                    <  4:324713/149‑1 (MQ=255)
                                                                                                                                     cAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCtgtg                >  4:207925/1‑149 (MQ=255)
                                                                                                                                     cAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCtgt                 <  3:262394/148‑1 (MQ=255)
                                                                                                                                      aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                             >  8:39978/1‑87 (MQ=38)
                                                                                                                                      aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                             <  7:314496/87‑1 (MQ=38)
                                                                                                                                      aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                             >  8:314496/1‑87 (MQ=38)
                                                                                                                                      aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                             <  7:39978/87‑1 (MQ=38)
                                                                                                                                      aGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCTGCCAATGCGCCCAGGTTTGGGCCACGGCATCATTCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTg                                                  <  7:247362/114‑1 (MQ=255)
                                                                                                                                      aGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTgtgt               <  2:167483/149‑1 (MQ=255)
                                                                                                                                            aGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCAc         <  2:231211/149‑1 (MQ=255)
                                                                                                                                               tCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACgg       >  2:11712/1‑148 (MQ=255)
                                                                                                                                                 cGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGGCaa    >  8:335321/1‑149 (MQ=255)
                                                                                                                                                   gcgcAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGGCAATg  <  8:282991/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGGCAATG  >  NC_002947/5747307‑5747602

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: ‑