Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,935,907 A→G 5.9% M228T (ATG→ACG)  flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,935,9070AG 5.9% 204.3 / 2.5 68M228T (ATG→ACG) flhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base G (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

CATGTCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGC  >  NC_002947/4935759‑4936044
                                                                                                                                                    |                                                                                                                                           
cATGTCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTgcagca                                                                                                                                             >  2:237500/1‑149 (MQ=255)
 aTGTCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCAt                                                                                                                                            >  2:114424/1‑149 (MQ=255)
    tCTTCAGAGCCTGAAGCACGGGTTACCATGATCGCGGCGGCGGAGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGcc                                                                                                                                         >  7:28163/1‑149 (MQ=255)
    tCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGcc                                                                                                                                         >  5:1890/1‑149 (MQ=255)
    tCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGc                                                                                                                                          <  3:201175/148‑1 (MQ=255)
    tCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGc                                                                                                                                          <  1:271000/148‑1 (MQ=255)
        cAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGa                                                                                                                                       <  8:195209/147‑1 (MQ=255)
             ccTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCAt                                                                                                                                >  1:272848/1‑149 (MQ=255)
                  aGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGCCCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCAc                                                                                                                           >  6:158914/1‑149 (MQ=255)
                     aCGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCTCATCGGCATGctgcagc                                                                                                                        <  4:302591/149‑7 (MQ=255)
                              cATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCCCCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTtga                                                                                                             >  7:33277/1‑149 (MQ=255)
                              cATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTtga                                                                                                             >  6:206230/1‑149 (MQ=255)
                               aTGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTtgat                                                                                                            <  5:206230/149‑1 (MQ=255)
                                tGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTtgatg                                                                                                           <  5:348658/149‑1 (MQ=255)
                                   tCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAAc                                                                                                        <  2:10785/149‑1 (MQ=255)
                                     gcggcggcgGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACa                                                                                                       <  1:237500/148‑1 (MQ=255)
                                        gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                           <  1:43405/141‑1 (MQ=255)
                                        gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                           <  1:37901/141‑1 (MQ=255)
                                        gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                           >  2:43405/1‑141 (MQ=255)
                                        gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                           >  2:37901/1‑141 (MQ=255)
                                          ggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCa                                                                                                 <  7:164621/149‑1 (MQ=255)
                                             ggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGc                                                                                               >  4:342026/1‑148 (MQ=255)
                                              gTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGc                                                                                               >  7:352422/1‑147 (MQ=255)
                                                              aTGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTcg                                                                             >  2:299903/1‑149 (MQ=255)
                                                                  cAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGAt                                                                                    <  6:31433/138‑1 (MQ=255)
                                                                  cAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGAt                                                                                    >  5:31433/1‑138 (MQ=255)
                                                                  cAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTcgccgc                                                                         >  4:4224/1‑149 (MQ=255)
                                                                   aaTTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAACAAGCCGGCGATGGCGTCGCcgcg                                                                        <  4:169224/149‑1 (MQ=255)
                                                                       gcgcCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCAc                                                                     <  3:342026/148‑1 (MQ=255)
                                                                       gcgcCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCAc                                                                     <  3:4224/148‑1 (MQ=255)
                                                                        cgcCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGa                                                                   <  2:347442/149‑1 (MQ=255)
                                                                            accaAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACtt                                                               >  2:246244/1‑149 (MQ=255)
                                                                              caAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGc                                                             >  4:93268/1‑149 (MQ=255)
                                                                               aaaCCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCt                                                            <  5:142473/149‑1 (MQ=255)
                                                                                 aCCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTgg                                                          <  6:1890/149‑1 (MQ=255)
                                                                                  ccGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGc                                                         >  2:142991/1‑149 (MQ=255)
                                                                                   cGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGGAGCATGCCGATCAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCa                                                        <  1:343716/149‑1 (MQ=255)
                                                                                    gTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCAc                                                       >  8:237214/1‑149 (MQ=255)
                                                                                     tCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCAcc                                                      >  7:241734/1‑149 (MQ=255)
                                                                                     tCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCAcc                                                      >  3:147697/1‑149 (MQ=255)
                                                                                       aCCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGt                                                    <  1:114424/149‑1 (MQ=255)
                                                                                        ccGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTc                                                   <  2:159383/149‑1 (MQ=255)
                                                                                         cGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTc                                                   >  8:200698/1‑148 (MQ=255)
                                                                                           aTGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCAt                                                >  7:286934/1‑149 (MQ=255)
                                                                                           aTGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCAt                                                >  1:250834/1‑149 (MQ=255)
                                                                                              gTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGa                                             <  8:241734/149‑1 (MQ=255)
                                                                                               tCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAg                                            >  1:337841/1‑149 (MQ=255)
                                                                                               tCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAg                                            >  2:163726/1‑149 (MQ=255)
                                                                                                    agGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTa                                       <  4:147697/149‑1 (MQ=255)
                                                                                                    agGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTa                                       <  8:352422/149‑1 (MQ=255)
                                                                                                       gCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCTAGGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCg                                         >  1:115924/1‑146 (MQ=255)
                                                                                                       gCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATg                                                                                   >  7:232470/1‑102 (MQ=255)
                                                                                                       gCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATg                                                                                   <  8:232470/102‑1 (MQ=255)
                                                                                                        cATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAAc                                   >  1:263043/1‑149 (MQ=255)
                                                                                                          tAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTc                                 <  2:272848/149‑1 (MQ=255)
                                                                                                             aCCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGc                              >  8:49801/1‑149 (MQ=255)
                                                                                                             aCCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGc                              <  8:286934/149‑1 (MQ=255)
                                                                                                              ccTTGCCGGCCTCGGTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGc                              <  3:38626/148‑1 (MQ=255)
                                                                                                              ccTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGcc                             >  1:14218/1‑149 (MQ=255)
                                                                                                                   ccGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCTTCGAGCCGTAGAACTCGGCCTCCTg                        >  8:231406/1‑149 (MQ=255)
                                                                                                                       ccTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGcc                    >  2:176295/1‑149 (MQ=255)
                                                                                                                           gCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCt                <  8:33277/149‑1 (MQ=255)
                                                                                                                                   aCATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCAcgg        <  7:29304/149‑1 (MQ=255)
                                                                                                                                    cATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCAcggc       >  7:118719/1‑149 (MQ=255)
                                                                                                                                     atatTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCAcggcg      >  8:8735/1‑149 (MQ=255)
                                                                                                                                       atTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGt    >  4:250519/1‑149 (MQ=255)
                                                                                                                                       atTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGt    >  1:59757/1‑149 (MQ=255)
                                                                                                                                        ttGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTg   >  4:35568/1‑149 (MQ=255)
                                                                                                                                          gtgCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGc  <  2:263043/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CATGTCTTCAGAGCCTGAAGCACGGGTGACCATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGC  >  NC_002947/4935759‑4936044

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: