Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,286,684 T→A 9.2% Q351L (CAG→CTG)  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,6840TA 9.2% 185.8 / 7.4 65Q351L (CAG→CTG) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base A (2/4);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

GTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATT  >  NC_002947/4286536‑4286830
                                                                                                                                                    |                                                                                                                                                  
gTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCctgct                                                                                                                                                    <  3:292637/149‑1 (MQ=255)
 ttttGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCctgctg                                                                                                                                                   <  7:131503/149‑1 (MQ=255)
       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                                                                                                                                    >  5:210508/1‑100 (MQ=255)
                                                                 tttCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                                    <  6:210508/100‑1 (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                                                   >  2:336700/1‑134 (MQ=255)
                                                                                                                aTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                                   <  1:336700/134‑1 (MQ=255)
                                                                                                                    aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                              >  1:180139/1‑103 (MQ=255)
                                                                                                                    aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                              <  2:180139/103‑1 (MQ=255)
                                                                                                                    aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                             >  3:241093/1‑136 (MQ=255)
                                                                                                                    aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                             <  4:241093/136‑1 (MQ=255)
                                                                                                                     cAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCt                               >  4:359228/1‑149 (MQ=255)
                                                                                                                        ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                           >  2:139829/1‑102 (MQ=255)
                                                                                                                        ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                           <  1:139829/102‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
GTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGATCTGCTCGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAACGCATT  >  NC_002947/4286536‑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: ‑