Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,659,952 T→C 9.4% intergenic (‑44/+31) 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,9520TC 9.4% 271.3 / 13.7 96intergenic (‑44/+31)PP_3223/PP_3224ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base C (6/3);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659804‑3660098
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                                 <  2:244224/143‑1 (MQ=255)
                         gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAACAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAg                                                                                                                                 >  1:244224/1‑143 (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                                                                                                >  2:57910/1‑150 (MQ=255)
                                                   ggTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGccc                                                                                                >  1:2073533/1‑150 (MQ=255)
                                                     tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             >  1:352255/1‑135 (MQ=255)
                                                     tGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTg                                                                                                             <  2:352255/135‑1 (MQ=255)
                                                          cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAAt                                                                                                        >  1:361680/1‑135 (MQ=255)
                                                          cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAAt                                                                                                        <  2:361680/135‑1 (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                                                                                                                       <  2:1755049/90‑1 (MQ=255)
                                                                                        tGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       >  1:1755049/1‑90 (MQ=255)
                                                                                           tGCAGTGGGGGCATGGTATCGATCCATCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGtt                                                        <  1:2348900/150‑1 (MQ=255)
                                                                                           tGCAGTGGGGGCATGGGATTGGTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTACTCGAGGCCCTCGCGATCCAACTGCCGAATTAGCGCGTTCTAATGCCGtt                                                        >  2:23121/1‑150 (MQ=255)
                                                                                           tGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  2:1743267/1‑106 (MQ=255)
                                                                                           tGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    <  1:1743267/106‑1 (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                                                                          <  2:302913/115‑1 (MQ=255)
                                                                                                            aTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCa                                                                          >  1:302913/1‑115 (MQ=255)
                                                                                                            aTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTgg                                       >  2:1349197/1‑150 (MQ=255)
                                                                                                               gATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAgga                                                                                                                                      >  1:559417/1‑51 (MQ=17)
                                                                                                               gATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAgga                                                                                                                                      <  2:559417/52‑2 (MQ=14)
                                                                                                                  ccTTCAATGAATGGGGCGCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGACGTCAGTCCTGGTAATCTAGGCACTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 <  2:120257/150‑1 (MQ=255)
                                                                                                                  ccTGCAATGAATGGGGCTCTCCTGTTCGGGCCAGTTCGCCGTACCAGGCAGTAGTGGCGTTGGGTGGAGTCCTGGCAATCGAGGCCCTCGCGATCGGACTCAACAATCAGCGCATCCCAACGCCGTTGGATACCTCCCCCCTGGCCATCg                                 <  2:119471/150‑1 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGTAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:719479/1‑150 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     >  1:1503570/1‑96 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCaa                                                                                     <  2:1503570/96‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659804‑3660098

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: