Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,659,949 C→G 9.1% intergenic (‑41/+34) PP_3223 ← / ← PP_3224 ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,659,9490CG 9.1% 294.6 / 16.6 99intergenic (‑41/+34)PP_3223/PP_3224ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (45/45);  new base G (6/3);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GAAGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659801‑3660098
                                                                                                                                                    |                                                                                                                                                     
gAAGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTTATCCTTCAATTAATGGGGCTCTCCTGTTCCGGCCa                                                                                                                                                      <  1:190951/150‑1 (MQ=255)
 aaGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGAATGGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAg                                                                                                                                                     >  1:726160/1‑150 (MQ=255)
   ggCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCGTTCAATGGATGGGGCGCTCCTGTTCCGGCCAGtt                                                                                                                                                   >  2:1591050/1‑150 (MQ=255)
   ggCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGtt                                                                                                                                                   >  1:1748287/1‑150 (MQ=255)
      gATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGc                                                                                                                                                >  2:2223212/1‑150 (MQ=255)
        tGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCg                                                                                                                                              >  1:1756841/1‑150 (MQ=255)
          ttGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGa                                                                                                                                            >  2:1944232/1‑150 (MQ=255)
           tGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGaa                                                                                                                                           <  1:1016362/150‑1 (MQ=255)
                     ggtggtggGGGTGTGGTCGGACGGGGAGAGCGCGGTGCGGCCCAGGGACGGGGCGCCGATAGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAg                                                                                                                                 <  2:500641/148‑1 (MQ=255)
                     gCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAg                                                                                                                                 >  1:72308/1‑150 (MQ=255)
                            gAGGTGTTGTCGGCCGGCGAGCACGCAGTGGGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          <  2:2133888/150‑1 (MQ=255)
                            gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          <  1:1504871/150‑1 (MQ=255)
                            gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAACAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAg                                                                                                                                 >  1:244224/1‑143 (MQ=255)
                            gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAACAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAg                                                                                                                                 <  2:244224/143‑1 (MQ=255)
                               gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCCGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  1:1591050/150‑1 (MQ=255)
                               gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  1:2223212/150‑1 (MQ=255)
                               gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCCCTTTATCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  1:555265/150‑1 (MQ=255)
                                tgttgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGag                                                                                                                      <  2:1748287/150‑1 (MQ=255)
                                       ggccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTcc                                                                                                               <  1:2152821/150‑1 (MQ=255)
                                         ccggcGAGTAAGCGGTGAGGCGCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  2:2398826/150‑1 (MQ=255)
                                         ccggcGAGCAAGCGGTGAGGCCCAGGGGCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCAAGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGCCAGTCCTg                                                                                                             <  1:1860484/150‑1 (MQ=255)
                                         ccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGTAATCTTTCCTg                                                                                                             <  1:211346/150‑1 (MQ=255)
                                         ccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  1:917869/150‑1 (MQ=255)
                                          cggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            >  2:1794497/1‑150 (MQ=255)
                                            gcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACa                                                                                                                                         >  1:628706/1‑119 (MQ=255)
                                            gcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACa                                                                                                                                         <  2:628706/119‑1 (MQ=255)
                                                 cAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGa                                                                                                     >  1:592875/1‑150 (MQ=255)
                                                  aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTg                                                                                                                         >  2:271236/1‑129 (MQ=255)
                                                  aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTg                                                                                                                         <  1:271236/129‑1 (MQ=255)
                                                  aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    <  1:180673/150‑1 (MQ=255)
                                                      ggTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGccc                                                                                                >  1:2073533/1‑150 (MQ=255)
                                                      ggTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGccc                                                                                                >  2:57910/1‑150 (MQ=255)
                                                        tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  2:352255/135‑1 (MQ=255)
                                                        tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             >  1:352255/1‑135 (MQ=255)
                                                             cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAAt                                                                                                        <  2:361680/135‑1 (MQ=255)
                                                             cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAAt                                                                                                        >  1:361680/1‑135 (MQ=255)
                                                                agggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCa                                                                                      >  2:478396/1‑150 (MQ=255)
                                                                    ccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTg                                                                                  >  2:2404848/1‑150 (MQ=255)
                                                                             cGATCGCGAGCAGTTGGTGCCGTGGGGGCATGGGATTGATCCTTCAATTAAGGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACa                                                                                                                                         >  2:1843901/1‑86 (MQ=25)
                                                                              gATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAgc                                                                        <  1:2404848/150‑1 (MQ=255)
                                                                                 cgcgAGCAGTTGGTGTAGTGGGGGCATGGGATTGATCCTTCAATTAATGGGGCTCGCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTCGGAGTCAGTCCTGGTAATCGCGGCCCTCGCGATCCAATTGCCGAATCAGCGCa                                                                     <  2:27075/150‑1 (MQ=255)
                                                                                      gCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCa                                                                <  1:2084131/150‑1 (MQ=255)
                                                                                        aGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              >  1:29366/1‑150 (MQ=255)
                                                                                          ttGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGc                                                            >  2:2136916/1‑150 (MQ=255)
                                                                                           tGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       >  1:1755049/1‑90 (MQ=255)
                                                                                           tGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  2:1755049/90‑1 (MQ=255)
                                                                                              tGCAGTGGGGGCATGGTATCGATCCATCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGtt                                                        <  1:2348900/150‑1 (MQ=255)
                                                                                              tGCAGTGGGGGCATGGGATTGGTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTACTCGAGGCCCTCGCGATCCAACTGCCGAATTAGCGCGTTCTAATGCCGtt                                                        >  2:23121/1‑150 (MQ=255)
                                                                                              tGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    <  1:1743267/106‑1 (MQ=255)
                                                                                              tGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  2:1743267/1‑106 (MQ=255)
                                                                                               gCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       <  1:1668206/150‑1 (MQ=255)
                                                                                                 aGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGa                                                     <  1:1944232/150‑1 (MQ=255)
                                                                                                    gggggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGc                                                                                                  <  2:2084804/102‑1 (MQ=255)
                                                                                                    gggggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGc                                                                                                  >  1:2084804/1‑102 (MQ=255)
                                                                                                     ggggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATAcc                                                 >  1:733394/1‑150 (MQ=255)
                                                                                                        gCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCAGAAAAAGCGCATCCCCATACCCATGGATAccggc                                              >  2:1959341/1‑147 (MQ=255)
                                                                                                           tGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTcccc                                           <  2:29366/150‑1 (MQ=255)
                                                                                                            gggATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTccccc                                          >  2:2401918/1‑150 (MQ=255)
                                                                                                              gATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTg                                        <  2:2073533/150‑1 (MQ=255)
                                                                                                               aTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCa                                                                          >  1:302913/1‑115 (MQ=255)
                                                                                                               aTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCa                                                                          <  2:302913/115‑1 (MQ=255)
                                                                                                               aTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTgg                                       >  2:1349197/1‑150 (MQ=255)
                                                                                                                  gATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAgga                                                                                                                                      <  2:559417/52‑2 (MQ=14)
                                                                                                                  gATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAgga                                                                                                                                      >  1:559417/1‑51 (MQ=17)
                                                                                                                     ccTTCAATGAATGGGGCGCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGACGTCAGTCCTGGTAATCTAGGCACTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 <  2:120257/150‑1 (MQ=255)
                                                                                                                     ccTGCAATGAATGGGGCTCTCCTGTTCGGGCCAGTTCGCCGTACCAGGCAGTAGTGGCGTTGGGTGGAGTCCTGGCAATCGAGGCCCTCGCGATCGGACTCAACAATCAGCGCATCCCAACGCCGTTGGATACCTCCCCCCTGGCCATCg                                 <  2:119471/150‑1 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGTAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:719479/1‑150 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     <  2:1503570/96‑1 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     >  1:1503570/1‑96 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  2:1759826/1‑150 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:2236400/1‑150 (MQ=255)
                                                                                                                       ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAACAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:512978/1‑150 (MQ=255)
                                                                                                                        tCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGaa                              >  2:2127827/1‑150 (MQ=255)
                                                                                                                          aatgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        <  1:706690/90‑1 (MQ=255)
                                                                                                                          aatgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        >  2:706690/1‑90 (MQ=255)
                                                                                                                           atgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGag                                                                                                                      <  2:95113/59‑1 (MQ=255)
                                                                                                                           atgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGag                                                                                                                      >  1:95113/1‑59 (MQ=255)
                                                                                                                            tgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAAcc                          <  1:2136916/150‑1 (MQ=255)
                                                                                                                                tgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCCCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCtt                      <  2:726160/150‑1 (MQ=255)
                                                                                                                                  gggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCAACGGCAAACAGCTTgg                    >  2:2089528/1‑150 (MQ=255)
                                                                                                                                   ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       <  1:721910/114‑1 (MQ=255)
                                                                                                                                   ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       >  2:721910/1‑114 (MQ=255)
                                                                                                                                    ttttgccTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGAc                  >  2:1770964/6‑150 (MQ=255)
                                                                                                                                    gCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAt                                                                            <  2:1543591/92‑1 (MQ=255)
                                                                                                                                    gCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAt                                                                            >  1:1543591/1‑92 (MQ=255)
                                                                                                                                       ctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGc               <  1:2401918/150‑1 (MQ=255)
                                                                                                                                        tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGGGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGCCTGCt              >  1:1293916/1‑150 (MQ=255)
                                                                                                                                        tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGa                   <  2:129429/145‑1 (MQ=255)
                                                                                                                                        tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGa                   >  1:129429/1‑145 (MQ=255)
                                                                                                                                         ccTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGc               <  2:2028619/148‑1 (MQ=255)
                                                                                                                                         ccTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGc               >  1:2028619/1‑148 (MQ=255)
                                                                                                                                         ccTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     <  2:454764/78‑1 (MQ=255)
                                                                                                                                         ccTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     >  1:454764/1‑78 (MQ=255)
                                                                                                                                          cTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGa                                                                                          >  1:1457832/1‑72 (MQ=25)
                                                                                                                                          cTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAGAACCGCTTGGACTGCTCa            >  1:2203539/1‑150 (MQ=255)
                                                                                                                                                 ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGa                                                     >  2:313869/1‑102 (MQ=255)
                                                                                                                                                 ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGa                                                     <  1:313869/102‑1 (MQ=255)
                                                                                                                                                  gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACGTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    >  2:437774/1‑150 (MQ=255)
                                                                                                                                                  gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTACCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    <  2:1756841/150‑1 (MQ=255)
                                                                                                                                                    cAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGa  >  1:1608509/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAGGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659801‑3660098

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