Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,414,212 C→G 8.8% S86T (AGC→ACC) ‡ PP_1236 ← glycine cleavage system transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,414,2120CG 8.8% 158.9 / 6.9 57S86T (AGC→ACC) ‡PP_1236glycine cleavage system transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base G (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCACCG  >  NC_002947/1414067‑1414359
                                                                                                                                                  |                                                                                                                                                   
cAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGc                                                                                                                                                  <  1:270392/149‑1 (MQ=255)
        tgtCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCagg                                                                                                                                          <  1:46705/149‑3 (MQ=255)
            gCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTc                                                                                                                                        >  7:203475/1‑147 (MQ=255)
              aGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGATCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAgg                                                                                                                                    <  5:202683/149‑1 (MQ=255)
              aGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCCGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAgg                                                                                                                                    <  3:67234/149‑1 (MQ=255)
                gTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGt                                                                                                                                  >  3:180460/1‑149 (MQ=255)
                        ttCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGt                                                                                                                          >  2:87047/1‑149 (MQ=255)
                                cGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGa                                                                                                                                               >  2:177253/1‑120 (MQ=255)
                                cGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGa                                                                                                                                               <  1:177253/120‑1 (MQ=255)
                                     gggtggtCAAGGAATAACTGGCACAGCTCGTTGTTGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCGGCGGGCGAACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCaa                                                                                                             <  4:180460/148‑1 (MQ=255)
                                          gtCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTg                                                                                                        >  8:218537/1‑149 (MQ=255)
                                            cAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTgg                                                                                                       >  1:93320/1‑149 (MQ=255)
                                             aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGc                                                                                                      <  5:126050/148‑1 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGTGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAgg                                                                                                   >  3:12244/1‑149 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCt                                                                                                         >  3:35746/1‑143 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCt                                                                                                         <  4:35746/143‑1 (MQ=255)
                                                 aagaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                                 <  1:291274/149‑1 (MQ=255)
                                                       tGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTc                                                                                           <  5:180840/149‑1 (MQ=255)
                                                              gCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGt                                                                                    <  2:2308/149‑1 (MQ=255)
                                                                  gtggatcatGGCCGGGCGTTAACCAGCGCTCACG‑AAGGCCACATCCGCATGGGCGTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCa                                                                                <  2:140246/146‑1 (MQ=255)
                                                                  gTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGc                                                                                  <  1:297065/147‑1 (MQ=255)
                                                                   ttGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAg                                                                               <  1:307735/149‑1 (MQ=255)
                                                                         atgTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTc                                                                         >  6:29322/1‑149 (MQ=255)
                                                                         atgTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTc                                                                         >  6:289572/1‑149 (MQ=255)
                                                                            tCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTcc                                                                        >  1:192686/1‑147 (MQ=255)
                                                                             ccGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAg                                                                     >  6:82648/1‑149 (MQ=255)
                                                                              cGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAg                                                                     <  5:82648/148‑1 (MQ=255)
                                                                                        gCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCAc                                                          >  3:295989/1‑149 (MQ=255)
                                                                                          aGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGAc                                                                                      <  5:327002/119‑1 (MQ=255)
                                                                                          aGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGAc                                                                                      >  6:327002/1‑119 (MQ=255)
                                                                                          aGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCt                                                        <  2:192686/149‑1 (MQ=255)
                                                                                             gcTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTg                                                       >  1:177484/1‑147 (MQ=255)
                                                                                              cTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcag                                                    <  5:289572/149‑1 (MQ=255)
                                                                                               tCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCATGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagc                                                   >  2:290380/1‑149 (MQ=255)
                                                                                               tCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcag                                                    >  7:161154/1‑148 (MQ=255)
                                                                                                       ggCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAgg                                             <  2:93320/148‑1 (MQ=255)
                                                                                                       ggCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAgg                                             <  1:87047/148‑1 (MQ=255)
                                                                                                          cacaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGgcgc                                         >  8:263067/1‑149 (MQ=255)
                                                                                                            caTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTg                                       <  8:203475/149‑1 (MQ=255)
                                                                                                               aCGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtc                                    >  3:41042/1‑149 (MQ=255)
                                                                                                                 cgAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCATGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtctc                                  <  1:290380/148‑1 (MQ=255)
                                                                                                                       gCCGGCGGGCGCACTTCCTGGCCCGCGCTGCGCACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTg                            <  1:167980/149‑1 (MQ=255)
                                                                                                                        ccTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTgg                           <  7:218537/149‑1 (MQ=255)
                                                                                                                         cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTgccgc                                                              >  6:180351/1‑113 (MQ=255)
                                                                                                                         cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTgccgc                                                              <  5:180351/113‑1 (MQ=255)
                                                                                                                         cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggc                          <  8:217112/149‑1 (MQ=255)
                                                                                                                                   tcgacgtGGTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcg                >  2:205375/7‑149 (MQ=255)
                                                                                                                                   aCTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcg                >  3:226821/1‑149 (MQ=255)
                                                                                                                                     ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcgc               >  4:27812/1‑148 (MQ=255)
                                                                                                                                      acgtGGTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcgc               >  5:125148/4‑147 (MQ=255)
                                                                                                                                       cgtGGTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCt                                                        <  3:102415/103‑1 (MQ=255)
                                                                                                                                       cgtGGTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCt                                                        >  4:102415/3‑105 (MQ=255)
                                                                                                                                       ccTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa             <  1:141557/148‑1 (MQ=255)
                                                                                                                                        cTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAgg           >  6:23478/1‑149 (MQ=255)
                                                                                                                                        cTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAgg           >  6:233787/1‑149 (MQ=255)
                                                                                                                                        cTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAg            <  1:208777/148‑1 (MQ=255)
                                                                                                                                              cgcgCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGacca     >  8:32960/1‑149 (MQ=255)
                                                                                                                                               gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcg                         >  2:77454/1‑128 (MQ=255)
                                                                                                                                               gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcg                         <  1:77454/128‑1 (MQ=255)
                                                                                                                                               gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGaccac    >  1:292522/1‑149 (MQ=255)
                                                                                                                                               gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGaccac    >  4:264319/1‑149 (MQ=255)
                                                                                                                                                 gcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCACcg  >  3:215448/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACG‑TAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCACCG  >  NC_002947/1414067‑1414359

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