Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,935,926 C→T 5.9% G222S (GGT→AGT)  flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,935,9260CT 5.9% 204.2 / 2.3 67G222S (GGT→AGT) flhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGG  >  NC_002947/4935789‑4936072
                                                                                                                                         |                                                                                                                                                    
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                                                                                                                                       >  2:43405/1‑141 (MQ=255)
          gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                                                       >  2:37901/1‑141 (MQ=255)
          gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                                                       <  1:37901/141‑1 (MQ=255)
          gcggcggTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCGTGCCGATGAGCATGCCACTGATCAGGTtgatg                                                                                                                                       <  1:43405/141‑1 (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                                                                                                                >  5:31433/1‑138 (MQ=255)
                                    cAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGAt                                                                                                                <  6:31433/138‑1 (MQ=255)
                                    cAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTcgccgc                                                                                                     >  4:4224/1‑149 (MQ=255)
                                     aaTTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAACAAGCCGGCGATGGCGTCGCcgcg                                                                                                    <  4:169224/149‑1 (MQ=255)
                                         gcgcCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCAc                                                                                                 <  3:4224/148‑1 (MQ=255)
                                         gcgcCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCAc                                                                                                 <  3:342026/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                                                                        >  2:163726/1‑149 (MQ=255)
                                                                 tCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAg                                                                        >  1:337841/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                                >  1:59757/1‑149 (MQ=255)
                                                                                                         atTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGt                                >  4:250519/1‑149 (MQ=255)
                                                                                                          ttGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTg                               >  4:35568/1‑149 (MQ=255)
                                                                                                            gtgCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGc                              <  2:263043/148‑1 (MQ=255)
                                                                                                                          cGATGAGCATGCGACCGCTCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCACGGCACGGCGGTGCTTGgcctgggcctgg               <  8:28163/149‑1 (MQ=255)
                                                                                                                           gATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGt              <  2:115924/149‑1 (MQ=255)
                                                                                                                           gATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGt              <  3:35568/149‑1 (MQ=255)
                                                                                                                            aTGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTc             <  5:158914/149‑1 (MQ=255)
                                                                                                                               aGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGa           <  1:163726/148‑1 (MQ=255)
                                                                                                                                  aTGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCGTGGCCTGGGCCTGGTCGATCAg       <  7:237214/149‑1 (MQ=255)
                                                                                                                                      cACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCATGGCCTGGGCCTGGTCGATCAGGCCg   >  5:374163/1‑149 (MQ=255)
                                                                                                                                       aCCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGTCCTGGTCGATCAGGCCgg  <  7:200698/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CATGATCGCGGCGGCGGTGGACAGCAGCAGTGATGGCAATTGCGCCACCAAACCGTCACCGATGGTCAGCAGGGCATAGACCTTGCCGGCCTCGCTGAACGACATATTGTGCTGCAGCATGCCGATGAGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGG  >  NC_002947/4935789‑4936072

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: