Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,470,854 Δ1 bp 6.5% coding (1678/2460 nt) PP_3084 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,470,8540T. 6.5% 194.5 / 7.1 62coding (1678/2460 nt)PP_3084ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (34/24);  new base . (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GGTGCCGTTGGCCCCTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAGGTA  >  NC_002947/3470709‑3470999
                                                                                                                                                 |                                                                                                                                                  
ggTGCCGTTGGCCCCTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTc                                                                                                                                                 >  3:49102/1‑149 (MQ=255)
  tGCCGTTGGCCCCTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGa                                                                                                                                               <  1:248421/149‑1 (MQ=255)
    ggcctGGCCCAGGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACg                                                                                                                                            >  1:211147/5‑149 (MQ=255)
     gcctGGCCCAGGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGc                                                                                                                                           >  4:187306/4‑149 (MQ=255)
          gCCCCTGGGCCAGGCCCTCGATGCGGGTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATc                                                                                                                                      <  8:77318/149‑1 (MQ=255)
            cccTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGa                                                                                                                                    <  7:214278/149‑1 (MQ=255)
              cTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGAGCGACt                                                                                                                                  >  1:194033/1‑149 (MQ=255)
               tGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTg                                                                                                                                 >  8:96739/1‑149 (MQ=255)
                   ccAGGCCCTCGATGCGGTTGTCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGATCAGCGTCGACAGCTGGGCGACTGccac                                                                                                                             <  6:261769/149‑3 (MQ=38)
                          cTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGAc                                                                                                                                              <  5:32919/126‑1 (MQ=255)
                          cTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGAc                                                                                                                                              >  6:32919/1‑126 (MQ=255)
                                  ggTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTc                                                                                                              <  4:49102/149‑1 (MQ=255)
                                   gTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTc                                                                                                              >  1:133155/1‑148 (MQ=255)
                                    ttGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCaa                                                                                                            <  2:211147/149‑1 (MQ=255)
                                      ggCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAAc                                                                                                           <  6:58006/148‑1 (MQ=255)
                                       gCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACa                                                                                                          <  3:187306/148‑1 (MQ=255)
                                             aaCGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGATCAGCGTCGACAGCTGGGCGACTGc                                                                                                                                <  8:190927/120‑1 (MQ=38)
                                             aaCGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGATCAGCGTCGACAGCTGGGCGACTGc                                                                                                                                >  7:190927/1‑120 (MQ=38)
                                              aCGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCgg                                                                                                  >  1:378307/1‑149 (MQ=255)
                                                gACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGcc                                                                                                                               >  3:350411/1‑118 (MQ=255)
                                                gACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGcc                                                                                                                               <  4:350411/118‑1 (MQ=255)
                                                     aTGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCca                                                                                           <  7:96739/149‑1 (MQ=255)
                                                         tCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCcagcag                                                                                       >  4:205176/1‑149 (MQ=255)
                                                          cccAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCagcagc                                                                                      <  6:48933/149‑1 (MQ=255)
                                                                   tGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGgtg                                                                             >  1:176691/1‑149 (MQ=255)
                                                                     aTCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGAGCGTc                                                                           >  2:145073/1‑149 (MQ=255)
                                                                         ccAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGAGTCGGa                                                                        >  4:40017/1‑148 (MQ=255)
                                                                              ccATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAAcgc                                                                  <  2:133155/149‑1 (MQ=255)
                                                                                aTACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGa                                                                <  2:176691/149‑1 (MQ=255)
                                                                                    acCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGct                                                            >  8:255883/1‑149 (MQ=255)
                                                                                        gggTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCttg                                                        >  5:141120/1‑149 (MQ=255)
                                                                                             gtgATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCg                                                                          <  6:94178/126‑1 (MQ=255)
                                                                                               gATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCg                                                   <  7:81442/147‑1 (MQ=255)
                                                                                               gATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCg                                                   >  8:81442/1‑147 (MQ=255)
                                                                                               gATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGt                                                 >  7:21509/1‑149 (MQ=255)
                                                                                               gATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGt                                                 >  4:247435/1‑149 (MQ=255)
                                                                                                    aCTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTg                                            >  8:338303/1‑149 (MQ=255)
                                                                                                     cTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGc                                           <  3:247435/149‑1 (MQ=255)
                                                                                                           gTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGCCGCTCTTGTTGCGGTTCTTGCcgacga                                     >  1:141706/1‑149 (MQ=255)
                                                                                                            tAAGGCGTGTTGCGCACACGTGTGGCAGTCGCCCAGCTGTCGACGCGGATCGACTGCCGGGGCACTCCCAATGTCAACAGAACGCGGTCATCAAGCAGCCGCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGAt                                    <  8:355027/149‑1 (MQ=255)
                                                                                                                  gtgtTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCgg                              >  1:335827/1‑149 (MQ=255)
                                                                                                                       gcgcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAAcgc                                                                  >  7:162028/1‑108 (MQ=255)
                                                                                                                       gcgcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAAcgc                                                                  <  8:162028/108‑1 (MQ=255)
                                                                                                                       gcgcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGt                         >  1:55129/1‑149 (MQ=255)
                                                                                                                        cgcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTg                        >  5:60434/1‑149 (MQ=255)
                                                                                                                        cgcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGAACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTg                        >  1:289964/1‑149 (MQ=255)
                                                                                                                         gcACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTg                        >  8:68303/1‑148 (MQ=255)
                                                                                                                           acacCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCa                      >  7:305162/1‑148 (MQ=255)
                                                                                                                               cTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGCCTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCg                 >  7:23709/1‑149 (MQ=255)
                                                                                                                               cTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCg                 >  7:246987/1‑149 (MQ=255)
                                                                                                                                 gtgGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCGAGCAGCCCGAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTATTGCCGACGATGCGCCGGTCGTGCAAATCGCt               >  6:89035/1‑149 (MQ=255)
                                                                                                                                  tgGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCt               <  3:205176/148‑1 (MQ=255)
                                                                                                                                  tgGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGAGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCt               <  3:40017/148‑1 (MQ=255)
                                                                                                                                    gCAGTCGATCAGCGTCGACAGCTGGGCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGcc            >  4:336086/1‑149 (MQ=255)
                                                                                                                                       gTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGa          >  5:367708/1‑148 (MQ=255)
                                                                                                                                        tCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGaaa        <  7:338303/149‑1 (MQ=255)
                                                                                                                                         cGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAAt       <  8:246987/149‑1 (MQ=255)
                                                                                                                                           cccAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAg     >  5:256534/1‑149 (MQ=255)
                                                                                                                                            ccAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAgg    >  4:177972/1‑149 (MQ=255)
                                                                                                                                              aGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCTCGCGGTCATCCAGCAGCCCCATGGTGTCGGACAACGCGACGCTCTTGTTGCTGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAGGTa  <  5:44448/149‑1 (MQ=255)
                                                                                                                                              aGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAGGTa  >  3:380693/1‑149 (MQ=255)
                                                                                                                                              aGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAGGTa  <  2:7364/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTGCCGTTGGCCCCTGGGCCAGGCCCTCGATGCGGTTGGCATAGAACGACACATGCTCCCAGGGTTTGATCACCAGGCCATACACCGGGGTTGTGATCGACTCATCGTAAGGCGTGTTGCGCACACCTGTGGCAGTCGCCCAGCTGTCGACGCTGATCGACTGCCGGCGCACGCCCACGGTCAACAGCACGCGGTCATCCAGCAGCCCCAGGGTGTCGGACAACGCGACGCTCTTGTTGCGGTTCTTGCCGACGATGCGCGGGTCGTGCAAATCGCTGCCGAAATAGGTA  >  NC_002947/3470709‑3470999

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: