Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,053,266 T→C 9.5% intergenic (+24/+22) 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,2660TC 9.5% 223.9 / 12.2 84intergenic (+24/+22)PP_2664/agmRtwo‑component system sensor histidine kinase/response regulator/glycerol metabolism activator
Reads supporting (aligned to +/- strand):  ref base T (34/42);  new base C (4/4);  total (38/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 = 2.66e-01

TCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTGGGCCTTGACGG  >  NC_002947/3053118‑3053397
                                                                                                                                                    |                                                                                                                                   
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                                                                                                                          >  1:1590215/1‑150 (MQ=255)
          cAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCgg                                                                                                                          >  2:145912/1‑150 (MQ=255)
           aaCGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGaa                                                                                                                              <  2:1334342/145‑1 (MQ=255)
           aaCGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGaa                                                                                                                              >  1:1334342/1‑145 (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:1066865/150‑1 (MQ=255)
                                 gcgAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGa                                                                                                   >  1:71744/1‑150 (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                                                                                      <  2:54487/150‑1 (MQ=255)
                                              cacGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAAt                                                                                      >  1:1050514/1‑150 (MQ=255)
                                               acGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAg                                                                                            <  1:1544315/143‑1 (MQ=255)
                                               acGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAg                                                                                            >  2:1544315/1‑143 (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                                                                                     <  1:1618848/130‑1 (MQ=255)
                                                                   aaGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATc                                                                                     >  2:1618848/1‑130 (MQ=255)
                                                                     gACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGc                                                                                           >  2:1136532/1‑122 (MQ=255)
                                                                     gACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGc                                                                                           <  1:1136532/122‑1 (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:1492755/150‑1 (MQ=255)
                                                                                   aaGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGtt                                                 <  2:1647187/150‑1 (MQ=255)
                                                                                    aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCa                                                                           <  1:1239257/123‑1 (MQ=255)
                                                                                    aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCa                                                                           >  2:1239257/1‑123 (MQ=255)
                                                                                    aGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGAAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtg                                                >  2:1357507/1‑150 (MQ=255)
                                                                                     gACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgt                                               <  1:1064712/150‑1 (MQ=255)
                                                                                      aCGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGAAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTtgtg                                              <  1:1357507/150‑1 (MQ=255)
                                                                                       cGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGATCAACCCGCGAAGCGGccaac                                                                                                                     <  1:460120/78‑4 (MQ=21)
                                                                                       cGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGATCAACCCGCGAAGCGGccaac                                                                                                                     >  2:460120/1‑75 (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)
                                                                                                                                                    |                                                                                                                                   
TCAACTACAGCAACGATCTCAAGCAGCAGATCCGCGAACTGGGTCACACGCTGATGCACAAGCCCGTAAGACCGATGAAACTCAAGACGGCGATGAGCCATCTGCTGGCCAGTGGCTTGGCATAGCGGTTGGCCGCTTCGCGGGTTGATCCGCGAAGCGGCCGGAACAGCTCAACGGCGCAGATAATCAGCAAAATCGATATCGCCAGCACTGAGGATGGCTTGTACCCGGTTGTGCACATTGAGTTTGCGCAAAATCGCCGAAACGTGGGCCTTGACGG  >  NC_002947/3053118‑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: