Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,709 T→A 7.3% R597R (CGT→CGA rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,7090TA 7.3% 768.2 / 18.6 268R597R (CGT→CGArsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (86/150);  new base A (17/6);  total (106/162)
Fisher's exact test for biased strand distribution p-value = 6.42e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

GCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAAC  >  NC_000913/1708560‑1708857
                                                                                                                                                     |                                                                                                                                                    
    caACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                                                                                                                  <  2:1643174/150‑1 (MQ=255)
     aaCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTAGATCCGCGCAAAGCGGAAGTCGAATCGTCTATTGACCGTGCaaa                                                                                                                                                 >  1:222618/1‑150 (MQ=17)
     aaCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                                                                                 <  1:3311640/150‑1 (MQ=255)
      aCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                                                <  1:560896/150‑1 (MQ=255)
         gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                                                                                                             <  1:1468258/150‑1 (MQ=255)
         gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgca                                                                                                                                             >  1:2441069/1‑149 (MQ=255)
             cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                                                                                                         <  1:463942/150‑1 (MQ=255)
             cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                                                                                                         >  2:2749210/1‑150 (MQ=255)
             cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                                                                                                         >  1:2150981/1‑150 (MQ=255)
             cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgcgc                                                                                                                                          >  1:617397/1‑149 (MQ=255)
              tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCaa                                                                                                                                        >  2:3958266/1‑150 (MQ=255)
                  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                                    <  2:934848/150‑1 (MQ=255)
                  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                                    <  2:1735840/150‑1 (MQ=255)
                   cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCACAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                                                                                                   <  1:1433449/150‑1 (MQ=255)
                    aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGa                                                                                                                                  >  2:2253731/1‑150 (MQ=255)
                     cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                                                >  1:3190732/1‑135 (MQ=255)
                      gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                                <  2:1144067/150‑1 (MQ=255)
                       tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACTGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                                                                               <  1:1164940/150‑1 (MQ=14)
                         aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                             <  1:3893695/150‑1 (MQ=255)
                         aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                             <  1:3153848/150‑1 (MQ=255)
                             tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                         <  1:2315263/150‑1 (MQ=255)
                             tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                         <  1:2315477/150‑1 (MQ=255)
                             tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                         <  1:1265861/150‑1 (MQ=255)
                              gCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAccg                                                                                                                        <  2:876686/150‑3 (MQ=255)
                               ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                                                                                                       <  2:752066/150‑1 (MQ=255)
                               ccGGTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGCGCGCTCCCCGCCCAGCGGCCGGCGAAGCCGCTCGGCCCCGGGCCACAGCGCGCACGCGGGCCCAGCACCAgc                                                                                                                       >  2:2916739/1‑149 (MQ=9)
                                 gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCGACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCa                                                                                                                                                  <  2:460403/121‑1 (MQ=255)
                                  ttGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                                                                                    <  2:659210/150‑1 (MQ=255)
                                     aaGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATg                                                                                                                 <  1:1025538/150‑1 (MQ=255)
                                       gcagcTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                               <  2:2375860/150‑1 (MQ=255)
                                       gcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTCGGGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                               >  2:1502332/1‑150 (MQ=255)
                                         agcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                                                                                                             <  1:878621/150‑1 (MQ=255)
                                          gcTATCGCCCTCGCCAAAGCTCGCATGCTGGAACAGCAACAGGATAATGCCGACGCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTGGAACAGCAACAGGCTAATGCGGag                                                                                                            <  2:2938782/150‑2 (MQ=11)
                                            tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAccc                                                                                                          >  2:3600404/1‑150 (MQ=255)
                                              tCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCag                                                                                                        <  1:1324846/150‑1 (MQ=255)
                                                   cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCTATCGGGAACAGCAACAGGCTAATGCGGAGCCAgaagaa                                                                                                   >  2:2406735/1‑150 (MQ=255)
                                                   cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAgaagaa                                                                                                   >  2:1436464/1‑150 (MQ=255)
                                                   cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGATACTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAgaagaa                                                                                                   <  2:3806941/150‑1 (MQ=255)
                                                   cgcgACAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAgaagaa                                                                                                   <  2:1254876/150‑1 (MQ=255)
                                                        cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                              <  1:3577262/150‑1 (MQ=255)
                                                         aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGAAACATGATAAATCCGAACCAGAATAACAggaa                                                                                             >  1:210069/1‑148 (MQ=11)
                                                          aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCAAGACGAACAGGTCg                                                                                            >  1:1629984/1‑150 (MQ=255)
                                                            gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                                                                          <  1:1864195/150‑1 (MQ=255)
                                                            gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                                                                          <  1:427885/150‑1 (MQ=255)
                                                             ctcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                         <  2:1287630/148‑1 (MQ=255)
                                                             cgctcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACTTCATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCTCTATTGCCCGTGCCAAAGCGCGCATGCTGGTACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                         <  2:1851820/150‑1 (MQ=14)
                                                             cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                         <  2:3289232/150‑1 (MQ=255)
                                                             cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                         <  2:2232523/150‑1 (MQ=255)
                                                                  aaGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCgcgc                                                                                    >  1:2722084/1‑150 (MQ=255)
                                                                  aaGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCgcgc                                                                                    <  2:1385956/150‑1 (MQ=255)
                                                                   aGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACa                                                                                                 >  2:3043011/1‑136 (MQ=12)
                                                                        gAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAgct                                                                              >  2:1033748/1‑149 (MQ=255)
                                                                        gAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcg                                                                              >  2:1946453/1‑149 (MQ=255)
                                                                        gAACAGCAACAGGATAATGAGGAACCAGAACCCCAGGTCGATCCGCGCACAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                              <  2:2658209/150‑1 (MQ=14)
                                                                         aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTGGAGCCGCGCAGAgcgg                                                                             >  2:3400573/1‑148 (MQ=255)
                                                                         aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                             <  1:4032590/150‑1 (MQ=255)
                                                                          aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCTGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTCGATCCGCGCAAAgccgc                                                                            >  1:3667262/1‑150 (MQ=255)
                                                                           cAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTAGATCCGCGCAAAgcggcc                                                                           >  1:2596581/1‑150 (MQ=255)
                                                                           cAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTCATGCGGAACCAGAAGAACAGGTGGATCCGCGCAAAgcggcc                                                                           >  2:3232625/1‑150 (MQ=255)
                                                                             tcatcAGGCTAATGCGGAACCAGAACAACATGTCGATCCGCTCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt                                                                         <  2:3269301/146‑1 (MQ=255)
                                                                              cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                        >  1:2889841/1‑150 (MQ=255)
                                                                                     cTAATGCGGAACCAGAACAACAGGACGATCCGCACAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                 <  1:2361349/150‑1 (MQ=255)
                                                                                      tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                <  1:3022670/150‑1 (MQ=255)
                                                                                      tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                <  2:522943/150‑1 (MQ=255)
                                                                                          gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAtt                                                            >  2:2602678/1‑150 (MQ=255)
                                                                                             gAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGcc                                                         >  2:1704207/1‑150 (MQ=255)
                                                                                             gAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                         <  2:704301/150‑1 (MQ=255)
                                                                                             gAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTCGCAGCGGCGATTGca                                                         >  2:2042661/1‑149 (MQ=11)
                                                                                             gAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGcc                                                         >  1:2770239/1‑150 (MQ=255)
                                                                                                ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGACGCGGCTATTGCCCga                                                      >  2:3678517/1‑149 (MQ=255)
                                                                                                 cAGAACCCCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                     <  2:3629410/150‑1 (MQ=255)
                                                                                                  aGAAAAACTGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                    <  2:964510/150‑1 (MQ=255)
                                                                                                    aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                  <  2:246053/150‑1 (MQ=255)
                                                                                                     acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                 <  1:1720716/150‑1 (MQ=255)
                                                                                                                aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGc                                      >  1:716880/1‑150 (MQ=255)
                                                                                                                aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGc                                      >  1:3543953/1‑150 (MQ=255)
                                                                                                                 tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                     <  1:3814079/150‑1 (MQ=255)
                                                                                                                  cctcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTg                                    <  2:318151/147‑1 (MQ=255)
                                                                                                                  ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGACAAGGTCGAGCCGCGCGAAGCCGCCGTCGAAGCCGCTCTTGCCCGAGCCACTGCGCGCAAGCTg                                    >  2:3973953/1‑150 (MQ=11)
                                                                                                                  ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTg                                    <  1:2465041/150‑1 (MQ=255)
                                                                                                                   cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                   <  1:1993686/150‑1 (MQ=255)
                                                                                                                   cgcgCAAAGCCGCCGTCGAAGCCGCTATGGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                   <  1:423494/150‑1 (MQ=255)
                                                                                                                     cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                 <  1:4143257/150‑1 (MQ=255)
                                                                                                                        aaaGCCGCCGTAGAAGCCGCTATTGACCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCTCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAg                              <  2:4042026/150‑1 (MQ=255)
                                                                                                                              gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAg                        >  2:3682505/1‑150 (MQ=255)
                                                                                                                              gccgTCGAAGCCGCTATTGCCCGTGACAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAg                        <  1:1066298/150‑1 (MQ=255)
                                                                                                                                cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCGGGAACAGCAACAGGc                      >  1:945579/1‑150 (MQ=255)
                                                                                                                                 gTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                     <  2:2121221/150‑1 (MQ=255)
                                                                                                                                       gCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGACAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcg               >  2:339889/1‑149 (MQ=255)
                                                                                                                                       gCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGACGTCGAAGCCGCTATTGCCAGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGcc               >  1:189410/1‑150 (MQ=255)
                                                                                                                                        ccGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCg              <  2:842670/150‑1 (MQ=255)
                                                                                                                                         cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa             >  1:1755080/1‑150 (MQ=255)
                                                                                                                                          gCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcggta            >  2:3203193/1‑146 (MQ=11)
                                                                                                                                            tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc          <  2:3729963/150‑1 (MQ=255)
                                                                                                                                            tATGGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc          <  2:2082602/150‑1 (MQ=255)
                                                                                                                                             aTTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa         <  2:4053331/150‑1 (MQ=255)
                                                                                                                                                gCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAGGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa      <  1:34228/150‑1 (MQ=255)
                                                                                                                                                gCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAGGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa      <  1:2716864/150‑1 (MQ=255)
                                                                                                                                                gCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa      <  1:3306933/150‑1 (MQ=255)
                                                                                                                                                 cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac     >  2:4098720/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAAC  >  NC_000913/1708560‑1708857

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