Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,286,687 G→C 9.1% A350G (GCC→GGC)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,286,6870GC 9.1% 197.4 / 8.5 66A350G (GCC→GGC) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (2/4);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

TCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATT  >  NC_002947/4286543‑4286830
                                                                                                                                                |                                                                                                                                               
tCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGc                                                                                                                                             >  2:70721/1‑149 (MQ=255)
 cTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGc                                                                                                                                             <  7:101058/148‑1 (MQ=255)
  tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGCCCTCCGCTGGCCTGCTGGGCAGCGa                                                                                                                                           >  8:59360/1‑149 (MQ=255)
  tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGa                                                                                                                                           >  3:313179/1‑149 (MQ=255)
    tGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                         >  2:226411/1‑146 (MQ=255)
    tGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgg                                                                                                                                           <  1:70721/147‑2 (MQ=255)
     gAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                         <  1:226411/148‑4 (MQ=255)
     gAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                         <  5:139839/148‑4 (MQ=255)
         gTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                                    >  3:257629/1‑149 (MQ=255)
             ccAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                          >  2:319608/1‑139 (MQ=255)
             ccAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                          <  1:319608/139‑1 (MQ=255)
              cAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTggg                                                                                                                                                 <  8:71152/131‑1 (MQ=255)
               aGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCAACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActg                                                                                                                              >  4:374093/1‑149 (MQ=255)
                gCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTCCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgc                                                                                                                             >  6:47228/1‑149 (MQ=255)
                 cgtattttGCTGGGTTATTGTCCCACTCTTGCCCTCCTCTCTATCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgct                                                                                                                            <  1:136502/143‑1 (MQ=255)
                  ggCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgctg                                                                                                                           <  4:81970/149‑1 (MQ=255)
                          cAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                                   >  6:363389/1‑149 (MQ=255)
                          cAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                                   >  5:108114/1‑149 (MQ=255)
                                        cACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTggcg                                                                                                     <  3:374093/149‑1 (MQ=255)
                                            cTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGgcgc                                                                                                 <  2:282314/149‑1 (MQ=255)
                                               gCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCt                                                                                                <  3:154401/147‑1 (MQ=255)
                                                cccTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgat                                                                                             <  7:119838/149‑1 (MQ=255)
                                                 cctcctCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgata                                                                                            >  7:320342/1‑149 (MQ=255)
                                                  ctcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgatac                                                                                           <  4:257629/149‑1 (MQ=255)
                                                   tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGAAAGATCCGCTGCTGGTTGGCGGCGCTgatag                                                                                           >  5:267058/1‑147 (MQ=255)
                                                   tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGAAAGATCCGCTGCTGGTTGGCGGCGCTga                                                                                              <  6:267058/145‑1 (MQ=255)
                                                   tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCAAAGCCCAGACCTCCGCTGGCCTCCTGGGCACCGACCCGCACGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCCATgatacc                                                                                          >  6:35229/1‑149 (MQ=255)
                                                     ctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccga                                                                                        <  6:108114/149‑1 (MQ=255)
                                                      tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGa                                                                                                                         >  5:289874/1‑115 (MQ=255)
                                                      tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGa                                                                                                                         <  6:289874/115‑1 (MQ=255)
                                                        gCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccgatac                                                                                     >  7:128945/1‑149 (MQ=255)
                                                          tttCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                                    <  6:210508/100‑1 (MQ=255)
                                                          tttCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                                    >  5:210508/1‑100 (MQ=255)
                                                          tttCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGataccgataccg                                                                                   >  3:123921/1‑149 (MQ=255)
                                                            tCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGataccgataccg                                                                                   <  5:363389/147‑1 (MQ=255)
                                                             cAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTg                                                                                 <  7:214021/148‑1 (MQ=255)
                                                              aaCCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGAt                                                                               >  3:357906/1‑149 (MQ=255)
                                                                 cTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcct                                                                            >  2:370921/1‑149 (MQ=255)
                                                                    tATCCGGCCACCGTGCGCCTGGACGATGTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgca                                                                         <  6:40439/149‑1 (MQ=255)
                                                                               cGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                              <  1:370921/149‑1 (MQ=255)
                                                                               cGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                              <  7:59360/149‑1 (MQ=255)
                                                                               cGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                              <  4:123921/149‑1 (MQ=255)
                                                                                      cTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgc                                                       <  4:313179/149‑1 (MQ=255)
                                                                                          actaTCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                                   <  4:34917/82‑1 (MQ=255)
                                                                                          aCGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                                   >  3:34917/1‑85 (MQ=255)
                                                                                            gATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCCTGCCGCCCt                                                 >  8:218831/1‑149 (MQ=255)
                                                                                            gATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCt                                                 >  7:16126/1‑149 (MQ=255)
                                                                                                  cTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCAt                                                         <  6:155535/135‑1 (MQ=255)
                                                                                                  cTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCAt                                                         >  5:155535/1‑135 (MQ=255)
                                                                                                  cTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTg                                                <  8:128945/144‑1 (MQ=255)
                                                                                                  cTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCt                                              <  2:342128/146‑1 (MQ=255)
                                                                                                   tCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGaa                                          >  5:157038/1‑149 (MQ=255)
                                                                                                         aTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                                   <  1:336700/134‑1 (MQ=255)
                                                                                                         aTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                                   >  2:336700/1‑134 (MQ=255)
                                                                                                             aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                              >  1:180139/1‑103 (MQ=255)
                                                                                                             aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                              <  2:180139/103‑1 (MQ=255)
                                                                                                             aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                             <  4:241093/136‑1 (MQ=255)
                                                                                                             aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                             >  3:241093/1‑136 (MQ=255)
                                                                                                              cAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCt                               >  4:359228/1‑149 (MQ=255)
                                                                                                                 ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                           <  1:139829/102‑1 (MQ=255)
                                                                                                                 ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                           >  2:139829/1‑102 (MQ=255)
                                                                                                                        cccAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGAtt                     <  5:47228/149‑1 (MQ=255)
                                                                                                                             aCCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCt                >  5:358401/1‑149 (MQ=255)
                                                                                                                              ccTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAATATGCTGTGATGTTTGGCGGCGCTGATACAGATACAGTGATGCTGGTACTGCTACCAGGCCTTGCAGCCCTTCTCGAACTCCTGCACGTGCACCGGATTGctctt               <  8:320342/149‑2 (MQ=255)
                                                                                                                                       gCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAAc       <  8:16126/148‑1 (MQ=255)
                                                                                                                                         ctgctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTg                              >  6:375868/1‑123 (MQ=255)
                                                                                                                                            ctgGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACgcatt  <  6:157038/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATT  >  NC_002947/4286543‑4286830

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