Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,426,616 T→C 6.1% E93G (GAG→GGG)  PP_2127 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,426,6160TC 6.1% 188.3 / 2.6 66E93G (GAG→GGG) PP_2127sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CGTCGCAGTGCCCCTGCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGT  >  NC_002947/2426470‑2426762
                                                                                                                                                  |                                                                                                                                                   
cGTCGCAGTGCCCCTGCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCg                                                                                                                                                   <  1:218563/149‑1 (MQ=255)
               gCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCggg                                                                                                                                   >  4:163009/1‑149 (MQ=255)
                cAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCgg                                                                                                                                    <  2:189275/147‑1 (MQ=255)
                      gtcgGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCCGAGCT‑GGCCGAGCGCat                                                                                                                              <  1:25376/147‑2 (MQ=255)
                          ggCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGAt                                                                                                                        <  3:278329/149‑1 (MQ=255)
                            cATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATc                                                                                                                       >  5:337541/1‑148 (MQ=255)
                                cgcCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCt                                                                                                                  <  4:71322/149‑1 (MQ=255)
                                cgcCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCt                                                                                                                  <  3:163009/149‑1 (MQ=255)
                                       gCAGCCACTGCAGTTGCAGCTTGGCAGCTGTGAGCATTTGACCGATCTCGTCTTGCTGTTCCCGGCCAATATGCTGGCGCTCGTCTTCTTGCACCTTGATCATGCGCTCGTCC‑ATCTCGGCCGGTCTCAGGTGTATCGACCTGTCgcc                                                                                                            <  6:198779/148‑1 (MQ=255)
                                           ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                       <  2:367125/149‑1 (MQ=255)
                                           ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                       >  5:42782/1‑149 (MQ=255)
                                           ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGAGCTCGTCTTCTTGCACCTAGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggc                                                                                                       >  7:149955/1‑149 (MQ=255)
                                            cACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                      >  1:194770/1‑149 (MQ=255)
                                            cACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                      >  3:68572/1‑149 (MQ=255)
                                              cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggccca                                                                                                    <  2:378388/149‑1 (MQ=255)
                                              cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcc                                                                                                      >  2:107515/1‑147 (MQ=255)
                                               tGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCgcccaggcccag                                                                                                   >  5:257582/1‑149 (MQ=255)
                                                 cAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGtt                                                                                                 <  3:20887/149‑1 (MQ=255)
                                                  aGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTg                                                                                                >  7:322712/1‑149 (MQ=255)
                                                      gCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTAGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTac                                                                                             >  1:24631/1‑148 (MQ=255)
                                                       cAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCAGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaa                                                                                           >  4:295267/1‑149 (MQ=255)
                                                         gCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaaca                                                                                         <  7:164486/149‑1 (MQ=255)
                                                         gCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaaca                                                                                         <  1:107515/149‑1 (MQ=255)
                                                                aGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCa                                                                                  <  7:309798/149‑1 (MQ=255)
                                                                  cGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGc                                                                                >  1:260243/1‑149 (MQ=255)
                                                                    gTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGt                                                                              >  4:104465/1‑149 (MQ=255)
                                                                     tGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTg                                                                             <  5:136156/149‑1 (MQ=255)
                                                                       aGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                           >  3:359714/1‑149 (MQ=255)
                                                                         cAGTTGCCCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGa                                                                         <  8:346873/149‑1 (MQ=255)
                                                                         cAGGTGACCGATGGCGTCGTGCAGTTCCCGGCCCAGATGCGGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGa                                                                         <  4:27256/149‑1 (MQ=255)
                                                                          aGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                        <  4:306264/149‑1 (MQ=255)
                                                                          aGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                        >  8:256664/1‑149 (MQ=255)
                                                                            ttGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAcc                                                                       >  4:14147/1‑148 (MQ=255)
                                                                              gACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCaa                                                                     <  6:257582/148‑1 (MQ=255)
                                                                                  gATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTCCAACCGCACCCCGCGGGGCGCCCAACTgcc                                                                >  6:258495/1‑148 (MQ=255)
                                                                                  gATGTCGTCGTGCAGTTCCCCGCCCAGATGCTGGCGCTCGTCTTCTTTCACC‑TGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTgcag                                                               <  6:249418/149‑1 (MQ=255)
                                                                                   aTGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTgcag                                                               <  8:322712/149‑1 (MQ=255)
                                                                                      tcgtcgTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCAGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCAc                                                            <  3:295267/149‑1 (MQ=255)
                                                                                         tcgtGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCCGCACcag                                                         >  2:189611/1‑149 (MQ=255)
                                                                                                    ccGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGcc                                               <  4:359714/148‑1 (MQ=255)
                                                                                                     cGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCATCGTGGCGACCAACTGCAGCACCAGCAGTCTGGCCgg                                             <  6:326532/149‑1 (MQ=255)
                                                                                                     cGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCgg                                             <  7:256664/149‑1 (MQ=255)
                                                                                                          cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTg                                        >  8:122311/1‑149 (MQ=255)
                                                                                                                   gcgcTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                 <  7:41733/147‑1 (MQ=255)
                                                                                                                   gcgcTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTg                                <  6:257028/148‑1 (MQ=255)
                                                                                                                         gTCTTCTTGCACCTTGCGCATGCGCGCGGCC‑CGCTCGGCCGGGCGCAGGGGGCTCGACCTGGCGCCCAGGCCCAGTTGTACACCAGCACCCCGCGTGGCGACCAACTGCAGCACCCGCAGGCTGGCCGGGATTGCCTGGTTGGCCAgg                          <  5:258495/148‑1 (MQ=255)
                                                                                                                         gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAgg                          >  6:280233/1‑148 (MQ=255)
                                                                                                                         gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                           <  3:14147/147‑1 (MQ=255)
                                                                                                                         gTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                           <  2:260243/147‑1 (MQ=255)
                                                                                                                               ttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgcg                   >  3:133543/1‑149 (MQ=255)
                                                                                                                                 gCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgcg                   <  6:334970/147‑1 (MQ=255)
                                                                                                                                 gCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCCGGATTGCCTGGTTGGCCAGGTACAgcg                   >  5:334970/1‑147 (MQ=255)
                                                                                                                                  cACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGcc                 <  6:337541/148‑1 (MQ=255)
                                                                                                                                   aCCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCAt               >  6:312435/1‑149 (MQ=255)
                                                                                                                                    ccTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCa                            >  6:6359/1‑135 (MQ=255)
                                                                                                                                    ccTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCa                            <  5:6359/135‑1 (MQ=255)
                                                                                                                                     cTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGt             >  2:321490/1‑149 (MQ=255)
                                                                                                                                       tccacctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTg                                    <  8:94005/118‑1 (MQ=38)
                                                                                                                                        gAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGtt            <  6:42782/147‑1 (MQ=255)
                                                                                                                                          acctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGTAGGCTGGCCGGGa                                           >  3:87398/4‑114 (MQ=38)
                                                                                                                                          acctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGTAGGCTGGCCGGGa                                           <  4:87398/111‑1 (MQ=38)
                                                                                                                                            ctGCGCCCGGCCGAGCT‑GGCCGAGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAAt      >  7:271365/2‑149 (MQ=255)
                                                                                                                                            aTGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGcc                 >  8:157681/1‑138 (MQ=255)
                                                                                                                                              gcgcTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAAt      >  6:106742/1‑147 (MQ=255)
                                                                                                                                                 cTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGt  >  2:238715/1‑148 (MQ=255)
                                                                                                                                                 cTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGt  >  4:20804/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTCGCAGTGCCCCTGCAACTCGTCGGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACCAACTGCAGCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGT  >  NC_002947/2426470‑2426762

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