Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,947,722 T→C 31.1% K251R (AAG→AGG)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,947,7220TC31.1% 73.4 / 45.0 61K251R (AAG→AGG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base C (10/9);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

AGCGCTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGA  >  NC_002947/5947576‑5947857
                                                                                                                                                  |                                                                                                                                       
agCGCTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTg                                                                                                                                       <  1:227341/149‑1 (MQ=255)
   gcTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGcc                                                                                                                                    <  5:334111/149‑1 (MQ=255)
   gcTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGcc                                                                                                                                    <  5:339458/149‑1 (MQ=255)
    cTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCt                                                                                                                                   >  3:218410/1‑149 (MQ=255)
       gggCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTtct                                                                                                                                >  7:132280/1‑149 (MQ=255)
       gggCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTtc                                                                                                                                 <  4:218410/148‑1 (MQ=255)
       gggCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGA‑CAGGCGCTTGGCCTtc                                                                                                                                 <  1:254468/147‑1 (MQ=255)
         gCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTtctc                                                                                                                               >  5:365415/1‑148 (MQ=255)
             ctcGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATc                                                                                                                           <  4:44475/148‑1 (MQ=255)
                gAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCAc                                                                                                                         <  5:11854/147‑1 (MQ=255)
                  gggcgtgggcgtCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCAtc                                                                                                                     >  6:45922/1‑149 (MQ=255)
                    gcgtgggcgtCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCAtctc                                                                                                                   <  8:132280/149‑1 (MQ=255)
                          gcgtCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAAc                                                                                                             >  5:129307/1‑149 (MQ=255)
                          gcgtCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAAc                                                                                                             >  5:13999/1‑149 (MQ=255)
                                acacCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCAc                                                                                                       <  4:218415/149‑1 (MQ=255)
                                acacCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCAc                                                                                                       >  3:84930/1‑149 (MQ=255)
                                       cGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCgg                                                                                                <  6:13999/149‑1 (MQ=255)
                                         gTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGt                                                                                              >  8:88387/1‑149 (MQ=255)
                                           cAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTg                                                                                            >  7:264179/1‑149 (MQ=255)
                                             ggatcTTTCCGGAGCTCGGGATCACCGTGTTGTAGGCACGCTCCAGTCGAGTTATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCTACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCTGGGTggc                                                                                          <  4:28208/145‑1 (MQ=255)
                                                 cTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCt                                                                                      >  8:198314/1‑149 (MQ=255)
                                                 cTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCt                                                                                      <  2:151377/149‑1 (MQ=255)
                                                     ccGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCtcgtc                                                                                  >  4:196273/1‑149 (MQ=255)
                                                        gAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGaa                                                                               <  2:176062/149‑1 (MQ=255)
                                                        gAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGaa                                                                               <  4:60219/149‑1 (MQ=255)
                                                             cGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTg                                                                           >  6:260519/1‑148 (MQ=255)
                                                              gggATAACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTgg                                                                          >  4:142098/1‑148 (MQ=255)
                                                               ggATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCATCAGGATAACCACATCTTTCTTGTGCTCGCCCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGaa                                                                        <  3:211287/149‑1 (MQ=255)
                                                                gATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGCCCAGGCGCCTGGCCTTCTCGAACACCATCTCGGCACCCTGCACGTGGCGGGGTGGCGGCTCGTCGAAGGTGGAAg                                                                       >  1:294966/1‑149 (MQ=255)
                                                                 aTCACGGGGGTGTAGGCACGAGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGc                                                                      <  8:146342/149‑1 (MQ=255)
                                                                    aCGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAAc                                                                   <  3:101850/149‑1 (MQ=255)
                                                                    aCGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAAc                                                                   <  7:88387/149‑1 (MQ=255)
                                                                    aCGGTGTGGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAAc                                                                   <  2:289198/149‑1 (MQ=255)
                                                                            gtAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCg                                                            <  4:84930/148‑1 (MQ=255)
                                                                                                   tCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCg                                    >  6:328659/1‑149 (MQ=255)
                                                                                                   tCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTAGGCAACCTGCACGTGGCGGGTAGGCGGCTCGTCGAAGGTGGACGCACCCACTTCGCCCCGCACGGTGCGCTGCATTTCg                                    >  8:331294/1‑149 (MQ=255)
                                                                                                           agcagGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCActtc                            >  5:201719/1‑149 (MQ=255)
                                                                                                                 aTAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCgg                      <  3:196273/149‑1 (MQ=255)
                                                                                                                     ccACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGgcgc                  <  5:45922/149‑1 (MQ=255)
                                                                                                                     ccACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGgcgc                  >  4:297490/1‑149 (MQ=255)
                                                                                                                     ccACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGgcgc                  >  5:145832/1‑149 (MQ=255)
                                                                                                                      cacaTCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCt                 >  6:133576/1‑149 (MQ=255)
                                                                                                                      cacaTCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCt                 >  7:324019/1‑149 (MQ=255)
                                                                                                                       acaTCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtc                <  6:145832/149‑1 (MQ=255)
                                                                                                                       acaTCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCgg                      <  2:87351/143‑1 (MQ=255)
                                                                                                                       acaTCTTTCTTGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCgg                      >  1:87351/1‑143 (MQ=255)
                                                                                                                        caTCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcg               <  5:260519/149‑1 (MQ=255)
                                                                                                                          tctttcttGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcgtc             >  2:339783/1‑149 (MQ=255)
                                                                                                                            tttcttGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGa           >  6:50013/1‑149 (MQ=255)
                                                                                                                            tttcttGTGCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTAGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGa           >  6:307959/1‑149 (MQ=255)
                                                                                                                              tcttGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcgtcg            <  8:324019/146‑1 (MQ=255)
                                                                                                                              tcttGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcgtc             <  6:129307/145‑1 (MQ=255)
                                                                                                                              tcttGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcgtc             <  6:365415/145‑1 (MQ=255)
                                                                                                                                  gtgCTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATcagc      <  1:339783/148‑1 (MQ=255)
                                                                                                                                    gCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATcagcag    >  1:265870/1‑148 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGa  <  5:133576/149‑1 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGa  <  2:265870/149‑1 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGa  >  2:213734/1‑149 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGa  <  3:297490/149‑1 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGt                                  <  5:186810/117‑1 (MQ=255)
                                                                                                                                     cTCGACCAGGCGCCTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGt                                  >  6:186810/1‑117 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
AGCGCTTGGGCTTCTCGAGGGCGTGGGCGTCGACACCACCGGTCAGCACCTTGCCGGAGCTCGGGATCACGGTGTTGTAGGCACGCGCCAGTCGAGTGATCGAGTCCAGCAGGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGA  >  NC_002947/5947576‑5947857

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