Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,233,690 T→C 5.9% V293A (GTG→GCG)  GbpR → galactose‑binding protein regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,233,6900TC 5.9% 198.5 / 2.5 68V293A (GTG→GCG) GbpRgalactose‑binding protein regulator
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/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 = 8.62e-01

CGCTTCGATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCAACTGGATCCGCTACAA  >  NC_002947/3233544‑3233837
                                                                                                                                                  |                                                                                                                                                    
cGCTTCGATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGc                                                                                                                                                    <  1:306503/149‑1 (MQ=255)
cGCTTCGATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGc                                                                                                                                                    >  6:345854/1‑149 (MQ=255)
     cGATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAACGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGAcc                                                                                                                                               >  3:217379/1‑149 (MQ=255)
      gATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGAcc                                                                                                                                               <  6:179761/148‑1 (MQ=255)
                  cACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTcgcg                                                                                                                                   <  3:88640/147‑1 (MQ=255)
                  cACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTcgcg                                                                                                                                   >  4:88640/1‑147 (MQ=255)
                        cATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACt                                                                                                                           >  3:148785/1‑149 (MQ=255)
                        cATGCTCGAAGCGTCCGAAACGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACt                                                                                                                           <  4:217379/149‑1 (MQ=255)
                             tCGTAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGGCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGtt                                                                                                                      <  4:350447/149‑1 (MQ=255)
                             tCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGtt                                                                                                                      >  2:213077/1‑149 (MQ=255)
                             tCGAAGCGGCCGAAATGATTGCGGTGGGGCGCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGGCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTgc                                                                                                                                           <  5:72248/129‑1 (MQ=255)
                               gAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGtt                                                                                                                      <  1:213077/147‑1 (MQ=255)
                                   cGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCtt                                                                                                                <  1:334836/149‑1 (MQ=255)
                                      ccGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGa                                                                                                              >  1:327573/1‑148 (MQ=255)
                                         aaaTGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATAcc                                                                                                          >  6:261985/1‑149 (MQ=255)
                                              aTTGCGGTGGTGCCCAATAACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTtgct                                                                                                     <  6:238747/149‑1 (MQ=255)
                                                  cGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCgg                                                                                                                                  <  6:27569/116‑1 (MQ=255)
                                                  cGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCgg                                                                                                                                  >  5:27569/1‑116 (MQ=255)
                                                           ccAATGACGTGGCCCAGCATTATGACCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTCGCGATGGCCAACCTCGGCGTGCTGACCCTCCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCTTAAACTGCgg                                                                                        >  4:290952/1‑149 (MQ=255)
                                                            cAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCggg                                                                                       >  5:80534/1‑149 (MQ=255)
                                                                    tGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGcactctta                                                                               >  6:87654/1‑144 (MQ=255)
                                                                       cccAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACtt                                                                            >  6:297981/1‑149 (MQ=255)
                                                                             aTTATGCCCGCTATGGCATGGCCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGc                                                                                                                                                    >  5:369151/1‑72 (MQ=255)
                                                                             aTTATGCCCGCTATGGCATGGCCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGc                                                                                                                                                    <  6:369151/72‑1 (MQ=255)
                                                                              ttATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGTCACAATAACTTTCACCTg                                                                     <  5:345854/149‑1 (MQ=255)
                                                                                 tGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCa                                                                   <  3:290952/148‑1 (MQ=255)
                                                                                 tGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCAt                                                                  >  2:50029/1‑149 (MQ=255)
                                                                                   cctgCTATGGCATGGTCGCCGTCCTGCCGGTGTAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGa                                                                <  8:163309/146‑1 (MQ=255)
                                                                                   cccGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGa                                                                <  5:297981/149‑1 (MQ=255)
                                                                                    ccGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAg                                                               >  2:345666/1‑149 (MQ=255)
                                                                                      gCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGc                                                             >  4:230021/1‑149 (MQ=255)
                                                                                          tGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGACCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTg                                                         >  6:3450/1‑149 (MQ=255)
                                                                                            gCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGccc                                                                                                                               <  5:39557/77‑1 (MQ=255)
                                                                                            gCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGccc                                                                                                                               >  6:39557/1‑77 (MQ=255)
                                                                                             cATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGt                                                      >  3:109373/1‑149 (MQ=255)
                                                                                             cATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGt                                                      <  1:345666/149‑1 (MQ=255)
                                                                                              aTGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGACTCGCGGCCCAACTCGGTTGCCCTTGATa                                                                                                            <  1:27389/95‑1 (MQ=255)
                                                                                              aTGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATa                                                                                                            >  2:27389/1‑94 (MQ=255)
                                                                                               tGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCt                                                    <  3:230021/149‑1 (MQ=255)
                                                                                                 gTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTa                                                   <  6:80534/148‑1 (MQ=255)
                                                                                                 gTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTAt                                                  <  2:327573/149‑1 (MQ=255)
                                                                                                   cGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgc                                                >  5:72160/1‑149 (MQ=255)
                                                                                                   cGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGCGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgc                                                >  7:41372/1‑149 (MQ=255)
                                                                                                   cGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGCGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgc                                                >  7:35790/1‑149 (MQ=255)
                                                                                                    gCCGTCCTGCCGGTGGAGTTGCGACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgcc                                               >  7:225760/1‑149 (MQ=255)
                                                                                                    gCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgcc                                               >  2:89953/1‑149 (MQ=255)
                                                                                                    gCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATgc                                                <  4:109373/148‑1 (MQ=255)
                                                                                                       gTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGa                                                                >  7:36054/1‑129 (MQ=255)
                                                                                                       gTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGa                                                                <  8:36054/129‑1 (MQ=255)
                                                                                                          cTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGa                                         >  2:16266/1‑149 (MQ=255)
                                                                                                             ccGGTCGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACAt                                      >  4:165933/1‑149 (MQ=255)
                                                                                                                 tGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGAcc                                  <  1:89953/149‑1 (MQ=255)
                                                                                                                      ttGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCt                                                                      >  2:172769/1‑108 (MQ=255)
                                                                                                                      ttGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCt                                                                      <  1:172769/108‑1 (MQ=255)
                                                                                                                         ccACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGccgccg                                           >  8:132983/1‑132 (MQ=255)
                                                                                                                         ccACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGccgccg                                           <  7:132983/132‑1 (MQ=255)
                                                                                                                           aCTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGAc                                        >  4:344929/1‑133 (MQ=255)
                                                                                                                           aCTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGAc                                        <  3:344929/133‑1 (MQ=255)
                                                                                                                           aCTGGAGATGGCCAACCTCGGCGTGCTGAACCTGCGCTCGCGGCCCAACTAGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGAATGAGGCCCTGCGTCTATGCCGCCGGACATGACGCATTCCAGCt                        >  5:269384/1‑149 (MQ=255)
                                                                                                                               gCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCttt                      <  4:148785/147‑1 (MQ=255)
                                                                                                                               gCGATGGCCAACCTCGGCGCGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCttt                      <  8:35790/147‑1 (MQ=255)
                                                                                                                               gCGATGGCCAACCTCGGCGCGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCttt                      <  8:41372/147‑1 (MQ=255)
                                                                                                                                cGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGc                   <  7:268327/149‑1 (MQ=255)
                                                                                                                                 gATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGc                   >  1:20179/1‑148 (MQ=255)
                                                                                                                                  aTGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCaa                 <  1:50029/149‑1 (MQ=255)
                                                                                                                                   tGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCAAc                <  2:20179/149‑1 (MQ=255)
                                                                                                                                     gCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCAACt               <  3:165933/148‑1 (MQ=255)
                                                                                                                                                 gTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCAACTGGATCCGCTACaa  <  5:261985/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTTCGATCCTGGCGACCACTGCCATGCTCGAAGCGTCCGAAATGATTGCGGTGGTGCCCAATGACGTGGCCCAGCATTATGCCCGCTATGGCATGGTCGCCGTCCTGCCGGTGGAGTTGCCACTGGCGATGGCCAACCTCGGCGTGCTGACCCTGCGCTCGCGGCCCAACTCGGTTGCCCTTGATACCTTGCTGCGTTACCTGCGGGCACAATAACTTTCACCTGCATGAGGCCCTGCGTCTATGCCGCCGGACATGACCCATTCCAGCTTTTGCAACTGGATCCGCTACAA  >  NC_002947/3233544‑3233837

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: