Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,980,807 A→G 6.0% I382T (ATC→ACC)  flgE ← flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,980,8070AG 6.0% 192.2 / 2.1 66I382T (ATC→ACC) flgEflagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base G (2/2);  total (29/37)
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.54e-01

TCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACAGGTAGCCGTTCTG  >  NC_002947/4980666‑4980955
                                                                                                                                             |                                                                                                                                                    
tCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCgg                                                                                                                                               >  5:102179/1‑149 (MQ=255)
tCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCAGAACGCCGATCACCg                                                                                                                                                <  1:358631/148‑1 (MQ=255)
 cAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCgg                                                                                                                                               <  8:287835/148‑1 (MQ=255)
    tGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACg                                                                                                                                           >  2:188557/1‑149 (MQ=255)
      gAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGc                                                                                                                                          <  4:240674/148‑1 (MQ=255)
       aaaTGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGGTGATCGGcgctccg                                                                                                                                        <  1:193202/149‑7 (MQ=255)
            gTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGATCGGCGCTCCGGACACCGGTAcc                                                                                                                                           >  1:65073/1‑140 (MQ=255)
            gTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGATCGGCGCTCCGGACACCGGTAcc                                                                                                                                           <  2:65073/141‑2 (MQ=255)
              cTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCa                                                                                                                                  >  7:23620/1‑148 (MQ=255)
                  gCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGa                                                                                                                             <  1:188557/149‑1 (MQ=255)
                       ggCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGActct                                                                                                                         <  2:353656/148‑1 (MQ=255)
                         ccTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTc                                                                                                                      <  7:11629/149‑1 (MQ=255)
                         ccTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTc                                                                                                                      <  3:233197/149‑1 (MQ=255)
                              tAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGtt                                                                                                                 <  8:206757/149‑1 (MQ=255)
                              tAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGtt                                                                                                                 <  6:30977/149‑1 (MQ=255)
                                gTTGCTCTGCGCCTTGATAAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTgg                                                                                                               >  6:33759/1‑149 (MQ=255)
                                   gCTCTGCGCCTTGATCAGGTTTACCAGCTCGCAGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCCCCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTAc                                                                                                            >  8:317220/1‑149 (MQ=255)
                                     tctGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATa                                                                                                                                >  5:304830/1‑127 (MQ=255)
                                     tctGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATa                                                                                                                                <  6:304830/127‑1 (MQ=255)
                                      ctGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGc                                                                                                         >  7:196594/1‑149 (MQ=255)
                                         cgcCTTGATCAGGTTTACCAGCTCGCCGGTCAGGGCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCaa                                                                                                      <  8:164611/149‑1 (MQ=255)
                                          gcCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGa                                                                                                                             >  5:101097/1‑125 (MQ=255)
                                          gcCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGa                                                                                                                             <  6:101097/125‑1 (MQ=255)
                                          gcCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAAt                                                                                                     >  6:143783/1‑149 (MQ=255)
                                           ccTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATc                                                                                                    <  1:314518/149‑1 (MQ=255)
                                            cTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCg                                                                                                   >  6:44692/1‑149 (MQ=255)
                                             ttGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCgg                                                                                                  >  1:14718/1‑149 (MQ=255)
                                               gATCAGGTTTACCAGCTCGCCGGTCCGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGc                                                                                                 <  8:196594/148‑1 (MQ=255)
                                               gATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGc                                                                                                 <  7:317220/148‑1 (MQ=255)
                                                 tCAGGTTTACCAGCTCGCCGGTCAGGGCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTc                                                                                              <  1:105413/149‑1 (MQ=255)
                                                     gTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGc                                                                                          >  8:358508/1‑149 (MQ=255)
                                                      tttACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGcc                                                                                         <  7:358508/149‑1 (MQ=255)
                                                          ccAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTg                                                                                      >  5:175521/1‑148 (MQ=255)
                                                           cAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTg                                                                                      >  6:299267/1‑147 (MQ=255)
                                                           cAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTg                                                                                      >  4:226208/1‑147 (MQ=255)
                                                              cccGCCGGTCAGGGCGACGTTGGAGTCCTCCAGGGCAGCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATGACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGt                                                                                 <  7:306676/147‑1 (MQ=255)
                                                                    ggTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGa                                                                                                                                      <  5:252488/90‑1 (MQ=255)
                                                                    ggTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGa                                                                                                                                      >  6:252488/1‑90 (MQ=255)
                                                                              aCGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAAt                                                                 >  3:292959/1‑149 (MQ=255)
                                                                                    gAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCAc                                                           <  5:33759/149‑1 (MQ=255)
                                                                                         ctccAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGAc                                                                                                                            >  7:68446/1‑79 (MQ=255)
                                                                                         ctccAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGAc                                                                                                                            <  8:68446/79‑1 (MQ=255)
                                                                                         ctccAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGc                                                      >  7:226074/1‑149 (MQ=255)
                                                                                          tccAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGcc                                                     >  8:125663/1‑149 (MQ=255)
                                                                                                 caccacCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCAcc                                              <  5:44692/149‑1 (MQ=255)
                                                                                                 caccacCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCAcc                                              <  2:315958/149‑1 (MQ=255)
                                                                                                   ccaccGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGTTCACCtt                                            <  3:226208/149‑1 (MQ=255)
                                                                                                     accGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGt                                          <  6:149326/149‑1 (MQ=255)
                                                                                                          aaaTCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTcc                                        <  4:362938/146‑1 (MQ=255)
                                                                                                          aaaTCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTcc                                        >  3:362938/1‑146 (MQ=255)
                                                                                                            aTCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGcc                                   >  5:368333/1‑149 (MQ=255)
                                                                                                             tCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGAt                                                  <  3:147187/133‑1 (MQ=255)
                                                                                                             tCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGAt                                                  >  4:147187/1‑133 (MQ=255)
                                                                                                                 gCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTg                               <  2:14718/148‑1 (MQ=255)
                                                                                                                 gCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTg                               <  8:226074/148‑1 (MQ=255)
                                                                                                                    gAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTg                            >  7:65554/1‑148 (MQ=255)
                                                                                                                      gggTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAg                         <  5:143783/149‑1 (MQ=255)
                                                                                                                       cgtgccGGTGATCGGCGCTCCGGACACCGGTACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGc                        >  1:133822/5‑149 (MQ=255)
                                                                                                                          aCCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTgg                     >  7:31204/1‑149 (MQ=255)
                                                                                                                           ccGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGc                    <  7:125663/149‑1 (MQ=255)
                                                                                                                            cGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGt                             >  4:19189/1‑139 (MQ=255)
                                                                                                                            cGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGt                             <  3:19189/139‑1 (MQ=255)
                                                                                                                             ggTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCa                   >  6:147549/1‑148 (MQ=255)
                                                                                                                               tgtCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAAc                <  6:175521/149‑1 (MQ=255)
                                                                                                                                 tCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACa               <  8:23620/148‑1 (MQ=255)
                                                                                                                                 tCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACAg              <  4:292959/149‑1 (MQ=255)
                                                                                                                                           ctaTCTCCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTTGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCTATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACAGGTAGCCGtt     <  3:205837/146‑1 (MQ=255)
                                                                                                                                             aTCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACAGGTAGCCGTTCTg  <  5:299267/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCCCGAGGCATAGGACTCTTTCCAGTTGGTACCGCCAATCGGCGTCAGGCCCTGCAGGTTGGCGAAGTTGGCAATCGCCACCTGGCCGATCACCTTGTCCTTGCCGTTGGTGAAGCTGGCAAACAGGTAGCCGTTCTG  >  NC_002947/4980666‑4980955

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