Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,809 T→A 6.8% D695V (GAC→GTC)  PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,684,8090TA 6.8% 171.9 / 2.9 59D695V (GAC→GTC) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base A (2/2);  total (28/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 = 8.86e-01

GCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATG  >  NC_002947/5684666‑5684957
                                                                                                                                               |                                                                                                                                                    
gCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACg                                                                                                                                                 <  4:87936/149‑1 (MQ=255)
 cTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACgg                                                                                                                                                >  5:32589/1‑149 (MQ=255)
       ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTgg                                                                                                                                           <  5:150671/148‑1 (MQ=255)
       ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                                          >  4:68031/1‑149 (MQ=255)
       ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                                          <  3:68031/149‑1 (MQ=255)
       ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTgg                                                                                                                                           <  2:303260/148‑1 (MQ=255)
       ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGc                                                                                                                                          <  2:50534/149‑1 (MQ=255)
        tgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGc                                                                                                                                          >  7:133015/1‑148 (MQ=255)
           cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                      <  6:242692/149‑1 (MQ=255)
           cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCCTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                      <  2:202534/149‑1 (MQ=255)
            aGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                            >  6:63912/1‑142 (MQ=255)
            aGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                            <  5:63912/142‑1 (MQ=255)
             gcgcGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGc                                                                                                                                          <  1:328032/143‑1 (MQ=255)
             gcgcGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGc                                                                                                                                          >  2:328032/1‑143 (MQ=255)
                    aacaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                             >  8:20480/1‑149 (MQ=255)
                    aacaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                             <  8:133015/149‑1 (MQ=255)
                       aacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAg                                                                                                                          <  7:20480/149‑1 (MQ=255)
                         cGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                             >  3:219589/1‑144 (MQ=255)
                         cGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                             <  4:219589/144‑1 (MQ=255)
                            ctctGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCcg                                                                                                                     <  3:139504/149‑1 (MQ=255)
                              ctGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCAGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGAGGTCCACGGCCTGGCCAACCAGCTCGATCAGTACcgcg                                                                                                                   >  4:289036/1‑149 (MQ=255)
                                   tCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                      >  5:128403/1‑125 (MQ=255)
                                   tCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                      <  6:128403/125‑1 (MQ=255)
                                           gATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATgg                                                                                                      >  6:256782/1‑149 (MQ=255)
                                              gTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                            <  3:252349/108‑1 (MQ=255)
                                              gTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                            >  4:252349/1‑108 (MQ=255)
                                                ccATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCt                                                                                                                                 <  3:59514/117‑1 (MQ=255)
                                                ccATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCt                                                                                                                                 >  4:59514/1‑117 (MQ=255)
                                                 catggcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGATCCGCGGCACCGGGGATggcggc                                                                                                  >  5:206065/1‑147 (MQ=255)
                                                      catggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCgg                                                                                                                  <  6:86751/126‑1 (MQ=255)
                                                                gAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                 <  6:32589/149‑1 (MQ=255)
                                                                         gACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCTCCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                 <  6:363359/140‑1 (MQ=255)
                                                                         gACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                 >  5:363359/1‑140 (MQ=255)
                                                                         gACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGcc                                                                         >  6:187167/1‑148 (MQ=255)
                                                                              aCTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGTGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGcc                                                                    >  7:330634/1‑148 (MQ=255)
                                                                               cTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGcc                                                                    >  8:31844/1‑147 (MQ=255)
                                                                                 gCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGcc                                                                    >  1:307204/1‑145 (MQ=255)
                                                                                  cGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGAt                                                               <  2:338348/149‑1 (MQ=255)
                                                                                   ggATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                                              <  5:187167/149‑1 (MQ=255)
                                                                                    gATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATgg                                                             <  5:256782/149‑1 (MQ=255)
                                                                                    gATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATgg                                                             <  7:31844/149‑1 (MQ=255)
                                                                                        aCGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaac                                                          <  3:95052/148‑1 (MQ=255)
                                                                                          gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCTCAGATGGaacaac                                                       >  1:89702/1‑149 (MQ=255)
                                                                                            ttGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGa                                                     >  3:16805/1‑149 (MQ=255)
                                                                                                tgcGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                 <  2:80127/117‑1 (MQ=255)
                                                                                                tgcGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                 >  1:80127/1‑117 (MQ=255)
                                                                                                          aTGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGCTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCa                                       <  6:206065/149‑1 (MQ=255)
                                                                                                                cGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccatt                                 <  4:16805/149‑1 (MQ=255)
                                                                                                                     cTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccatt                                 >  6:251672/1‑144 (MQ=255)
                                                                                                                     cTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccatt                                 <  5:251672/144‑1 (MQ=255)
                                                                                                                         aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCgg                        >  4:127432/1‑149 (MQ=255)
                                                                                                                         aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCg                         <  2:307204/148‑1 (MQ=255)
                                                                                                                              ggttGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAtt                   >  5:158321/4‑149 (MQ=38)
                                                                                                                                   gCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAg              >  8:131340/1‑149 (MQ=255)
                                                                                                                                          ggcggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGt       <  7:131340/149‑1 (MQ=255)
                                                                                                                                          ggcggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGt       >  6:150523/1‑149 (MQ=255)
                                                                                                                                          ggcggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGt       >  2:127380/1‑149 (MQ=255)
                                                                                                                                          ggcggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGt       <  3:127432/149‑1 (MQ=255)
                                                                                                                                            cggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGtt      >  8:143908/1‑148 (MQ=255)
                                                                                                                                               tCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTc                                         <  2:322115/110‑1 (MQ=255)
                                                                                                                                               tCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATg  <  7:103120/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATG  >  NC_002947/5684666‑5684957

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: ‑