Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,717 Δ1 bp 7.4% 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. 7.4% 210.0 / 11.5 68coding (703/867 nt)yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base . (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571568‑3571866
                                                                                                                                                     |                                                                                                                                                      
tCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCAtgctgc                                                                                                                                                        >  2:1820481/1‑150 (MQ=255)
   aCGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgc                                                                                                                                                     >  2:1977786/1‑150 (MQ=255)
     gCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCgc‑gc                                                                                                                                                  >  2:267765/1‑150 (MQ=255)
     gCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCgc‑gc                                                                                                                                                  <  1:267765/150‑1 (MQ=255)
         tCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑gcagcc                                                                                                                                              >  2:992427/1‑149 (MQ=255)
           cgcgCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCG‑GCAGCATg                                                                                                                                            <  1:720099/150‑1 (MQ=255)
            gcgcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGc                                                                                                                                           >  1:263285/1‑150 (MQ=255)
                      ttCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAatat                                                                                                                                 <  2:1419850/150‑1 (MQ=255)
                                cATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTTGGTGCAGTTGCTGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgc                                                                                                                       <  2:809925/150‑1 (MQ=255)
                                cATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgc                                                                                                                       <  2:428351/150‑1 (MQ=255)
                                 aTCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGTTgca                                                                                                                      <  2:221450/150‑1 (MQ=255)
                                  tCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGcc                                                                                                                                                       <  2:2269352/117‑1 (MQ=255)
                                  tCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGcc                                                                                                                                                       >  1:2269352/1‑117 (MQ=255)
                                  tCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTgcag                                                                                                                     >  2:1949248/1‑150 (MQ=255)
                                      ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTG‑CGCGGCAGCATg                                                                                                                                            <  1:1610456/123‑1 (MQ=255)
                                      ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTG‑CGCGGCAGCATg                                                                                                                                            >  2:1610456/1‑123 (MQ=255)
                                          aaaaCTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGa                                                                                                                                       <  2:1514733/124‑1 (MQ=255)
                                          aaaaCTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGa                                                                                                                                       >  1:1514733/1‑124 (MQ=255)
                                                 gTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                      >  1:2024956/1‑150 (MQ=255)
                                                   tCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGc                                                                                                    >  1:1219543/1‑150 (MQ=255)
                                                      tGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGTGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCaa                                                                                                 <  1:865318/150‑1 (MQ=255)
                                                       gCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAAt                                                                                                <  1:1442567/150‑1 (MQ=255)
                                                           gACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                            >  1:2386172/1‑150 (MQ=255)
                                                             ccAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                          >  1:340771/1‑150 (MQ=255)
                                                             ccAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                          <  2:340771/150‑1 (MQ=255)
                                                               aGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                          >  2:1711474/1‑148 (MQ=255)
                                                               aGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                          <  1:1711474/148‑1 (MQ=255)
                                                                ggCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGa                                                                                                                                       <  1:53605/102‑1 (MQ=255)
                                                                ggCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGa                                                                                                                                       >  2:53605/1‑102 (MQ=255)
                                                                     tACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTcc                                                                                  <  1:313120/150‑1 (MQ=255)
                                                                         gggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTa                                                                              >  1:2237147/1‑150 (MQ=255)
                                                                          ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCATCCGTAACTAGAGGCGCCATCCTGGCCATGCGGCCGC‑GCAGCATGCTCGACAAGATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                             >  1:648073/1‑150 (MQ=255)
                                                                           gCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATc                                                                            <  2:1459625/150‑1 (MQ=255)
                                                                               gATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGaa                                                                        >  1:1322539/1‑150 (MQ=255)
                                                                                  tCACGAAATGGAGTCGTACCACGGGATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCATCATGCTCGACAATATGCCCGGGt                                                                                                                         <  2:1185941/98‑1 (MQ=255)
                                                                                  tCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGt                                                                                                                         >  1:1185941/1‑98 (MQ=255)
                                                                                    aCGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGGCATGCTGCGCG‑GCAGCATggcca                                                                                                                                        <  1:573342/81‑5 (MQ=37)
                                                                                    aCGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGGCATGCTG‑CGCGGCAGCATggcca                                                                                                                                        >  2:573342/1‑77 (MQ=37)
                                                                                      gAAATGGAGTCGTACCACGGCCTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctgg                                                                 <  1:1820481/150‑1 (MQ=255)
                                                                                      gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctgg                                                                 <  2:495549/150‑1 (MQ=255)
                                                                                          tGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggctgg                                                             <  1:1538995/150‑1 (MQ=255)
                                                                                               tCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATgg                                                        >  1:152590/1‑150 (MQ=255)
                                                                                                cGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGc                                                       <  2:2237147/150‑1 (MQ=255)
                                                                                                 gTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCATGGCCAATGTCGGAGGACGTCCGCTATCTGAAAACCTGGCTGGTATGGcg                                                      >  2:1249288/1‑150 (MQ=255)
                                                                                                  tACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                      <  1:774776/101‑1 (MQ=255)
                                                                                                  tACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                      >  2:774776/1‑101 (MQ=255)
                                                                                                     cACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcgacga                                                  >  1:965789/1‑150 (MQ=255)
                                                                                                     cACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcgacga                                                  >  1:1626954/1‑150 (MQ=255)
                                                                                                      aCGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAg                                                 <  1:37677/150‑1 (MQ=255)
                                                                                                        ggCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGt                                               <  2:2024956/150‑1 (MQ=255)
                                                                                                          cATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTAc                                             <  2:965789/150‑1 (MQ=255)
                                                                                                           aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GACGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACg                                            <  1:1968304/150‑1 (MQ=255)
                                                                                                            tGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACgg                                           <  2:1626954/150‑1 (MQ=255)
                                                                                                                 gCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAatat                                                                                                                                 >  2:1358764/1‑59 (MQ=255)
                                                                                                                 gCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAatat                                                                                                                                 <  1:1358764/59‑1 (MQ=255)
                                                                                                                   cTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcgcg                                    >  1:1518229/1‑150 (MQ=255)
                                                                                                                         cAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTAc                                             >  2:1289788/1‑135 (MQ=255)
                                                                                                                         cAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTAc                                             <  1:1289788/135‑1 (MQ=255)
                                                                                                                           gTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGaa                                                                        >  2:869394/1‑106 (MQ=255)
                                                                                                                           gTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGaa                                                                        <  1:869394/106‑1 (MQ=255)
                                                                                                                               gggTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGANCAt                        <  2:1518229/150‑1 (MQ=255)
                                                                                                                                ggTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATgg                                                        >  1:1067712/1‑117 (MQ=255)
                                                                                                                                ggTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATgg                                                        <  2:1067712/117‑1 (MQ=255)
                                                                                                                                   gCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTc                    <  1:1977786/150‑1 (MQ=255)
                                                                                                                                       gCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGa                <  1:992427/150‑1 (MQ=255)
                                                                                                                                        ccTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGaa               >  1:822093/1‑150 (MQ=255)
                                                                                                                                         cTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAg              >  2:119654/1‑150 (MQ=255)
                                                                                                                                         cTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTGTGAGCATGTTCGTGAAg              >  1:213173/1‑150 (MQ=255)
                                                                                                                                              cATGCTG‑CGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt         <  1:1366338/150‑1 (MQ=255)
                                                                                                                                              cATGCTG‑CGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt         >  2:1366338/1‑150 (MQ=255)
                                                                                                                                               aTGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCtt        >  1:1201583/1‑150 (MQ=255)
                                                                                                                                                tgctgcCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTg       >  2:101957/1‑150 (MQ=255)
                                                                                                                                                     ccGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaagaa  <  2:1201583/150‑1 (MQ=255)
                                                                                                                                                     ccGC‑CCAGCATGATCTACAATAGGCTCGGGTGCAGCACCTTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGTCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaagaa  <  2:1322539/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGC‑GCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAA  >  NC_002947/3571568‑3571866

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