Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,725,646 T→G 9.5% K440N (AAA→AAC)  opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,725,6460TG 9.5% 175.9 / 7.7 62K440N (AAA→AAC) opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base G (3/3);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTA  >  NC_002947/5725498‑5725774
                                                                                                                                                     |                                                                                                                                
ggTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCt                                                                                                                                  >  7:158844/1‑149 (MQ=255)
ggTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCt                                                                                                                                  <  6:172743/149‑1 (MQ=255)
  tCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCttt                                                                                                                                <  6:25988/149‑1 (MQ=255)
       cTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                           <  8:158844/149‑1 (MQ=255)
       cTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                           >  3:196959/1‑149 (MQ=255)
            gCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAAc                                                                                                                      >  8:88755/1‑149 (MQ=255)
              tCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTACCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACgg                                                                                                                    >  5:101561/1‑149 (MQ=255)
                gTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACGCACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGgcg                                                                                                                  >  8:253989/1‑149 (MQ=255)
                 tCAAGCACGTTTGGCGGCAGTGCTTCGTATCTGACCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCTTCAGCGCAGCTTCGA‑CGAAGTCGTTGGAAAGGATTGCACCGGCGAACGCGCATTTTACCGTCAACGgcgc                                                                                                                 <  2:88322/149‑1 (MQ=255)
                   aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGc                                                                                                                              >  6:86161/1‑134 (MQ=255)
                   aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGc                                                                                                                              <  5:86161/134‑1 (MQ=255)
                   aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGgcgc                                                                                                                 >  7:151161/1‑147 (MQ=255)
                     gCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                             >  8:56766/1‑149 (MQ=255)
                       aGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGgc                                                                                                           <  1:174370/149‑1 (MQ=255)
                          tttGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCa                                                                                                        <  7:88755/149‑1 (MQ=255)
                            tGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAg                                                                                                       <  7:56766/148‑1 (MQ=255)
                                  cAGTTCTTCGTAGCTACCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGa                                                                                                <  6:101561/149‑1 (MQ=255)
                                   aGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGAt                                                                                               <  8:151161/149‑1 (MQ=255)
                                        ttcTTAGATGCACTGCAAGTCGTAGGCAATCTACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgcgt                                                                                                               <  6:94865/128‑3 (MQ=255)
                                        ttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACg                                                                                          >  8:238936/1‑149 (MQ=255)
                                        ttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgcgt                                                                                                               >  5:94865/1‑126 (MQ=255)
                                           gTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGc                                                                                        <  6:4391/148‑1 (MQ=255)
                                           gTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCt                                                                                       <  1:183004/149‑1 (MQ=255)
                                                 gCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGGATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGc                                                                                 >  1:191652/1‑149 (MQ=255)
                                                 gCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                             >  5:126026/1‑121 (MQ=255)
                                                 gCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                             <  6:126026/121‑1 (MQ=255)
                                                  cccGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGc                                                                                 <  2:219158/148‑1 (MQ=255)
                                                           tCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGgatga                                                                           <  7:65141/145‑1 (MQ=255)
                                                           tCGTAGGCAATCCACACACCCCAGCCGGCACGGCCCATCAGCGCAGCTTCGA‑CGAACACCTGCGAAAGGATGGCACCGGCGAACCCGCCTTTGCCCGCCAACCGCGCCAGCGCGCAGAGCTCCATGAACGGCTTCATGCGgatgatac                                                                        >  8:65141/1‑146 (MQ=255)
                                                             gTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATcgc                                                                       <  4:39401/147‑1 (MQ=255)
                                                             gTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATcgc                                                                       >  3:39401/1‑147 (MQ=255)
                                                                     tCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                     <  6:254409/141‑1 (MQ=255)
                                                                     tCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                     >  5:254409/1‑141 (MQ=255)
                                                                             ccccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTg                                                     <  5:45067/149‑1 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAAc                                                                                           <  5:166120/110‑1 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAAc                                                                                           >  6:166120/1‑110 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATg                                                                                  <  3:81219/119‑1 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATg                                                                                  >  4:81219/1‑119 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGc                                                    <  1:147797/149‑1 (MQ=255)
                                                                                cAGCCGGCACGGCGCATCAGCGCAGCTTCGAGCGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGc                                                    >  1:28868/1‑148 (MQ=255)
                                                                                           gcgcATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACt                                       <  5:140830/149‑1 (MQ=255)
                                                                                               cacagcgcagcgTCGA‑CGAAGACGTCGGAACCGATGGCACCGGCGAACGCGCCTTTGCCCCGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGG‑CCCAGTAGTGCGACTCACCCAACTGCCAg                                  <  8:114946/147‑1 (MQ=255)
                                                                                               aTCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGc                                     >  5:238550/1‑147 (MQ=255)
                                                                                                             ga‑ctaAGTCGTGTGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCgg                     <  7:241870/149‑1 (MQ=255)
                                                                                                                cgaAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCt                                                                                       >  6:171119/1‑81 (MQ=255)
                                                                                                                cgaAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCt                                                                                       <  5:171119/81‑1 (MQ=255)
                                                                                                                 gaAGTCGGGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTGCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGAGTCACCCAACTGCCAGAAGAGCAGACCGGccg                  >  8:80717/1‑147 (MQ=255)
                                                                                                                       gTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAAc             <  7:201993/148‑1 (MQ=255)
                                                                                                                          ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGt          <  1:79328/148‑1 (MQ=255)
                                                                                                                          ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTg         >  2:62845/1‑149 (MQ=255)
                                                                                                                          ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTg         >  2:25244/1‑149 (MQ=255)
                                                                                                                          ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTg         <  6:238550/149‑1 (MQ=255)
                                                                                                                           gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGa                                                  <  7:254353/107‑1 (MQ=255)
                                                                                                                           gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGa                                                  >  8:254353/1‑107 (MQ=255)
                                                                                                                           gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                           <  3:121053/130‑1 (MQ=255)
                                                                                                                           gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                           >  4:121053/1‑130 (MQ=255)
                                                                                                                             aaGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACACTGGGc      <  7:238936/149‑1 (MQ=255)
                                                                                                                                 aTGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTa  <  2:28868/149‑1 (MQ=255)
                                                                                                                                  tGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTa  >  5:150902/1‑148 (MQ=255)
                                                                                                                                             ctggcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCa                                              <  1:92911/90‑1 (MQ=38)
                                                                                                                                             ctggcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCa                                              >  2:92911/4‑93 (MQ=38)
                                                                                                                                             aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCAcc                                            <  4:18048/95‑1 (MQ=255)
                                                                                                                                             aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCAcc                                            >  3:18048/1‑95 (MQ=255)
                                                                                                                                                gcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCg    >  4:111996/1‑132 (MQ=255)
                                                                                                                                                gcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCg    <  3:111996/132‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                
GGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGA‑CGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTA  >  NC_002947/5725498‑5725774

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