Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 497,697 G→T 7.2% L523L (CTG→CTT PP_0409 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947497,6970GT 7.2% 215.6 / 3.6 68L523L (CTG→CTTPP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base T (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGC  >  NC_002947/497553‑497845
                                                                                                                                                |                                                                                                                                                     
tCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGAcctggccctggc                                                                                                                                                     >  2:347366/1‑147 (MQ=255)
  tGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGc                                                                                                                                                 >  2:188868/1‑149 (MQ=255)
   ggAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTTCGCAACCTGGTGCACCCCGAGGACCGTGCTGATAATCACAGGGTGCACGACCTGGCCCTTGCCAGc                                                                                                                                                 >  5:77297/1‑148 (MQ=255)
    gAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCg                                                                                                                                                <  3:199169/148‑1 (MQ=255)
      gTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGAc                                                                                                                                              >  3:96816/1‑148 (MQ=255)
           ccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                        <  7:218081/149‑1 (MQ=255)
           ccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                        >  8:218081/1‑149 (MQ=255)
           ccTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCCTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGAc                                                                                                                                        >  7:5397/1‑149 (MQ=255)
            cTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAgtcctg                                                                                                                                      >  6:261760/1‑144 (MQ=255)
               tatCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTg                                                                                                                                      <  4:164570/147‑1 (MQ=255)
                atCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCATCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTc                                                                                                                                   >  5:184190/1‑149 (MQ=255)
                        gggTAGCGCGGCGTTGGTGGGGTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAGTCCTGGTCGCTGGccag                                                                                                                          <  6:343255/149‑4 (MQ=255)
                               gAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCAt                                                                                                                   >  5:162070/1‑149 (MQ=255)
                                aGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAGTCCTGGTCGCTGGCCAGGgccaggt                                                                                                                  <  5:261760/149‑6 (MQ=37)
                                  gCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCt                                                                                                                >  3:175212/1‑149 (MQ=255)
                                  gCTTTGGGGGGGTAAGCACCCAGCGAGAACCAGCACCCGCAGCTGGGCCTGCGCAACCGGGGGCACCCCGACGACCGTGCGGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAGGACCAGTCCTGGTCGCTGGCCAGGgccaggtcc                                                                                                                <  3:157683/149‑8 (MQ=16)
                                    tttGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCa                                                                                                              >  5:344573/1‑149 (MQ=255)
                                      tggtggGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCAc                                                                                                             <  6:344573/148‑1 (MQ=255)
                                               aTGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAggg                                                                                                    <  1:347366/148‑1 (MQ=255)
                                                   aCCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGGGAAc                                                                                               >  4:302629/1‑149 (MQ=255)
                                                   aCCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGAACAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGacc                                                                                               >  5:268228/1‑147 (MQ=255)
                                                     ccAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAg                                                                                             >  3:106520/1‑149 (MQ=255)
                                                      cAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCAccc                                                                                                           >  2:314577/1‑134 (MQ=255)
                                                      cAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCAccc                                                                                                           <  1:314577/134‑1 (MQ=255)
                                                      cAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAgc                                                                                            <  4:106520/149‑1 (MQ=255)
                                                       aGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAgcg                                                                                           >  8:289900/1‑149 (MQ=255)
                                                             aTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                     >  8:30792/1‑149 (MQ=255)
                                                                cAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTggcctt                                                                                   >  7:53882/1‑144 (MQ=255)
                                                                cAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGa                                                                                  <  3:95670/149‑1 (MQ=255)
                                                                  gCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGata                                                                                <  1:222852/149‑1 (MQ=255)
                                                                    aCCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGatat                                                                               >  8:206507/1‑148 (MQ=255)
                                                                      ccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCt                                                                                         >  1:202903/1‑136 (MQ=255)
                                                                      ccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCt                                                                                         <  2:202903/136‑1 (MQ=255)
                                                                       cGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAg                                                                           <  7:206507/149‑1 (MQ=255)
                                                                       cGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAg                                                                           <  7:289900/149‑1 (MQ=255)
                                                                                 ccTGAGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGccc                                                                 <  6:77297/149‑1 (MQ=255)
                                                                                  cTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCg                                                                >  4:190614/1‑149 (MQ=255)
                                                                                   tGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                     <  5:82246/127‑1 (MQ=255)
                                                                                   tGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                     >  6:82246/1‑127 (MQ=255)
                                                                                     cgcAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGc                                                              <  1:188868/148‑1 (MQ=255)
                                                                                     cgcAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGc                                                              <  8:66990/148‑1 (MQ=255)
                                                                                        aaCCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCCGCGCCCGGCATCt                                                          >  5:244136/1‑149 (MQ=255)
                                                                                          ccTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTgg                                                        <  3:302629/149‑1 (MQ=255)
                                                                                              gTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCaa                                                    >  8:160743/1‑149 (MQ=255)
                                                                                                 cACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACgg                                                 >  2:48690/1‑149 (MQ=255)
                                                                                                   ccccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGcc                                               <  3:190614/149‑1 (MQ=255)
                                                                                                     ccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                     >  5:206391/1‑109 (MQ=255)
                                                                                                     ccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGc                                                                                     <  6:206391/109‑1 (MQ=255)
                                                                                                     ccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGccc                                                                 <  5:160678/129‑1 (MQ=255)
                                                                                                     ccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGccc                                                                 >  6:160678/1‑129 (MQ=255)
                                                                                                         ggACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTg                                                                                        >  8:4504/1‑102 (MQ=255)
                                                                                                         ggACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTg                                                                                        <  7:4504/102‑1 (MQ=255)
                                                                                                            ccGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAgtgt                                      <  4:175212/149‑1 (MQ=255)
                                                                                                            ccGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAgtgt                                      <  7:30792/149‑1 (MQ=255)
                                                                                                               tGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                    >  5:140851/1‑148 (MQ=255)
                                                                                                                  tGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACg                                >  8:197666/1‑149 (MQ=255)
                                                                                                                       aTCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGgtg                             >  7:9963/1‑147 (MQ=255)
                                                                                                                        tCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAgtg                                       <  3:327314/136‑1 (MQ=255)
                                                                                                                        tCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAgtg                                       >  4:327314/1‑136 (MQ=255)
                                                                                                                           caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                                 >  7:162156/1‑91 (MQ=255)
                                                                                                                           caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                                 <  8:162156/91‑1 (MQ=255)
                                                                                                                           caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCt                       <  1:48690/149‑1 (MQ=255)
                                                                                                                             gggTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTgg                     <  3:346348/149‑1 (MQ=255)
                                                                                                                               gTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGACAGGTAGTGTGGGACGGTGTAGTCTGGGa                   >  2:252907/1‑149 (MQ=255)
                                                                                                                                tGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGa                   >  3:122107/1‑148 (MQ=255)
                                                                                                                                  cAGGACCTGGCCCTGGCCAGCGACCAGGACTGTTCCTGGGAGGGGCGCATCCTCACGCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACAt                <  8:5397/149‑1 (MQ=255)
                                                                                                                                  cAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACAt                >  7:5944/1‑149 (MQ=255)
                                                                                                                                    ggACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTTGGACGGTGTAGTCTGGGACATca              <  6:184190/149‑1 (MQ=255)
                                                                                                                                    ggACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATca              >  6:222088/1‑149 (MQ=255)
                                                                                                                                       cctggccctggcCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATcaca            >  4:51833/1‑148 (MQ=255)
                                                                                                                                           gccctggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggg       <  3:51833/149‑1 (MQ=255)
                                                                                                                                            ccctggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGGGTAGTCTGGGACGTCACACAggg       >  7:270556/1‑148 (MQ=255)
                                                                                                                                             cctggccAGCGACCATGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACCGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTCGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGCACATCACACAGGGCa     >  7:181046/1‑149 (MQ=255)
                                                                                                                                                ggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                    >  2:332869/1‑115 (MQ=255)
                                                                                                                                                ggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                    <  1:332869/115‑1 (MQ=255)
                                                                                                                                                ggccAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGc  >  7:216728/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTCCTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGCAAGC  >  NC_002947/497553‑497845

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: