Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,717 Δ1 bp 8.3% coding (703/867 nt) yneJ → LysR family DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,571,7170C. 8.3% 217.4 / 15.0 72coding (703/867 nt)yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base . (3/3);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

CTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571569‑3571866
                                                                                                                                                    |                                                                                                                                                       
cTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGTCATGCTGcc                                                                                                                                                        <  2:654504/150‑1 (MQ=255)
   cGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgc‑g                                                                                                                                                    >  1:1762454/1‑150 (MQ=255)
    gCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCgc‑gc                                                                                                                                                   >  1:898248/1‑150 (MQ=255)
     ccTTCCGCGCCAACTGTTCATACCGCCATCACTTAGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑gca                                                                                                                                                  <  2:1065647/150‑1 (MQ=255)
        tCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑gcagca                                                                                                                                               <  2:360666/150‑1 (MQ=255)
           gcgcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTG‑CGCGGCAGCATgg                                                                                                                                            <  2:224079/150‑2 (MQ=255)
             gcCAACTGTTCATACCGCCACCTCTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTc                                                                                                                                          <  1:574362/150‑1 (MQ=255)
                aaCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGAc                                                                                                                                       >  1:1291263/1‑150 (MQ=255)
                aaCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGAc                                                                                                                                       >  2:1175649/1‑150 (MQ=255)
                 aCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACa                                                                                                                                      >  2:237195/1‑150 (MQ=255)
                   tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCt                                                                                                                                           <  2:1241851/143‑1 (MQ=255)
                   tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCt                                                                                                                                           >  1:1241851/1‑143 (MQ=255)
                         tACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGccc                                                                                                                              >  2:1733588/1‑150 (MQ=255)
                         tACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGccc                                                                                                                              >  2:1899662/1‑150 (MQ=255)
                          aCCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCg                                                                                                                             >  2:358156/1‑150 (MQ=255)
                             gCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGt                                                                                                                          <  1:1175649/150‑1 (MQ=255)
                                aTCTCTTCGAAAACTGGTTCGTCTAGGACCTGGCGGTACCGGGCAAAATTCGCGAAATGGGGTCGTACCACGGCCTGTTGGCCTGCGTCAGTGCGGGTTCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTTGACAATATGCCCGGGTgcc                                                                                                                       <  2:893379/150‑2 (MQ=255)
                                     ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGTCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgcagagat                                                                                                                  >  1:893379/1‑146 (MQ=255)
                                     ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTG‑CGCGGCAGCATg                                                                                                                                             <  2:901524/123‑1 (MQ=255)
                                     ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTG‑CGCGGCAGCATg                                                                                                                                             >  1:901524/1‑123 (MQ=255)
                                        gAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTc                                                                                                               >  2:333366/1‑150 (MQ=255)
                                          aaaCTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAg                                                                                                             <  1:195891/150‑1 (MQ=255)
                                             cTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAgcgc                                                                                                          <  1:358156/150‑1 (MQ=255)
                                             cTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAgcgc                                                                                                          >  2:58230/1‑150 (MQ=255)
                                              tGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGcc                                                                                                         >  2:1633349/1‑150 (MQ=255)
                                                 ttCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTgg                                                                                                      <  1:58230/150‑1 (MQ=255)
                                                     tGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATg                                                                                                                                             <  2:213888/107‑1 (MQ=255)
                                                     tGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATg                                                                                                                                             >  1:213888/1‑107 (MQ=255)
                                                          gACCAGGCGGTACCGGGCAAGATTCACGGAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             <  2:1087402/150‑1 (MQ=255)
                                                            ccAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           <  1:237195/150‑1 (MQ=255)
                                                             cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCggag                                                                                          >  2:1330983/1‑150 (MQ=255)
                                                               ggCGGTACCAGGCAAGATTCACGAAAGGGGGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGCCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTCCAGCACCGTCAGCGCCTGGCCAATGTCggagga                                                                                        <  2:186652/150‑1 (MQ=255)
                                                                gCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGAc                                                                                       >  2:798666/1‑150 (MQ=255)
                                                                gCGGTACCGGGCAAAATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGGCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGAc                                                                                       >  2:71126/1‑150 (MQ=255)
                                                                     aCCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCg                                                                                  >  1:835403/1‑150 (MQ=255)
                                                                      ccGGGCGAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGc                                                                                 >  2:1294131/1‑150 (MQ=255)
                                                                      ccGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgca                                                                                                                       >  1:213270/1‑112 (MQ=255)
                                                                      ccGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgca                                                                                                                       <  2:213270/112‑1 (MQ=255)
                                                                      ccGGGCAAGATTCACGAAATGGAGTCGTACCACAGCGCTTTGTCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGc                                                                                 <  1:98122/150‑1 (MQ=255)
                                                                        gggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGATGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTa                                                                               >  2:838188/1‑150 (MQ=255)
                                                                         ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCTGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             <  1:1067248/135‑1 (MQ=255)
                                                                         ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             >  2:1067248/1‑135 (MQ=255)
                                                                           cAAGATTCCCTAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCt                                                                            <  1:1516807/150‑1 (MQ=255)
                                                                           cAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCt                                                                            <  1:1516766/150‑1 (MQ=255)
                                                                                   aCGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACct                                                                   <  2:1560729/150‑1 (MQ=255)
                                                                                    cGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGTCCATGCTGCCGC‑GCTGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACctg                                                                  >  1:1382114/1‑150 (MQ=255)
                                                                                     gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACctgg                                                                 <  1:1294131/150‑1 (MQ=255)
                                                                                     gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACctgg                                                                 <  1:1633349/150‑1 (MQ=255)
                                                                                       aaTGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCGTGCTCGACAATATGCCCGGGGGCGGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTGGCCTATCTGA‑AAACctggct                                                               >  1:1038741/1‑150 (MQ=255)
                                                                                           gAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTa                                                           <  2:835403/150‑1 (MQ=255)
                                                                                            aGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTAt                                                          <  1:1109158/150‑1 (MQ=255)
                                                                                                gTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGcg                                                      >  1:1840372/1‑150 (MQ=255)
                                                                                                  aCCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGcgac                                                    >  1:1838736/1‑150 (MQ=255)
                                                                                                   ccACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCAGAGGACTTCCGCTATCTGA‑AAACCTGGCTGTTCAGGcgaca                                                   <  2:470941/150‑2 (MQ=255)
                                                                                                       ggCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATTTGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGt                                               >  2:417905/1‑150 (MQ=255)
                                                                                                          aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACg                                            <  1:1733588/150‑1 (MQ=255)
                                                                                                               ggCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGtctc                                       <  1:1330983/150‑1 (MQ=255)
                                                                                                                gCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTcg                                      >  2:75461/1‑150 (MQ=255)
                                                                                                                  cTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGACAAACCT‑GCTTATATGGCGACGAGGTACGGTCTcgcg                                    <  1:1483240/150‑1 (MQ=255)
                                                                                                                          gTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGa                            >  2:125165/1‑150 (MQ=255)
                                                                                                                                    cGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCGAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGt                  >  2:1000376/1‑150 (MQ=255)
                                                                                                                                       ccTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGaa               <  2:1838736/150‑1 (MQ=255)
                                                                                                                                        tcgagcATGCTGCGCG‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTAGTATGGCGACGAGGTACGGTCTAGCGCAGTTTGAGCATGTTCGTGAAg              >  2:806550/6‑150 (MQ=255)
                                                                                                                                        cTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGCCTTCCGCTATCTGA‑AAACCTGGCTGGTAGGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGGGAAg              >  2:245561/1‑150 (MQ=255)
                                                                                                                                         tGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTg                             <  2:409531/134‑1 (MQ=255)
                                                                                                                                         tGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTg                             >  1:409531/1‑134 (MQ=255)
                                                                                                                                          gagcATGCTG‑CGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           <  1:1521920/69‑1 (MQ=37)
                                                                                                                                          gagcATGCTG‑CGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           >  2:1521920/4‑72 (MQ=37)
                                                                                                                                          gagcATGCTG‑CGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGGGGAGTTTAAACAGCTGA‑AAACCTGGCTGGTATGGCGACGAGGAACGGAATCGCGCAGGTTGAGCATGTTCGTGAAGtt            >  1:512339/4‑150 (MQ=255)
                                                                                                                                           gCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTg           >  1:1896882/1‑150 (MQ=255)
                                                                                                                                             cATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt         <  1:71126/150‑1 (MQ=255)
                                                                                                                                             cATGCTGCCGC‑GCAGCATGCTCCACACTATGCCCGGTTGCAGCACCGTCAGCGCCTGGCCATTGTCGAATACCACCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt         <  1:798666/150‑1 (MQ=255)
                                                                                                                                              aTGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCtt        <  1:125165/150‑1 (MQ=255)
                                                                                                                                                  tgcCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACctggc                                                                >  1:149071/1‑90 (MQ=255)
                                                                                                                                                  tgcCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACctggc                                                                <  2:149071/90‑1 (MQ=255)
                                                                                                                                                    ccGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaagaa  <  1:758385/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571569‑3571866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑