Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,747,447 A→T 6.0% L176Q (CTG→CAG)  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,4470AT 6.0% 198.6 / 2.6 67L176Q (CTG→CAG) PP_5043PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base T (2/2);  total (35/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 = 7.95e-01

TCATCGGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCAC  >  NC_002947/5747300‑5747595
                                                                                                                                                   |                                                                                                                                                    
tCATCGGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCCTGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAg                                                                                                                                                     >  8:228349/1‑149 (MQ=255)
tCATCGGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACAGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAg                                                                                                                                                     <  4:51397/149‑1 (MQ=255)
     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                                                                                                       <  1:248857/149‑1 (MQ=255)
                                              cACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATgg                                                                                                       <  2:23178/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                                          <  3:142575/148‑1 (MQ=255)
                                                                                                            ccccTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTgg                                          <  5:266211/148‑1 (MQ=255)
                                                                                                             cccTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTggg                                         >  6:60471/1‑148 (MQ=255)
                                                                                                             ccatgctgGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTggg                                         >  4:262394/4‑148 (MQ=255)
                                                                                                                     gCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGc                                 <  5:210359/148‑1 (MQ=255)
                                                                                                                     gCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGc                                 <  6:52165/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                                                                                  <  4:168326/83‑1 (MQ=255)
                                                                                                                                     gCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCgcgc                                                                                  >  3:168326/1‑83 (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:314496/87‑1 (MQ=38)
                                                                                                                                             aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                      <  7:39978/87‑1 (MQ=38)
                                                                                                                                             aGCGGCTGCGCCGGAGCTGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTtgt                                                                      >  8:314496/1‑87 (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)
                                                                                                                                                   |                                                                                                                                                    
TCATCGGCCATGGCCGTCAGCCATGCACTCCAACCGCGTTTCGAGGCACCGGTGATGATGAACTGCTCGAGGCGTAGCTCGGGCAACGCGTCAGCCATGAGGTCCATGCCCCTGCTGGCGACAGCGGCCATTGGCACATGCAGCGGCAGCTCCGGCGCAGCCGCTTCGTTACGCAGCCAATGCGCCCAGGTATGGGCCACGGCATCATCCTCGCGCATCGGCGCTTGTGCATCACGGAAGGTCACATACTGGTTGGGCACATCGCTGATCGATGCCACAACAGTCTGTGTCTGCAC  >  NC_002947/5747300‑5747595

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