Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,747,983 A→C 10.7% F175V (TTC→GTC)  PP_2403 ← chemotaxis protein CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,747,9830AC10.7% 146.3 / 9.4 56F175V (TTC→GTC) PP_2403chemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base A (27/23);  new base C (4/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CTGGGCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGG  >  NC_002947/2747838‑2748124
                                                                                                                                                 |                                                                                                                                             
cTGGGCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAacac                                                                                                                                            >  4:287443/1‑149 (MQ=255)
   ggCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCg                                                                                                                                         <  1:151833/149‑1 (MQ=255)
   ggCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCg                                                                                                                                         <  2:149527/149‑1 (MQ=255)
     cTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCgca                                                                                                                                       <  3:327702/149‑1 (MQ=255)
      tCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAaca                                                                                                                                             <  8:156972/142‑1 (MQ=255)
      tCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAaca                                                                                                                                             >  7:156972/1‑142 (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                                                                                                                      <  3:324339/148‑1 (MQ=255)
                       cAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGgcg                                                                                                                      >  4:324339/1‑148 (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                                                                                                  <  3:287443/149‑1 (MQ=255)
                                          cTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGt                                                                                                  <  7:241147/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                                                >  5:353119/1‑149 (MQ=255)
                                                                                            ggTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGAc                                                >  4:329261/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        <  2:275306/149‑1 (MQ=255)
                                                                                                                                    gggTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGtt        <  1:379871/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                                            >  3:370951/1‑104 (MQ=255)
                                                                                                                                             aGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGaa                                            <  4:370951/104‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CTGGGCTCGCGGGCATCGTCACCCAGGCGCTGGCGCAGGCGGCTGACCAGCAAGTCGATGGAGCGGTCGAAGATGTCGGCTTCGCGGCCCTGGGTCAGGTTAAGCAGCTGTTCGCGGCTGAGCACCCGCTGCGGGTGGTCGAGGAACACTCGCAGCAGGCGGTACTCGGCGCCGCTCAGGGCCACCAGGGTGCCTTCGCTGTCGAGCAGGTGGCGCGCGGTGGTGTCCAGGCGCCATTGACCGAAACCCAGCAGGCGGCTGCTTTCACTGATGGTCAGGTTGGGCGG  >  NC_002947/2747838‑2748124

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