Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,747,450 T→G 7.1% E175A (GAG→GCG)  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,4500TG 7.1% 196.7 / 4.4 66E175A (GAG→GCG) PP_5043PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base G (3/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGG  >  NC_002947/5747305‑5747597
                                                                                                                                                 |                                                                                                                                                   
ggCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTct                                                                                                                                                   <  4:56595/148‑2 (MQ=255)
  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:209969/1‑148 (MQ=255)
                                                                               aaCGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCAt                                                                    >  7:215674/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                                                  >  3:324713/1‑149 (MQ=255)
                                                                                                gTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACt                                                  >  5:52165/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                                                                        <  7:39978/87‑1 (MQ=38)
                                                                                                                                        aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                        >  8:314496/1‑87 (MQ=38)
                                                                                                                                        aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                        <  7:314496/87‑1 (MQ=38)
                                                                                                                                        aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                        >  8:39978/1‑87 (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)
                                                                                                                                                 |                                                                                                                                                   
GGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCACGG  >  NC_002947/5747305‑5747597

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