Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,086 A→G 7.8% T854A (ACC→GCC)  PP_4704 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,0860AG 7.8% 176.5 / 4.8 64T854A (ACC→GCC) PP_4704hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

AGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGTT  >  NC_002947/5349939‑5350234
                                                                                                                                                   |                                                                                                                                                     
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)
                           aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACggg‑g                                                                                                                          <  3:352229/149‑1 (MQ=255)
                             gCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑Gt                                                                                                                         >  4:287025/1‑148 (MQ=255)
                               tGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCAACGGG‑GTGa                                                                                                                       >  5:313864/1‑148 (MQ=255)
                                 gCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGAc                                                                                                                      <  5:342503/147‑1 (MQ=255)
                                  cATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAg                                                                                                                   <  3:287025/149‑1 (MQ=255)
                                    tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGCGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGt                                                                                                                 >  1:187576/1‑149 (MQ=255)
                                      cgcgCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGc                                                                                                               >  4:234086/1‑149 (MQ=255)
                                              aaGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCCACGGG‑GTGACAAGGTGCCGGCCAAt                                                                                                       <  6:313864/149‑1 (MQ=255)
                                                  gcgcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCt                                                                                                   <  8:148663/149‑1 (MQ=255)
                                                    gcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTc                                                                                                  >  5:268782/1‑148 (MQ=255)
                                                           aCTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtgg                                                                                          <  2:187576/149‑1 (MQ=255)
                                                               cgcCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggt                                                                                      <  1:272387/149‑1 (MQ=255)
                                                                 ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtg                                                                                     >  3:197703/1‑148 (MQ=255)
                                                                 ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑CTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                    >  2:135929/1‑149 (MQ=255)
                                                                  cTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGAAAGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGga                                                                                   >  1:297675/1‑149 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                             >  2:182855/1‑106 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                             <  1:182855/106‑1 (MQ=255)
                                                                    ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgac                                                                                  >  5:286901/1‑148 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                           >  3:296300/1‑135 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                           <  4:296300/135‑1 (MQ=255)
                                                                        aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAg                                                                             <  6:286901/149‑1 (MQ=255)
                                                                                 aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                         <  1:284814/144‑1 (MQ=255)
                                                                                 aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                         >  2:284814/1‑144 (MQ=255)
                                                                                  ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCaa                                                                                                        <  2:306699/112‑1 (MQ=255)
                                                                                  ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCaa                                                                                                        >  1:306699/1‑112 (MQ=255)
                                                                                     gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                                  >  4:159427/1‑147 (MQ=255)
                                                                                     gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                                  <  3:159427/147‑1 (MQ=255)
                                                                                          cAGGGCACCGTCGGGCCGCTGATTGCCCAACAGCGCCTTGAGCGTTATTTCCACGGCACCTTGTAAAATATCCTGCATCGTCGGGTGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGc                                                            <  1:106800/149‑1 (MQ=255)
                                                                                              ccAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑CTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGaa                                                       <  1:135929/149‑1 (MQ=255)
                                                                                                gtcgctggCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAgg                                                     <  1:12861/143‑1 (MQ=255)
                                                                                                   gACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAt                                                  <  2:297675/149‑1 (MQ=255)
                                                                                                    aCGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAtt                                                 >  2:176238/1‑149 (MQ=255)
                                                                                                                   cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                             <  4:265938/138‑1 (MQ=255)
                                                                                                                   cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                             >  3:265938/1‑138 (MQ=255)
                                                                                                                    ccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTCACCGTACCCGCCTg                                  >  8:43274/1‑148 (MQ=255)
                                                                                                                     cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                          <  6:161007/107‑1 (MQ=255)
                                                                                                                     cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                          >  5:161007/1‑107 (MQ=255)
                                                                                                                      aaCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                               <  6:268782/149‑1 (MQ=255)
                                                                                                                       aCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                               >  6:99179/1‑148 (MQ=255)
                                                                                                                          gcgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTTTATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagc                           <  1:112155/149‑1 (MQ=255)
                                                                                                                           cgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagca                          <  5:99179/149‑1 (MQ=255)
                                                                                                                              cTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgcc                       >  5:102305/1‑149 (MQ=255)
                                                                                                                                     gTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGcc                >  5:146704/1‑149 (MQ=255)
                                                                                                                                        tgTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGaa             >  4:170463/2‑149 (MQ=255)
                                                                                                                                         gTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt            >  1:225603/1‑149 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGCTCCTGCAGCGACGGG‑GTGACAACGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt            <  8:260976/149‑1 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                               >  3:227441/1‑130 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                               <  4:227441/130‑1 (MQ=255)
                                                                                                                                         gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt            <  6:146704/149‑1 (MQ=255)
                                                                                                                                          tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg           >  1:236915/1‑149 (MQ=255)
                                                                                                                                          tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg           <  1:176238/149‑1 (MQ=255)
                                                                                                                                            ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATgg          <  4:197703/148‑1 (MQ=255)
                                                                                                                                              aGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAg       <  2:236915/149‑1 (MQ=255)
                                                                                                                                               ggCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTCCCGTCCCCGCCTGAACCAGCAGCCGCCAGCCGAATGGGAg       >  7:312541/1‑148 (MQ=255)
                                                                                                                                                gCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                    >  3:108577/1‑70 (MQ=255)
                                                                                                                                                gCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                    <  4:108577/70‑1 (MQ=255)
                                                                                                                                                  cACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGt   >  4:160582/1‑149 (MQ=255)
                                                                                                                                                   aCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGa        <  3:80129/143‑1 (MQ=255)
                                                                                                                                                   aCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGa        >  4:80129/1‑143 (MQ=255)
                                                                                                                                                   aCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGtt  >  3:109671/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAGTT  >  NC_002947/5349939‑5350234

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