Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 457,616 T→G L64R (CTA→CGA)  clpX → ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913457,6160TG97.6% 280.4 / ‑1.1 82L64R (CTA→CGA) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (41/39);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCTCTTTTTGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAA  >  NC_000913/457469‑457762
                                                                                                                                                   |                                                                                                                                                  
gCTCTTTTTGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGAcc                                                                                                                                                 <  2:220131/151‑1 (MQ=255)
   ctTTTTGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGcgaccgac                                                                                                                                              >  1:103814/1‑151 (MQ=255)
    tttttGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGAcg                                                                                                                                             >  1:1167104/1‑151 (MQ=255)
     ttttGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGcgaccgac                                                                                                                                              <  1:971561/149‑1 (MQ=255)
       ttGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGAcgccg                                                                                                                                          >  1:1190132/1‑151 (MQ=255)
          cGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCAt                                                                                                                                       <  1:1055436/151‑1 (MQ=255)
           ggCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATg                                                                                                                                      >  1:961758/1‑151 (MQ=255)
           ggCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATg                                                                                                                                      >  2:1061832/1‑151 (MQ=255)
            gCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGa                                                                                                                                     <  2:1169350/151‑1 (MQ=255)
             cAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGaa                                                                                                                                    <  1:875036/151‑1 (MQ=255)
              aaaaGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGaaa                                                                                                                                   <  1:1068637/151‑1 (MQ=255)
              aaaaGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGaaa                                                                                                                                   <  1:909082/151‑1 (MQ=255)
                    cAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGc                                                                                                                              >  2:483621/1‑150 (MQ=255)
                     aGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCa                                                                                                                             >  2:640603/1‑150 (MQ=255)
                     aGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCa                                                                                                                             >  2:817380/1‑150 (MQ=255)
                          gAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAaccac                                                                                                                        <  2:1028959/150‑1 (MQ=255)
                           aaGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCt                                                                                                                      <  1:591512/151‑1 (MQ=255)
                             gTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTgg                                                                                                                    <  2:153305/151‑1 (MQ=255)
                                 gcAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGAc                                                                                                                  >  2:500172/1‑149 (MQ=255)
                                   aaGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                              >  1:1200027/1‑151 (MQ=255)
                                     gCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                            <  1:899997/151‑1 (MQ=255)
                                        gATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGt                                                                                                          <  1:759564/150‑1 (MQ=255)
                                          ttGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTAt                                                                                                       >  1:126113/1‑151 (MQ=255)
                                           ttccGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                      >  2:793358/3‑151 (MQ=255)
                                            gCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                      <  2:1020696/150‑1 (MQ=255)
                                               ggTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCGCCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGc                                                                                                   <  2:912366/150‑1 (MQ=255)
                                               ggTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                  >  2:1083984/1‑151 (MQ=255)
                                                gTCCATCCGTGTATATCTGCGACGATTGTGTTGATTTATGGAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                  >  2:1208263/1‑150 (MQ=255)
                                                gTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCa                                                                                                 <  2:229198/151‑1 (MQ=255)
                                                gTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCa                                                                                                 <  2:878569/151‑1 (MQ=255)
                                                   cATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAgg                                                                                               >  1:1037902/1‑150 (MQ=255)
                                                    aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGa                                                                                              >  1:384010/1‑150 (MQ=255)
                                                     tCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAAc                                                                                            >  2:720839/1‑151 (MQ=255)
                                                        gtgtATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAgg                                                                                         >  1:161116/1‑151 (MQ=255)
                                                            atatCTGCGACGAATGCGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCg                                                                                       >  1:670464/1‑149 (MQ=255)
                                                                   cgacgaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTg                                                                              <  2:446188/151‑1 (MQ=255)
                                                                    gacgaATGTGTTGAATTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGAGACCGAAGCCGCATGTAATCCGCAACCACCTGGACGATTCCGTTATCGGCCAGGAACAGGAGAAAAAATTg                                                                              >  2:945747/1‑150 (MQ=255)
                                                                     acgaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                            >  1:928857/1‑151 (MQ=255)
                                                                     acgaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                            >  2:1133810/1‑151 (MQ=255)
                                                                       gaATGTGTTGATTTATGTAACGACATCATTCGCGAAGTGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCTCCTGGACGATTACGTTATCGGCCGGGAACAGGCGAAAAAAGTGCTgg                                                                          <  2:887084/151‑1 (MQ=255)
                                                                       gaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTg                                                                           <  2:1125557/150‑1 (MQ=255)
                                                                       gaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGCGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTgg                                                                          >  2:723716/1‑151 (MQ=255)
                                                                           gtgtTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCgg                                                                       <  2:2072/150‑1 (MQ=255)
                                                                              ttGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                   >  1:586935/1‑151 (MQ=255)
                                                                                gATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCtggcggtggcgg                                                                 >  2:30369/1‑151 (MQ=255)
                                                                                 aTTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCgg                                                                 <  2:213292/150‑1 (MQ=255)
                                                                                     aTGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGtata                                                             <  2:618220/150‑1 (MQ=255)
                                                                                     aTGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGAACGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCAAAGAACAGGCGAAAAAAGTGCTGGCGGACGCGGtcta                                                             >  1:970914/1‑150 (MQ=255)
                                                                                       gTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACa                                                           <  1:939438/150‑1 (MQ=255)
                                                                                        tAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGGCGCGGTATACaa                                                          >  2:1034364/1‑150 (MQ=255)
                                                                                        tAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGGCGCGGCATACaa                                                          >  2:626942/1‑150 (MQ=255)
                                                                                          aCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACaa                                                          >  2:539885/1‑148 (MQ=255)
                                                                                             aCATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTa                                                    <  1:587514/151‑1 (MQ=255)
                                                                                                tcatTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATAAAACCATTaaa                                                  >  1:1125644/1‑148 (MQ=255)
                                                                                                  atTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGAGACAATCAtttaaaa                                                >  2:1068841/1‑145 (MQ=255)
                                                                                                      gcgAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCTACGCCCCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCt                                           <  2:941247/151‑1 (MQ=255)
                                                                                                         aaGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCCGGAACAGGCGAAAAAAGTGCTGGCGGGCGGGGGATACAACCATTACAAAccgctg                                          >  2:147429/1‑144 (MQ=255)
                                                                                                              aTTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACg                                   <  1:1132001/151‑1 (MQ=255)
                                                                                                               ttAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAACCATTTCAAACGTCTGCGCCAc                                    >  2:241148/1‑149 (MQ=255)
                                                                                                                 aaaGAAGTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGc                                 <  2:735593/150‑1 (MQ=255)
                                                                                                                  aagaagTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCg                                >  1:107613/1‑150 (MQ=255)
                                                                                                                     aagTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTGTACAACCATTACAAAAGTCTGCGCAACGGCGATa                             >  2:217858/1‑150 (MQ=255)
                                                                                                                      agTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAACCATTACAAACGTCTGCGAAACGGc                                 >  1:971962/1‑145 (MQ=255)
                                                                                                                       gTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                          <  1:1061902/151‑1 (MQ=255)
                                                                                                                       gTTGCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAACCATTACAAACGTCTGCGCAACGGCGATAcc                           >  2:138296/1‑150 (MQ=255)
                                                                                                                          gCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGc                        <  2:1180453/150‑1 (MQ=255)
                                                                                                                          gCACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGc                        >  1:941122/1‑150 (MQ=255)
                                                                                                                           cACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                      <  2:561166/151‑1 (MQ=255)
                                                                                                                           cACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCAATACAAACGTCTGCGCAACGGCGATACCa                          >  1:1160688/1‑147 (MQ=255)
                                                                                                                           cACCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAAACCTTACAAACGTCTGCGCAACGGCGATAACAGCa                       >  2:1047496/1‑150 (MQ=255)
                                                                                                                            aCCGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                      >  1:200885/1‑150 (MQ=255)
                                                                                                                             ccGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATg                    <  1:793399/151‑1 (MQ=255)
                                                                                                                             ccGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATg                    <  2:283730/151‑1 (MQ=255)
                                                                                                                             ccGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGGATACAACCATTACAAACGTTTGCGCAACGGCGATACCAGc                        >  2:109256/1‑147 (MQ=255)
                                                                                                                              cGCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGCAAACAACCACTTCAAACATCTGCGCAACGGCGAGAACAGGCATgg                   >  2:68384/1‑151 (MQ=255)
                                                                                                                               gCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATgg                   >  2:502011/1‑150 (MQ=255)
                                                                                                                               gCATCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGc                  <  1:930093/151‑1 (MQ=255)
                                                                                                                                 aTCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGt                <  1:497091/151‑1 (MQ=255)
                                                                                                                                 aTCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGt                <  2:928796/151‑1 (MQ=255)
                                                                                                                                 aTCGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGt                <  1:719908/151‑1 (MQ=255)
                                                                                                                                   cGTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCg              <  1:77044/151‑1 (MQ=255)
                                                                                                                                    gTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGa             <  2:380593/151‑1 (MQ=255)
                                                                                                                                    gTGAACGCAGTGCGCGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGa             <  2:182691/151‑1 (MQ=255)
                                                                                                                                               gcgcGACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaaa  <  2:995005/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCTTTTTGCGGCAAAAGCCAGCATGAAGTGCGCAAGCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAA  >  NC_000913/457469‑457762

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