Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,488,073 T→A 7.4% Y86F (TAC→TTC)  degS ← serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,488,0730TA 7.4% 265.6 / 9.5 92Y86F (TAC→TTC) degSserine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base A (4/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

TCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATC  >  NC_002947/1487924‑1488219
                                                                                                                                                      |                                                                                                                                                  
taagcagGTAGCCCTCCAGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgt‑gt                                                                                                                                                    >  2:1166225/3‑150 (MQ=255)
 cagtagGTAGCTCTCTGGGCTCATGATTACCGCCGAACTCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTa                                                                                                                                                   <  2:495683/150‑1 (MQ=255)
   gcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACa                                                                                                                                                 >  2:1022160/1‑150 (MQ=255)
    cagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAg                                                                                                                                                <  2:484999/150‑1 (MQ=255)
    cagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAg                                                                                                                                                <  1:948342/150‑1 (MQ=255)
    cagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAg                                                                                                                                                >  1:1869098/1‑150 (MQ=255)
      ggTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCTCTGTTTGGGCAGGTTATCGCCTAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCAAGGTGGT‑GAACaaaa                                                                                                                                              <  2:1173646/150‑4 (MQ=255)
       gTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGtt                                                                                                                                             >  1:1081805/1‑150 (MQ=255)
       gTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTGATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGtt                                                                                                                                             >  2:1537346/1‑150 (MQ=255)
             ctctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCAc                                                                                                                                       >  1:2347453/1‑150 (MQ=255)
                   gCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCgg                                                                                                                                 <  1:2087770/150‑1 (MQ=255)
                   gCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCgg                                                                                                                                 <  2:150238/150‑1 (MQ=255)
                     tcgtGATTACAGCCGAACCCAGGCCTGATTCCCAACGGCGCTGCTTGGGCAGGTTATCGCCGAGGAACCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTAACCACCTTGGTGGT‑GTACAGGTTGGCCACCGGCgggg                                                                                                                               >  2:470700/4‑150 (MQ=255)
                      cATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGggc                                                                                                                              <  1:568377/150‑1 (MQ=255)
                        tGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGAAGGTTATCGCCGAAGACGAGGCGGAAATGAGGGTCTTCGAACAGCGGGTGGTCGCTGTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGggcgg                                                                                                                            <  2:1464282/150‑1 (MQ=255)
                         gATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCAccgccg                                                                                                                                  >  2:7855/1‑143 (MQ=255)
                         gATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCAccgccg                                                                                                                                  <  1:7855/143‑1 (MQ=255)
                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGAggcgacnatc                                                                                                                       <  2:1518722/150‑7 (MQ=255)
                                 ccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGGAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTcg                                                                                                                   >  1:1196461/1‑149 (MQ=255)
                                  cgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAc                                                                                                                  <  1:370206/150‑1 (MQ=255)
                                       ccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGAGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgt                                                                                                             >  2:1035641/1‑150 (MQ=255)
                                         aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                           >  2:215103/1‑150 (MQ=255)
                                         aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTT‑GCGCCCCTGGCGCTCGAACGGGTGGTCCACCCCCTGGCCGGCCAGGGCCAcggcgccg                                                                                                           >  2:1728996/1‑150 (MQ=255)
                                          ggCTTGATTCCCAGCGGCGCTGATTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                                          <  2:1201839/150‑1 (MQ=255)
                                             ttGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcggcgt                                                                                                       >  1:1591162/1‑150 (MQ=255)
                                                 ttCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGCACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGa                                                                                                   <  1:1728996/150‑1 (MQ=255)
                                                 ttCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGa                                                                                                   <  2:2044606/150‑1 (MQ=255)
                                                  tCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGac                                                                                                  >  2:191996/1‑150 (MQ=255)
                                                    ccAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                           >  2:2278534/1‑139 (MQ=255)
                                                    ccAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                           <  1:2278534/139‑1 (MQ=255)
                                                     cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGGCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACAcc                                                                                               >  1:2195235/1‑150 (MQ=255)
                                                      aGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCg                                                                                              <  1:321435/150‑1 (MQ=255)
                                                       gcggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCgg                                                                                             <  2:1668640/150‑1 (MQ=255)
                                                       gcggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCgg                                                                                             <  1:2271789/150‑1 (MQ=255)
                                                             cTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTg                                                                                       >  1:70857/1‑150 (MQ=255)
                                                             cTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTg                                                                                       >  2:987778/1‑150 (MQ=255)
                                                              tGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGc                                                                                      <  1:1791286/150‑1 (MQ=255)
                                                                tttGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCAt                                                                                    <  2:1081805/150‑1 (MQ=255)
                                                                tttGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGACAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCAt                                                                                    >  1:1753125/1‑150 (MQ=255)
                                                                    ggCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATc                                                                                >  2:1547269/1‑150 (MQ=255)
                                                                       aGGTTATCGCCGAAGAAGCGGCGGTACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                             <  1:1524113/150‑1 (MQ=255)
                                                                       aGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGACACGGGGACCTGTTAGACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATGCgg                                                                             >  2:1401876/1‑150 (MQ=255)
                                                                       aGGTTACCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                             <  1:1826005/150‑1 (MQ=255)
                                                                        ggTTATCTCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGACAGGGTCACGGCGTCGGCGTACGACACCGGACCTTGCATGATCCggg                                                                            <  2:1225149/150‑1 (MQ=255)
                                                                             tCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTc                                                                       >  1:377528/1‑150 (MQ=255)
                                                                                ccGAATAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCtgt                                                                    >  2:822278/1‑150 (MQ=255)
                                                                                      aagCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCt                                                              <  2:1021452/150‑1 (MQ=255)
                                                                                        gcggcgGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtg                                                            <  1:114337/150‑1 (MQ=255)
                                                                                        gcggcgGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtg                                                            >  1:2098172/1‑150 (MQ=255)
                                                                                           gcggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcg                                                                                                                 <  1:1868727/94‑1 (MQ=255)
                                                                                           gcggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcg                                                                                                                 >  2:1868727/1‑94 (MQ=255)
                                                                                             ggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCa                                                       >  1:2040274/1‑150 (MQ=255)
                                                                                             ggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCa                                                       >  1:1433847/1‑150 (MQ=255)
                                                                                              gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                      <  1:191996/150‑1 (MQ=255)
                                                                                              gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGACAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                      <  2:1753125/150‑1 (MQ=255)
                                                                                               aaCTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAgg                                                     <  1:987778/150‑1 (MQ=255)
                                                                                                   gCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGa                                                 >  2:1102122/1‑150 (MQ=255)
                                                                                                      ggTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGt                                              >  1:1252231/1‑150 (MQ=255)
                                                                                                       gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCggg                                                                            <  2:2043098/119‑1 (MQ=255)
                                                                                                       gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCggg                                                                            >  1:2043098/1‑119 (MQ=255)
                                                                                                       gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtg                                                            >  1:1746652/1‑135 (MQ=255)
                                                                                                       gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtg                                                            <  2:1746652/135‑1 (MQ=255)
                                                                                                       gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTc                                             >  1:2407833/1‑150 (MQ=255)
                                                                                                        tCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTcc                                            <  2:1252231/150‑1 (MQ=255)
                                                                                                          ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctg                                          >  2:2125953/1‑150 (MQ=255)
                                                                                                          ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctg                                          >  1:1220919/1‑150 (MQ=255)
                                                                                                           tCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTACACTGTGGCACATCCGGCAGGTTGACGTCctgg                                         >  1:1945640/1‑150 (MQ=255)
                                                                                                              aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctg                                      <  1:1537346/150‑1 (MQ=255)
                                                                                                              aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctg                                      <  1:1547269/150‑1 (MQ=255)
                                                                                                                 aGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGc                                   >  2:338910/1‑150 (MQ=255)
                                                                                                                  gCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTCCTGCAGGTTGTCGTCCTCGCTGGGCa                                  >  1:1000571/1‑150 (MQ=255)
                                                                                                                     ggTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGcc                               <  2:2347453/150‑1 (MQ=255)
                                                                                                                         ggCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGNGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGAGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAc                           >  1:599198/1‑150 (MQ=255)
                                                                                                                          gcgcTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAcc                          >  1:697347/1‑150 (MQ=255)
                                                                                                                              ttttGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                      >  1:2288800/1‑150 (MQ=255)
                                                                                                                              ttttGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                      >  2:136383/1‑150 (MQ=255)
                                                                                                                               cctgtacaccacCAAGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGCGGACAGAGGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTg                     >  1:1971581/7‑150 (MQ=255)
                                                                                                                                ctgtacaccacCAAGGTGGT‑GAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                                                >  1:1134441/6‑122 (MQ=255)
                                                                                                                                ctgtacaccacCAAGGTGGT‑GAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGc                    >  1:2185423/6‑150 (MQ=255)
                                                                                                                                ctGTACACCACCAAGGTGGT‑GAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                                                <  2:1134441/121‑1 (MQ=255)
                                                                                                                                 tgtacaccacCAAGGTGGT‑GAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCg                             <  2:418495/136‑1 (MQ=255)
                                                                                                                                 tGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCg                   <  2:2195235/150‑1 (MQ=255)
                                                                                                                                 tGTACACCACCAAGGTGGT‑GAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCg                             >  1:418495/1‑140 (MQ=255)
                                                                                                                                   ttcgacaccTTGGTGGT‑GTACAGGTTAGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGa                 >  2:691265/6‑150 (MQ=255)
                                                                                                                                    tCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaa                <  1:230210/150‑1 (MQ=255)
                                                                                                                                      accaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTTTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa              >  1:983179/1‑150 (MQ=255)
                                                                                                                                      accaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa              <  1:2125953/150‑1 (MQ=255)
                                                                                                                                      accaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa              >  2:891161/1‑150 (MQ=255)
                                                                                                                                       ccaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGTCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAccc                         >  2:1401495/1‑138 (MQ=255)
                                                                                                                                       ccaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                                    <  1:1983499/127‑1 (MQ=255)
                                                                                                                                       ccaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                                    >  2:1983499/1‑127 (MQ=255)
                                                                                                                                       ccaccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAccc                         <  1:1401495/138‑1 (MQ=255)
                                                                                                                                          ccTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt          >  1:1754227/1‑150 (MQ=255)
                                                                                                                                                  gt‑gtACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTAGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATc  <  1:1022160/150‑1 (MQ=255)
                                                                                                                                                  gt‑gtACAGGTTAGCAACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATc  >  2:700133/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGT‑GTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATC  >  NC_002947/1487924‑1488219

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