Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,720:1 +G 8.5% coding (706/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,7201.G 8.5% 269.9 / 19.0 71coding (706/867 nt)yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (33/32);  new base G (3/3);  total (36/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 = 9.33e-01

CGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571572‑3571866
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                       >  2:1175649/1‑150 (MQ=255)
             aaCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGAc                                                                                                                                       >  1:1291263/1‑150 (MQ=255)
              aCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACa                                                                                                                                      >  2:237195/1‑150 (MQ=255)
                tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCt                                                                                                                                           >  1:1241851/1‑143 (MQ=255)
                tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCt                                                                                                                                           <  2:1241851/143‑1 (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                                                                                                                                             >  1:213888/1‑107 (MQ=255)
                                                  tGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATg                                                                                                                                             <  2:213888/107‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
CGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGA‑AAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571572‑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: ‑