Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,085 C→G 7.9% A853A (GCC→GCG) ‡ PP_4704 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,0850CG 7.9% 191.2 / 6.4 63A853A (GCC→GCG) ‡PP_4704hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

AGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGT  >  NC_002947/5349939‑5350233
                                                                                                                                                  |                                                                                                                                                     
aGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCACCAGCGCCTTGAGCGTGAGTACCAGGCCAc                                                                                                                                                     >  2:112155/1‑149 (MQ=255)
          cccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActca                                                                                                                                           >  7:70241/1‑146 (MQ=255)
           ccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActcac                                                                                                                                          <  8:70241/149‑5 (MQ=255)
             aGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActcacgc                                                                                                                                        <  6:236377/149‑7 (MQ=255)
                  tcCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCACCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCa                                                                                                                                   <  8:252525/149‑1 (MQ=255)
                    cTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGc                                                                                                                                 >  7:105145/1‑149 (MQ=255)
                          cAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACggg                                                                                                                           <  8:105145/149‑1 (MQ=255)
                           aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgggg                                                                                                                         <  3:352229/149‑1 (MQ=255)
                             gCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGt                                                                                                                        >  4:287025/1‑148 (MQ=255)
                               tGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCAACGGGGTGa                                                                                                                      >  5:313864/1‑148 (MQ=255)
                                 gCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGAc                                                                                                                     <  5:342503/147‑1 (MQ=255)
                                  cATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAg                                                                                                                  <  3:287025/149‑1 (MQ=255)
                                    tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGCGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGt                                                                                                                >  1:187576/1‑149 (MQ=255)
                                      cgcgCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGc                                                                                                              >  4:234086/1‑149 (MQ=255)
                                              aaGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCCACGGGGTGACAAGGTGCCGGCCAAt                                                                                                      <  6:313864/149‑1 (MQ=255)
                                                  gcgcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCt                                                                                                  <  8:148663/149‑1 (MQ=255)
                                                    gcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTc                                                                                                 >  5:268782/1‑148 (MQ=255)
                                                           aCTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtgg                                                                                         <  2:187576/149‑1 (MQ=255)
                                                               cgcCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggt                                                                                     <  1:272387/149‑1 (MQ=255)
                                                                 ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtg                                                                                    >  3:197703/1‑148 (MQ=255)
                                                                 ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGCTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                   >  2:135929/1‑149 (MQ=255)
                                                                  cTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGAAAGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGga                                                                                  >  1:297675/1‑149 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                            >  2:182855/1‑106 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                            <  1:182855/106‑1 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgac                                                                                 >  5:286901/1‑148 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                          <  4:296300/135‑1 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                          >  3:296300/1‑135 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAg                                                                            <  6:286901/149‑1 (MQ=255)
                                                                                 aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                        <  1:284814/144‑1 (MQ=255)
                                                                                 aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                        >  2:284814/1‑144 (MQ=255)
                                                                                  ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCaa                                                                                                       >  1:306699/1‑112 (MQ=255)
                                                                                  ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCaa                                                                                                       <  2:306699/112‑1 (MQ=255)
                                                                                     gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                                 <  3:159427/147‑1 (MQ=255)
                                                                                     gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                                 >  4:159427/1‑147 (MQ=255)
                                                                                          cAGGGCACCGTCGGGCCGCTGATTGCCCAACAGCGCCTTGAGCGTTATTTCCACGGCACCTTGTAAAATATCCTGCATCGTCGGGTGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGc                                                           <  1:106800/149‑1 (MQ=255)
                                                                                              ccAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGCTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGaa                                                      <  1:135929/149‑1 (MQ=255)
                                                                                                gtcgctggCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAgg                                                    <  1:12861/143‑1 (MQ=255)
                                                                                                   gACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAt                                                 <  2:297675/149‑1 (MQ=255)
                                                                                                    aCGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAtt                                                >  2:176238/1‑149 (MQ=255)
                                                                                                                   cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                            >  3:265938/1‑138 (MQ=255)
                                                                                                                   cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                            <  4:265938/138‑1 (MQ=255)
                                                                                                                    ccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTCACCGTACCCGCCTg                                 >  8:43274/1‑148 (MQ=255)
                                                                                                                     cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                         <  6:161007/107‑1 (MQ=255)
                                                                                                                     cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                         >  5:161007/1‑107 (MQ=255)
                                                                                                                      aaCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                              <  6:268782/149‑1 (MQ=255)
                                                                                                                       aCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                              >  6:99179/1‑148 (MQ=255)
                                                                                                                          gcgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTTTATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagc                          <  1:112155/149‑1 (MQ=255)
                                                                                                                           cgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagca                         <  5:99179/149‑1 (MQ=255)
                                                                                                                              cTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgcc                      >  5:102305/1‑149 (MQ=255)
                                                                                                                                     gTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGcc               >  5:146704/1‑149 (MQ=255)
                                                                                                                                        tgTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGaa            >  4:170463/2‑149 (MQ=255)
                                                                                                                                         gTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt           >  1:225603/1‑149 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGCTCCTGCAGCGACGGGGTGACAACGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt           <  8:260976/149‑1 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                              <  4:227441/130‑1 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                              >  3:227441/1‑130 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt           <  6:146704/149‑1 (MQ=255)
                                                                                                                                          tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg          <  1:176238/149‑1 (MQ=255)
                                                                                                                                          tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg          >  1:236915/1‑149 (MQ=255)
                                                                                                                                            ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATgg         <  4:197703/148‑1 (MQ=255)
                                                                                                                                              aGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAg      <  2:236915/149‑1 (MQ=255)
                                                                                                                                               ggCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTCCCGTCCCCGCCTGAACCAGCAGCCGCCAGCCGAATGGGAg      >  7:312541/1‑148 (MQ=255)
                                                                                                                                                gCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                   >  3:108577/1‑70 (MQ=255)
                                                                                                                                                gCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                   <  4:108577/70‑1 (MQ=255)
                                                                                                                                                  cACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGt  >  4:160582/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGT  >  NC_002947/5349939‑5350233

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