Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,720,690 Δ1 bp 6.2% coding (1215/1944 nt) PP_5021 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,720,6900C. 6.2% 206.2 / 6.5 64coding (1215/1944 nt)PP_5021methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base . (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CAGCATGCGCGAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCGGCCGTATCGAA  >  NC_002947/5720543‑5720837
                                                                                                                                                   |                                                                                                                                                   
cAGCATGCGCGAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGAc                                                                                                                                                     >  7:242297/1‑148 (MQ=255)
  gCATGCGCGAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAg                                                                                                                                                  <  8:242297/149‑1 (MQ=255)
  gCATGCGCGAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAg                                                                                                                                                  <  3:56408/149‑1 (MQ=255)
        gcgAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGc                                                                                                                                             >  4:216524/1‑148 (MQ=255)
          gagaTGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGAc                                                                                                                                          >  4:187099/1‑149 (MQ=255)
                  cTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTg                                                                                                                                        <  1:171638/143‑1 (MQ=255)
                  cTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTg                                                                                                                                        >  2:171638/1‑143 (MQ=255)
                        ccTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCAACAATCAGCAGCGCGAGGAAACCGA‑CAGGTTGCGACTGCCATGAACCAGAt                                                                                                                           >  1:185195/1‑149 (MQ=255)
                                 ccTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCa                                                                                                                   >  6:277752/1‑149 (MQ=255)
                                  cTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACAGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCAc                                                                                                                  <  5:264144/149‑1 (MQ=255)
                                      tCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGTCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTg                                                                                                              >  3:197108/1‑149 (MQ=255)
                                              aTTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGt                                                                                                      >  3:126182/1‑149 (MQ=255)
                                              aTTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCAACAATCAGCAGCGCGAGGAAACCGA‑CAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTg                                                                                                     <  2:185195/149‑1 (MQ=255)
                                               ttGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCAcc                                                                                                                 <  1:159570/137‑1 (MQ=255)
                                               ttGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCAcc                                                                                                                 >  2:159570/1‑137 (MQ=255)
                                                 gATCAGGGTGTCGGCCAGCTGTCACAGGTCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCTGATGAGCGCCACCGTGCAGGAAGTgg                                                                                                    <  4:197108/148‑1 (MQ=255)
                                                 gATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTgg                                                                                                    <  2:32870/148‑1 (MQ=255)
                                                  aTCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGca                                                                                                  >  7:270293/1‑149 (MQ=255)
                                                      gggTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAAc                                                                                                                                >  5:140101/1‑115 (MQ=255)
                                                      gggTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAAc                                                                                                                                <  6:140101/115‑1 (MQ=255)
                                                        gtgtCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAAc                                                                                            >  3:293283/1‑149 (MQ=255)
                                                             ggCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGa                                                                                       <  2:318109/149‑1 (MQ=255)
                                                              gCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAg                                                                                      >  2:172838/1‑149 (MQ=255)
                                                                 aGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCaa                                                                                   <  1:165781/149‑1 (MQ=255)
                                                                 aGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGA‑CAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGc                                                                                     >  1:172659/1‑146 (MQ=255)
                                                                 aGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGA‑CAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCaa                                                                                   <  2:172659/148‑1 (MQ=255)
                                                                  ttcgctacaGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAg                                                                                  >  3:30206/7‑149 (MQ=255)
                                                                   cTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGc                                                                                 <  4:95987/149‑1 (MQ=255)
                                                                        acaGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGGCTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTcc                                                                             <  2:14034/148‑1 (MQ=255)
                                                                          aGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTg                                                                          >  3:61115/1‑149 (MQ=255)
                                                                              cgccACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGa                                                                                       <  7:16019/132‑1 (MQ=255)
                                                                              cgccACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGa                                                                                       >  8:16019/1‑132 (MQ=255)
                                                                              cgccACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCcg                                                                      >  8:65174/1‑149 (MQ=255)
                                                                                ccACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCcgcg                                                                    <  1:172838/149‑1 (MQ=255)
                                                                                 cACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGa                                                                   <  5:106393/149‑1 (MQ=255)
                                                                                  aCGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGAc                                                                  <  8:270293/149‑1 (MQ=255)
                                                                                           gcggcCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGa                                                         >  6:267474/1‑149 (MQ=255)
                                                                                            cggcCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGAc                                                        <  3:216524/149‑1 (MQ=255)
                                                                                             ggcCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGAc                                                        >  2:343134/1‑148 (MQ=255)
                                                                                                     gCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGCCCAACGCTGACCAGCAGgcg                                               >  7:86325/1‑149 (MQ=255)
                                                                                                     gCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGgcg                                               >  8:176641/1‑149 (MQ=255)
                                                                                                       gAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCa                                             <  4:30206/149‑1 (MQ=255)
                                                                                                            cacCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgg                                        <  3:187099/149‑1 (MQ=255)
                                                                                                            cacCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgg                                        <  3:194196/149‑1 (MQ=255)
                                                                                                             acCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgg                                        >  2:203559/1‑148 (MQ=255)
                                                                                                             acCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgg                                        <  1:203559/148‑1 (MQ=255)
                                                                                                               cAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCATCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgggt                                      >  4:322681/1‑148 (MQ=255)
                                                                                                               cAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGc                                                                                 <  8:233307/105‑1 (MQ=255)
                                                                                                               cAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGc                                                                                 >  7:233307/1‑105 (MQ=255)
                                                                                                               cAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAggtgggt                                      >  7:46414/1‑148 (MQ=255)
                                                                                                               cAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATCAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACCCCGACCAAGCCTCACTGCCCGCGACCAACGCTGACCAGCATGCGCAcgtcggt                                      >  3:57335/1‑148 (MQ=255)
                                                                                                                aaGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGa                                    <  5:277752/149‑1 (MQ=255)
                                                                                                                 aGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGAc                                   <  1:343134/149‑1 (MQ=255)
                                                                                                                      cgcATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAgg                               <  4:293283/148‑1 (MQ=255)
                                                                                                                          tcaatcaGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCCACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGt                           >  7:259553/1‑148 (MQ=255)
                                                                                                                            aatcaGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGtctc                        <  5:267474/149‑1 (MQ=255)
                                                                                                                            aatcaGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGtctc                        >  7:17735/1‑149 (MQ=255)
                                                                                                                            aatcaGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGtctc                        <  7:65174/149‑1 (MQ=255)
                                                                                                                            aatcaGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGtct                         <  4:61115/148‑1 (MQ=255)
                                                                                                                                    gcgcGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAg                   <  1:138627/146‑1 (MQ=255)
                                                                                                                                        gAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCgg            >  6:49029/1‑149 (MQ=255)
                                                                                                                                               ccGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCGGCCGTATc     <  7:176641/149‑1 (MQ=255)
                                                                                                                                                cGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGACGACATCGGCCGTATc     >  4:49791/1‑148 (MQ=255)
                                                                                                                                                cGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCGGCCGTATCg    >  2:116589/1‑149 (MQ=255)
                                                                                                                                                  aCCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCGGCCGTATCGaa  >  5:168702/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCATGCGCGAGATGACCTTGAGCCTGCGCGACCTGATCGGCGGCATTGATCAGGGTGTCGGCCAGCTGTCACAGGCCGCCACGCAGTTGGCGGCCAGTAGCGAAGACACCAAGCTGCGCATCAATCAGCAGCGCGAGGAAACCGACCAGGTTGCGACTGCCATGAACCAGATGAGCGCCACCGTGCAGGAAGTGGCACAGAACGCCGAGCAAGCCTCCCTGGCCGCGACCAACGCTGACCAGCAGGCGCAGGTGGGTGACCAGGTAGTCTCCGAAGCCATCGGCCGTATCGAA  >  NC_002947/5720543‑5720837

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