Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,747,987 T→C 10.3% R173R (CGA→CGG)  PP_2403 ← chemotaxis protein CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,747,9870TC10.3% 149.4 / 7.3 58R173R (CGA→CGG) PP_2403chemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (4/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGGCAGCATGCGT  >  NC_002947/2747841‑2748134
                                                                                                                                                  |                                                                                                                                                   
ggCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCg                                                                                                                                                   <  1:151833/149‑1 (MQ=255)
ggCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCg                                                                                                                                                   <  2:149527/149‑1 (MQ=255)
  cTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCgca                                                                                                                                                 <  3:327702/149‑1 (MQ=255)
     gcgGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGAGTGTTCCTCGACCACCCgcagc                                                                                                                                               <  8:240588/148‑1 (MQ=255)
     gcgGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCACCAAGTCGATGGAGCGGTCGAAGATTTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCTCGCCTGAGCACCCGCTGCGGGTGGTAGAGGAACACTCgcag                                                                                                                                                >  2:379871/1‑147 (MQ=255)
       gggCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAgg                                                                                                                                            <  8:7627/149‑1 (MQ=255)
           aTCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCgg                                                                                                                                         >  4:171281/1‑148 (MQ=255)
              gTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACt                                                                                                                                     <  2:59573/149‑1 (MQ=255)
              gTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTAGCGGCCCTGGGTCAGGTTAAGCAGCTGTTAGCGGCTGAGCCCCCCCTGCGGGTGGTAGAGGAACACTCGCAGCAGGCGGTACt                                                                                                                                     >  1:225487/1‑149 (MQ=255)
                  cccAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGgc                                                                                                                                 >  5:177491/1‑149 (MQ=255)
                   ccAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGgcg                                                                                                                                >  5:323980/1‑149 (MQ=255)
                    cAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGgcg                                                                                                                                >  4:324339/1‑148 (MQ=255)
                    cAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGgcg                                                                                                                                <  3:324339/148‑1 (MQ=255)
                       gcgcTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGcgccgc                                                                                                                            >  1:188315/1‑149 (MQ=255)
                         gcTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTc                                                                                                                          <  5:10129/149‑1 (MQ=255)
                           tGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGAGTGTTCCTCGACCACCCgcagc                                                                                                                                               >  3:230666/1‑126 (MQ=255)
                           tGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGAGTGTTCCTCGACCACCCgcagc                                                                                                                                               <  4:230666/126‑1 (MQ=255)
                             gcgcAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAg                                                                                                                        <  2:320041/147‑1 (MQ=255)
                               gcAGGCGGTTGCCCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTAAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGcc                                                                                                                    >  7:151355/1‑149 (MQ=255)
                                       cTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGt                                                                                                            <  7:241147/149‑1 (MQ=255)
                                       cTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGt                                                                                                            <  3:287443/149‑1 (MQ=255)
                                                  gTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGt                                                                                                 <  2:188315/149‑1 (MQ=255)
                                                  gTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGt                                                                                                 >  4:150228/1‑149 (MQ=255)
                                                     gATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTc                                                                                                >  4:70355/1‑147 (MQ=255)
                                                          aGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAgg                                                                                         >  1:275306/1‑149 (MQ=255)
                                                            cGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTg                                                                                       <  6:177491/149‑1 (MQ=255)
                                                               tCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGgcg                                                                                    >  4:245785/1‑149 (MQ=255)
                                                                     aTGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCggt                                                                              >  5:227963/1‑149 (MQ=255)
                                                                        tCGGCTTCGCGGCCCTGGGTCATGTTAAGCAGCTGGTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGGAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCggtggt                                                                           <  3:104376/149‑1 (MQ=255)
                                                                        tCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCggtggt                                                                           >  6:18752/1‑149 (MQ=255)
                                                                          ggCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGgtgt                                                                         <  6:323980/149‑1 (MQ=255)
                                                                                gcgGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAAGACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCAGTCGAGCAGGTGGCTCGCCGTGGTGTCCAGgc                                                                   >  7:168841/1‑149 (MQ=255)
                                                                                     ccTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCAt                                                              >  4:358498/1‑149 (MQ=255)
                                                                                         ggTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGAc                                                          >  4:329261/1‑149 (MQ=255)
                                                                                         ggTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGAc                                                          >  5:353119/1‑149 (MQ=255)
                                                                                              ggTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGaaa                                                     <  6:227963/149‑1 (MQ=255)
                                                                                               gTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAAc                                                    >  1:363421/1‑149 (MQ=255)
                                                                                               gTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGCTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGCTGGTGTCCAGGCGCCATTGACCGAAAc                                                    >  1:311376/1‑149 (MQ=255)
                                                                                                  aaGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAAccc                                                  <  8:214277/148‑1 (MQ=255)
                                                                                                    gcagcTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCcagc                                               <  3:245785/149‑1 (MQ=255)
                                                                                                    gcagcTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCcag                                                >  4:372786/1‑148 (MQ=255)
                                                                                                               cgGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGgctgct                                    <  6:353119/149‑1 (MQ=255)
                                                                                                                         cctgCTGCGAGTGTTCCTCGACCACCCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATg                          >  6:237027/4‑149 (MQ=255)
                                                                                                                         cccGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATg                          >  4:367493/1‑149 (MQ=255)
                                                                                                                            gCTGCGAGTGTTCCTCGACCACCCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGt                        >  3:370229/1‑148 (MQ=255)
                                                                                                                            gCTGCGAGTGTTCCTCGACCACCCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGt                        >  5:44196/1‑148 (MQ=255)
                                                                                                                                cGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAg                     <  7:291992/147‑1 (MQ=255)
                                                                                                                                 gggTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGtt                  <  1:379871/149‑1 (MQ=255)
                                                                                                                                 gggTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGtt                  <  2:275306/149‑1 (MQ=255)
                                                                                                                                       tCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGAt                           >  7:342711/1‑134 (MQ=255)
                                                                                                                                       tCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGggcgg            >  1:3098/1‑149 (MQ=255)
                                                                                                                                       tCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGggcgg            >  1:157809/1‑149 (MQ=255)
                                                                                                                                       tCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGGGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGAt                           <  8:342711/134‑1 (MQ=255)
                                                                                                                                        cGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGggcgg            >  6:136731/1‑148 (MQ=255)
                                                                                                                                          aGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGaa                                                      <  4:370951/104‑1 (MQ=255)
                                                                                                                                          aGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGaa                                                      >  3:370951/1‑104 (MQ=255)
                                                                                                                                               cacTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGGCAGCATGc    <  2:363421/149‑1 (MQ=255)
                                                                                                                                                 cTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGGCAGCATGCgt  <  3:329261/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGGCAGCATGCGT  >  NC_002947/2747841‑2748134

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