Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,000,796 T→G 5.8% I110L (ATC→CTC)  ycaC‑IV ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,000,7960TG 5.8% 1166.7 / 38.9 343I110L (ATC→CTC) ycaC‑IVputative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (159/161);  new base G (9/11);  total (169/174)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

TGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGCTCGAAGC  >  NC_002947/6000647‑6000945
                                                                                                                                                     |                                                                                                                                                     
tgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatgat                                                                                                                                                       <  1:3143527/150‑1 (MQ=255)
 gcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATc                                                                                                                                                      >  2:1037219/1‑150 (MQ=255)
        tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                                                                                                               >  2:8025814/1‑150 (MQ=255)
             cGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTg                                                                                                                                          <  2:697428/150‑1 (MQ=255)
                tCGTGGGCGGCGTCGCTGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCa                                                                                                                                       <  1:2390732/150‑1 (MQ=255)
                tCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCa                                                                                                                                       >  2:447252/1‑150 (MQ=255)
                   tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTg                                                                                                                                    <  2:1569896/150‑1 (MQ=255)
                   tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGAACAGTTGCTTCTTGCCAGTg                                                                                                                                    <  2:3955410/150‑1 (MQ=255)
                   tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGAACAGTTGCTTCTTGCCAGTg                                                                                                                                    <  2:3180350/150‑1 (MQ=255)
                           gTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCAc                                                                                                                            >  2:141800/1‑150 (MQ=255)
                                   tCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGa                                                                                                                    <  2:8898191/150‑1 (MQ=255)
                                    cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                                                                                   <  2:858085/150‑1 (MQ=255)
                                    cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                                                                                   <  1:8406592/150‑1 (MQ=255)
                                    cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                                                                                   <  1:7567882/150‑1 (MQ=255)
                                      tGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGa                                                                                                                 >  1:5004311/1‑150 (MQ=255)
                                      tGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGa                                                                                                                 >  2:6485694/1‑150 (MQ=255)
                                       ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                                                                                                                >  1:1030459/5‑150 (MQ=255)
                                       ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                                <  1:3679610/150‑1 (MQ=255)
                                       ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                                <  1:2050050/150‑1 (MQ=255)
                                          gTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTc                                                                                                             <  2:3477271/150‑1 (MQ=255)
                                          attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                             >  1:2147562/2‑150 (MQ=255)
                                            tGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCtt                                                                                                           <  2:1844390/150‑1 (MQ=255)
                                             gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                          >  2:4220440/1‑150 (MQ=255)
                                             gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                          <  2:9110949/150‑1 (MQ=255)
                                              aaGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCg                                                                                                         >  1:2972637/1‑150 (MQ=255)
                                              aaGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCg                                                                                                         <  2:2645127/150‑1 (MQ=255)
                                               aGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGt                                                                                                        >  1:6256065/1‑150 (MQ=255)
                                                ggTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGtt                                                                                                       <  2:24068/150‑1 (MQ=255)
                                                 gTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTg                                                                                                      <  2:6328218/150‑1 (MQ=255)
                                                  tACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGt                                                                                                     <  2:2796374/150‑1 (MQ=255)
                                                     cGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTccc                                                                                                  >  2:7095481/1‑150 (MQ=255)
                                                        aGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAgg                                                                                               >  2:4946160/1‑150 (MQ=255)
                                                         ggCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGc                                                                                              <  2:407993/150‑1 (MQ=255)
                                                         cgCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGc                                                                                              <  1:6730216/149‑1 (MQ=255)
                                                            aTCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGtt                                                                                           <  2:2903351/150‑1 (MQ=255)
                                                            aTCAGTGACCACGAACACGTCGAACTCTTGTTCGAGGGCCGAGAGGGACTTTAAAGCTACGCTCCCTTGCGTCACAACTCCGGCGATGATGAGGTGGGTCTTGCGACGTGACTTCTCTTCCTTGACGTAGTCTTCGTTGTCCCAGGCGtt                                                                                           <  2:7644962/150‑1 (MQ=255)
                                                            aTCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                                                                                           >  1:7303127/1‑150 (MQ=255)
                                                               aGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt                                                                                        <  1:1429672/150‑1 (MQ=255)
                                                                  gACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTg                                                                                     <  1:3470463/150‑1 (MQ=255)
                                                                   aCCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTgg                                                                                    <  2:3715692/150‑1 (MQ=255)
                                                                      aCGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                                 >  1:1516623/1‑150 (MQ=255)
                                                                       cGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAg                                                                                >  1:5922893/1‑150 (MQ=255)
                                                                        gcacacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgg                                                                               <  2:6758949/148‑1 (MQ=255)
                                                                         aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg                                                                              <  1:8071427/150‑1 (MQ=255)
                                                                         aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg                                                                              <  1:1570470/150‑1 (MQ=255)
                                                                          ccctggtgAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                             <  1:6979276/143‑1 (MQ=255)
                                                                          acacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                             >  1:5746350/1‑150 (MQ=255)
                                                                           cacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcg                                                                            >  1:8187957/1‑150 (MQ=255)
                                                                              gTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                         >  2:1296680/1‑150 (MQ=255)
                                                                               tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                          >  1:1134809/1‑116 (MQ=255)
                                                                               tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                        >  2:569680/1‑150 (MQ=255)
                                                                               tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGAGGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                          <  2:1134809/116‑1 (MQ=255)
                                                                                 gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                                      <  2:2489128/150‑1 (MQ=255)
                                                                                  aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg                                                                     >  2:8837554/1‑150 (MQ=255)
                                                                                  aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg                                                                     >  2:6238579/1‑150 (MQ=255)
                                                                                    ctctTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTa                                                                   <  1:9133913/150‑1 (MQ=255)
                                                                                     tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAg                                                                                                               >  2:5842641/1‑105 (MQ=255)
                                                                                     tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGt                                                                                                              <  1:5492511/106‑1 (MQ=255)
                                                                                     tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                                  <  1:5899663/150‑1 (MQ=255)
                                                                                     acttcttCCATGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                                  <  2:9145419/149‑1 (MQ=255)
                                                                                      cttcttCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACg                                                                 >  1:1741333/1‑150 (MQ=255)
                                                                                      cttcttCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACg                                                                 <  2:2142047/150‑1 (MQ=255)
                                                                                       ttcttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                <  2:3312852/150‑1 (MQ=255)
                                                                                          ttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGc                                                             >  1:6093658/1‑150 (MQ=255)
                                                                                            ccAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGt                                                           >  2:5269114/1‑150 (MQ=255)
                                                                                             cAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTc                                                          >  2:4589065/1‑150 (MQ=255)
                                                                                                    gACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAAc                                                   >  1:4634146/1‑150 (MQ=255)
                                                                                                     aCAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACa                                                  <  2:6016204/150‑1 (MQ=255)
                                                                                                      cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                                 <  2:6283325/150‑1 (MQ=255)
                                                                                                      cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                                 <  2:3147801/150‑1 (MQ=255)
                                                                                                      cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                                 <  2:167789/150‑1 (MQ=255)
                                                                                                        gggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGTTCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAggg                                               <  1:7702963/150‑1 (MQ=255)
                                                                                                        gggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAggg                                               <  1:8320398/150‑1 (MQ=255)
                                                                                                         ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                              <  1:398151/150‑1 (MQ=255)
                                                                                                         ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                              <  1:2023206/150‑1 (MQ=255)
                                                                                                          gCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCt                                             >  1:4783930/1‑150 (MQ=255)
                                                                                                              ggAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCa                                         >  1:3347864/1‑150 (MQ=255)
                                                                                                                   gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCa                                                                                                                             >  1:8804389/1‑61 (MQ=255)
                                                                                                                   gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCt                                    <  1:1218996/150‑1 (MQ=255)
                                                                                                                     tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                                  <  2:7320625/150‑1 (MQ=255)
                                                                                                                     tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                                  >  1:4335519/1‑150 (MQ=255)
                                                                                                                     tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                                  <  1:4485728/150‑1 (MQ=255)
                                                                                                                     tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                                  >  1:5787113/1‑150 (MQ=255)
                                                                                                                        gCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCAc                               <  2:8356151/150‑1 (MQ=255)
                                                                                                                            acTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGt                           >  1:9071771/1‑150 (MQ=255)
                                                                                                                            acTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTGTCAGTTCTGGCACCAGt                           >  2:5052685/1‑150 (MQ=255)
                                                                                                                              ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg                         <  1:4102445/150‑1 (MQ=255)
                                                                                                                              ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg                         <  1:4936672/150‑1 (MQ=255)
                                                                                                                              ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg                         <  2:4390050/150‑1 (MQ=255)
                                                                                                                                cAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGc                       <  2:5536943/150‑1 (MQ=255)
                                                                                                                                   tCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGt                    >  1:5568352/1‑150 (MQ=255)
                                                                                                                                      caACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgg                 <  2:4983999/150‑1 (MQ=255)
                                                                                                                                       aaCGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTggg                >  1:6421483/1‑150 (MQ=255)
                                                                                                                                        aCGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgggg               >  1:868838/1‑150 (MQ=255)
                                                                                                                                          gCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGcc             >  1:798168/1‑150 (MQ=255)
                                                                                                                                              gCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGc         <  2:4655068/150‑1 (MQ=255)
                                                                                                                                              gCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGc         <  2:3660960/150‑1 (MQ=255)
                                                                                                                                               cGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGCt        >  1:7037881/1‑150 (MQ=255)
                                                                                                                                                gatgatCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGCTc       <  2:2192694/150‑1 (MQ=255)
                                                                                                                                                     tCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGCTCGAAGc  >  2:8986580/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTGCTCGAAGC  >  NC_002947/6000647‑6000945

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