Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,053,264 G→A 9.4% intergenic (+22/+24) PP_2664 → / ← agmR two‑component system sensor histidine kinase/response regulator/glycerol metabolism activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,053,2640GA 9.4% 235.5 / 12.2 85intergenic (+22/+24)PP_2664/agmRtwo‑component system sensor histidine kinase/response regulator/glycerol metabolism activator
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (4/4);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

GGTCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTGGGCCTTGACGG  >  NC_002947/3053116‑3053397
                                                                                                                                                    |                                                                                                                                     
ggTCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGa                                                                                                                                      >  2:279397/1‑150 (MQ=255)
  tCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATc                                                                                                                                    <  2:179974/150‑1 (MQ=255)
    aaCTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACAAGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCcg                                                                                                                                  <  2:1054242/150‑1 (MQ=255)
      cTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCcgcg                                                                                                                                >  1:731519/1‑150 (MQ=255)
        aCAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGaa                                                                                                                              <  1:721135/150‑1 (MQ=255)
         cAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAg                                                                                                                             <  1:88590/150‑1 (MQ=255)
          aGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGATCAACCCGCGAAGc                                                                                                                            <  2:110830/150‑1 (MQ=255)
           gCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCg                                                                                                                           >  1:1538941/1‑150 (MQ=255)
            cAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCgg                                                                                                                          >  2:145912/1‑150 (MQ=255)
            cAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCgg                                                                                                                          >  1:1590215/1‑150 (MQ=255)
             aaCGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGaa                                                                                                                              >  1:1334342/1‑145 (MQ=255)
             aaCGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGaa                                                                                                                              <  2:1334342/145‑1 (MQ=255)
              aCGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGcc                                                                                                                        >  2:892524/1‑150 (MQ=255)
                   ctcAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAAc                                                                                                                   <  2:160188/150‑1 (MQ=255)
                    tcAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACa                                                                                                                  <  1:1520559/150‑1 (MQ=255)
                      aaGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGc                                                                                                                >  1:828616/1‑150 (MQ=255)
                            cagATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACg                                                                                                          <  2:1646873/150‑1 (MQ=255)
                             agATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACgg                                                                                                         <  2:1097622/150‑1 (MQ=255)
                              gATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGgc                                                                                                        <  2:731519/150‑1 (MQ=255)
                                   gcgAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGa                                                                                                   >  1:71744/1‑150 (MQ=255)
                                   gcgAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGa                                                                                                   <  1:1066865/150‑1 (MQ=255)
                                     gAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATa                                                                                                 <  2:1590215/150‑1 (MQ=255)
                                         tGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCa                                                                                             <  2:71744/150‑1 (MQ=255)
                                             tCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGgcgc                                                                                                      >  1:1570230/1‑137 (MQ=255)
                                             tCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGgcgc                                                                                                      <  2:1570230/137‑1 (MQ=255)
                                               acacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTTGCTGGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCaaaa                                                                                       <  1:1005674/150‑1 (MQ=255)
                                               acacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCaaaa                                                                                       >  1:54487/1‑150 (MQ=255)
                                               acacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCACTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCaaaa                                                                                       <  1:1005673/150‑1 (MQ=255)
                                               acacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGACCCATCTGCTGGCCCGTCCCGGGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCaaaa                                                                                       <  1:1005474/150‑1 (MQ=255)
                                                cacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAAt                                                                                      >  1:1050514/1‑150 (MQ=255)
                                                cacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAAt                                                                                      <  2:54487/150‑1 (MQ=255)
                                                 acGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAg                                                                                            >  2:1544315/1‑143 (MQ=255)
                                                 acGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAg                                                                                            <  1:1544315/143‑1 (MQ=255)
                                                 acGCTGATCCACAAGCCCGTAAGACCGATGAAACTCAAAACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATc                                                                                     <  1:892524/150‑1 (MQ=255)
                                                           acaAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCa                                                                           >  2:1064712/1‑150 (MQ=255)
                                                             aaGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCa                                                                                                     >  2:1880844/1‑122 (MQ=255)
                                                             aaGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCa                                                                                                     <  1:1880844/122‑1 (MQ=255)
                                                                   gTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGGTTTTGGCAGCACTGAg                                                                   <  2:1050514/150‑1 (MQ=255)
                                                                     aaGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATc                                                                                     >  2:1618848/1‑130 (MQ=255)
                                                                     aaGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATc                                                                                     <  1:1618848/130‑1 (MQ=255)
                                                                       gACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGc                                                                                           <  1:1136532/122‑1 (MQ=255)
                                                                       gACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGc                                                                                           >  2:1136532/1‑122 (MQ=255)
                                                                        aCCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATgg                                                              <  2:1538941/150‑1 (MQ=255)
                                                                          cGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCt                                                            <  1:279397/150‑1 (MQ=255)
                                                                              gAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGa                                                                                   >  2:693059/1‑123 (MQ=255)
                                                                              gAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGa                                                                                   <  1:693059/123‑1 (MQ=255)
                                                                                aaCTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTAcc                                                      >  1:1132514/1‑150 (MQ=255)
                                                                                aaCTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGATTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTAcc                                                      >  1:1570690/1‑150 (MQ=255)
                                                                                   tCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGAATCGAGGTTGGTTGTTTCGCGGGTTGATCAGCTAAGCGGCCGGAACAGCTAAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCgg                                                   >  1:1492755/1‑150 (MQ=255)
                                                                                     aaGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGtt                                                 <  2:1647187/150‑1 (MQ=255)
                                                                                     aaGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGtt                                                 <  2:1492755/150‑1 (MQ=255)
                                                                                      aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCa                                                                           >  2:1239257/1‑123 (MQ=255)
                                                                                      aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCa                                                                           <  1:1239257/123‑1 (MQ=255)
                                                                                      aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGAAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtg                                                >  2:1357507/1‑150 (MQ=255)
                                                                                       gACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgt                                               <  1:1064712/150‑1 (MQ=255)
                                                                                        aCGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGAAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgtg                                              <  1:1357507/150‑1 (MQ=255)
                                                                                         cGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGATCAACCCGCGAAGCGGccaac                                                                                                                     >  2:460120/1‑75 (MQ=21)
                                                                                         cGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGATCAACCCGCGAAGCGGccaac                                                                                                                     <  1:460120/78‑4 (MQ=21)
                                                                                             gATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACAt                                         >  2:1318830/1‑150 (MQ=255)
                                                                                             gATGAGCCATCTGCTGGCCAGTGGCTAGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGa                                                                                                   <  1:787434/92‑1 (MQ=255)
                                                                                             gATGAGCCATCTGCTGGCCAGTGGCTAGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGa                                                                                                   >  2:787434/1‑92 (MQ=255)
                                                                                                gAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGGTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGa                                      <  2:355701/150‑1 (MQ=255)
                                                                                                   ccATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATc                                                                                              <  1:1809365/91‑1 (MQ=255)
                                                                                                   ccATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATc                                                                                              >  2:1809365/1‑91 (MQ=255)
                                                                                                   ccATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGtt                                   >  2:1112383/1‑150 (MQ=255)
                                                                                                   ccATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGtt                                   <  1:1318830/150‑1 (MQ=255)
                                                                                                   ccATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCCCGGGTTGATCCGCGAAGCGGCCGGAACAGCTAAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCGCATTGAGtt                                   >  1:298588/1‑150 (MQ=255)
                                                                                                      tCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTgc                                <  1:1112383/150‑1 (MQ=255)
                                                                                                               cAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCg                       >  1:675666/1‑150 (MQ=255)
                                                                                                                aGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCGGGTTGTGGACATTGAGTTTGGGCAAAATCGc                      >  1:1094359/1‑150 (MQ=255)
                                                                                                                  tGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACCGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGCTGTGCACCTTGAGTTTGCGCAACATCGCCg                    >  2:274059/1‑150 (MQ=255)
                                                                                                                  tGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGGGGATAACACCTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGcgg                    >  1:991990/1‑148 (MQ=255)
                                                                                                                        ggCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGt              <  2:1132514/150‑1 (MQ=255)
                                                                                                                         gCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATAGCCAGCACTGAGGACGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGAACAATCGCCGAAACggg             <  2:1529523/150‑3 (MQ=255)
                                                                                                                           aTAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTggg           >  1:999044/1‑150 (MQ=255)
                                                                                                                            tAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACAtt                                        <  2:722146/120‑1 (MQ=255)
                                                                                                                            tAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACAtt                                        >  1:722146/1‑120 (MQ=255)
                                                                                                                                   tGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAACTCGCCGAAACGTGGGCCTTGACg   <  1:881122/150‑1 (MQ=255)
                                                                                                                                   cggCCGCTTCGCGGATCAACCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGc                                             >  1:931820/2‑108 (MQ=38)
                                                                                                                                   cggCCGCTTCGCGGATCAACCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGc                                             <  2:931820/107‑1 (MQ=38)
                                                                                                                                    ggCCGCTTCGCGGATCAACCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgt                                               <  2:223019/105‑1 (MQ=38)
                                                                                                                                    ggCCGCTTCGCGGATCAACCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgt                                               >  1:223019/1‑105 (MQ=38)
                                                                                                                                    ggCCGCTTCGCGGATCAACCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTGGGCCTTGACgg  >  2:714762/1‑150 (MQ=255)
                                                                                                                                                 gTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCtt                                                           <  2:1645563/80‑1 (MQ=255)
                                                                                                                                                 gTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCtt                                                           >  1:1645563/1‑80 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
GGTCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTGGGCCTTGACGG  >  NC_002947/3053116‑3053397

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