Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,778 A→T 6.6% S113T (TCG→ACG)  PP_2203 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,7780AT 6.6% 155.6 / 2.6 61S113T (TCG→ACG) PP_2203hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/19);  new base T (2/2);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

GGTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGC  >  NC_002947/2508631‑2508921
                                                                                                                                                   |                                                                                                                                               
ggTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAg                                                                                                                                                <  8:104090/149‑1 (MQ=255)
 gTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGc                                                                                                                                               <  7:2014/149‑1 (MQ=255)
  tAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGcc                                                                                                                                              >  8:70297/1‑149 (MQ=255)
              ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGcc                                                                                                                                  >  5:250804/1‑149 (MQ=255)
              ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGcc                                                                                                                                  <  5:67943/149‑1 (MQ=255)
                cACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCg                                                                                                                                <  1:403813/149‑1 (MQ=255)
                 aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCgg                                                                                                                               >  6:315425/1‑149 (MQ=255)
                   cAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCa                                                                                                                             >  8:375716/1‑149 (MQ=255)
                          gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAg                                                                                                                      <  7:375716/149‑1 (MQ=255)
                            tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGc                                                                                                                    >  2:198379/1‑149 (MQ=255)
                            tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGccc                                                                                                                                 <  2:30138/136‑1 (MQ=255)
                            tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGccc                                                                                                                                 >  1:30138/1‑136 (MQ=255)
                              agaAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGccc                                                                                                                  >  5:135453/1‑149 (MQ=255)
                               gaAACAGCGCCCCTCGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCGTGGGTCACCCCGTTGACGTGGTAGGCATCCTTGTGGTCACCCTTCTCCTTGCCGGGCGAGCCCTTGGGCATGCCCGGCACCGCCAg                                                                                                                      <  5:369594/144‑1 (MQ=255)
                               gaAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAg                                                                                                                      >  6:369594/1‑144 (MQ=255)
                                 acccgcgcCCCTTGGCGGAGCATCTGTAGTTGGCCCCCCCGGTGGCCTGGCCATCCGGGGGCCCCCCTATTCCGTGGTTGGCATCCTTGTGGCCCCCCATCTCCATGCCGGGCGCCCCCTTGGGCATGCCCGGCACCGCCAGGCCCTTg                                                                                                               <  2:370189/145‑1 (MQ=255)
                                 aaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCtt                                                                                                                >  6:172870/1‑148 (MQ=255)
                                  aCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCA‑CCCGATGGCCTGGTAGGCATGCGTGTTGTGACCCATATCCATTCCGGGCGAGCCCATGGGCATTCCCGGCACCGCCAGGcccgttag                                                                                                             >  1:224848/1‑144 (MQ=255)
                                         cccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGtt                                                                                                       >  2:83312/1‑149 (MQ=255)
                                          ccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTg                                                                                                      >  2:375456/1‑149 (MQ=255)
                                           cTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGc                                                                                                     >  2:78538/1‑149 (MQ=255)
                                                 tGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTgg                                                                                               >  7:95411/1‑149 (MQ=255)
                                                   aGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGGCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGc                                                                                              <  2:129225/148‑1 (MQ=255)
                                                     cATCAGTAGTTGGCGACCCCGGGGTCCTGGCCACCCGGGGCCCCCCCGATGACCTGGTAGGCACCCTTGGGGTCACCCATCACCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcag                                                                                           <  5:315425/149‑1 (MQ=255)
                                                       tCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcag                                                                                           >  7:236443/1‑147 (MQ=255)
                                                         agtagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGTTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                       <  8:95411/149‑1 (MQ=255)
                                                           tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                             >  1:222911/1‑125 (MQ=255)
                                                                 gCGACCAGGGTGTCCTGGCCATCCTGGGTCACCCCGATGACGTGTTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                  <  8:92702/146‑1 (MQ=255)
                                                                 gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                  >  7:92702/1‑146 (MQ=255)
                                                                  cGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGCGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGc                                                                              <  1:241129/149‑1 (MQ=255)
                                                                   gACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCt                                                                             >  4:179193/1‑149 (MQ=255)
                                                                        cGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCTCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggcc                                                                        <  1:78538/149‑1 (MQ=255)
                                                                          gtgtCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccgg                                                                      >  8:287440/1‑149 (MQ=255)
                                                                          gtgtCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccgg                                                                      <  7:287440/149‑1 (MQ=255)
                                                                               cTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGt                                                                 >  3:282004/1‑149 (MQ=255)
                                                                                  gCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTg                                                                <  1:83312/147‑1 (MQ=255)
                                                                                    cATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCt                                                            >  5:129548/1‑149 (MQ=255)
                                                                                        cTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcga                                                        >  4:120156/1‑149 (MQ=255)
                                                                                          gggTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcg                                                                           <  6:297180/128‑1 (MQ=255)
                                                                                          gggTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcg                                                                           >  5:297180/1‑128 (MQ=255)
                                                                                           ggTCACCCCGATGACCTGGTAGGCATCCATGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgacg                                                      >  1:23328/1‑148 (MQ=255)
                                                                                             tCACCCCGATGACGTGGTAGGCATCATTGTGGTCACCAATCTCCATGCCGGGCTAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAAc                                                   <  6:135453/149‑1 (MQ=255)
                                                                                                   cGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccg                                             >  2:131687/1‑149 (MQ=255)
                                                                                                    gATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccgc                                            >  3:162121/1‑149 (MQ=255)
                                                                                                       gACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGa                                         >  7:333429/1‑149 (MQ=255)
                                                                                                                      cTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                           <  6:129548/148‑1 (MQ=255)
                                                                                                                       ttGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCa                         >  3:154118/1‑149 (MQ=255)
                                                                                                                          ggtagaCCCATCTCCATGCCGGGCGAGCCCATCGGCAAGCCCGGCACCGCGAGGCCCGAGAGCTCGATGCGTTTCGCCAGCAGGCGCACCAGCTCGGCCGGCGCGGGGACTTCGTCGAACAAGCCGCAGATCCCCCCGGGGTGGCacca                      >  6:40506/6‑147 (MQ=255)
                                                                                                                               aCCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCaa                 <  3:120156/149‑1 (MQ=255)
                                                                                                                               aCCCATCACCGGGCGGGGCGCCCCCTGGGGCCGGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCaa                 <  1:198379/149‑1 (MQ=255)
                                                                                                                                          tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccgc                                            >  7:99006/1‑111 (MQ=255)
                                                                                                                                          tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccgc                                            <  8:99006/111‑1 (MQ=255)
                                                                                                                                          tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc       >  2:45934/1‑148 (MQ=255)
                                                                                                                                            ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                    >  3:144146/1‑133 (MQ=255)
                                                                                                                                            ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCAc    >  4:145808/1‑149 (MQ=255)
                                                                                                                                            ccGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCAc    >  4:246884/1‑149 (MQ=255)
                                                                                                                                              gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc              <  2:107793/137‑1 (MQ=255)
                                                                                                                                              gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc              >  1:107793/1‑137 (MQ=255)
                                                                                                                                              gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgc         <  7:92275/142‑1 (MQ=255)
                                                                                                                                              gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgc         >  8:92275/1‑142 (MQ=255)
                                                                                                                                              gggCGAGCCCATGCGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGACGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc  >  6:392272/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGC  >  NC_002947/2508631‑2508921

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