Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,143 A→C 10.5% intergenic (+91/+154) mak → / ← araJ manno(fructo)kinase/putative transport protein AraJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913411,1430AC10.5% 242.0 / 20.6 86intergenic (+91/+154)mak/araJmanno(fructo)kinase/putative transport protein AraJ
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base C (4/5);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

CGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTA  >  NC_000913/410994‑411292
                                                                                                                                                     |                                                                                                                                                     
cGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATaa                                                                                                                                                       >  1:1515152/1‑150 (MQ=255)
    gCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCg                                                                                                                                                   <  2:426498/150‑1 (MQ=255)
     cACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTCGGCCGGATAAGGCGt                                                                                                                                                  >  2:696855/1‑150 (MQ=255)
     cACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGt                                                                                                                                                  >  1:1431899/1‑150 (MQ=255)
      aCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGtt                                                                                                                                                 <  2:1834025/150‑1 (MQ=255)
       cGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGtgc                                                                                                                                                >  2:1145864/1‑148 (MQ=255)
       cGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTc                                                                                                                                                <  1:880459/150‑1 (MQ=255)
        ggTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCa                                                                                                                                               <  2:1060313/150‑1 (MQ=255)
          tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                             >  1:1818646/1‑150 (MQ=255)
          tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                             <  2:1703540/150‑1 (MQ=255)
           gATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAGAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                            <  2:1438373/150‑1 (MQ=255)
            aTTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAAcgcc                                                                                                                                           >  1:1802537/1‑150 (MQ=255)
            aTTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAAcgcc                                                                                                                                           >  1:1795219/1‑150 (MQ=255)
             ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                          <  1:1242670/150‑1 (MQ=255)
             ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                          <  2:859112/150‑1 (MQ=255)
               ccAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                        >  1:167490/1‑150 (MQ=255)
                 aGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                      <  2:1293322/150‑6 (MQ=38)
                  gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  2:1795219/150‑1 (MQ=38)
                  gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  2:1802537/150‑1 (MQ=38)
                  gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  2:1008789/150‑1 (MQ=38)
                     gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                  <  2:221164/150‑1 (MQ=255)
                     gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                  <  2:659356/150‑1 (MQ=255)
                     gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                  <  1:1079152/150‑1 (MQ=255)
                       gTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                <  2:589124/150‑3 (MQ=37)
                           gcggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAg                                                                                                                            <  2:1279245/150‑1 (MQ=255)
                                gcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTc                                                                                                                       <  2:331438/150‑1 (MQ=255)
                                     gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATg                                                                                                                  >  2:943584/1‑150 (MQ=255)
                                        gTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCt                                                                                                               <  1:766292/150‑1 (MQ=255)
                                         ttATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTg                                                                                                              <  2:637911/150‑1 (MQ=255)
                                           aTGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGAGTGCTCCGATGCCTGAt                                                                                                            <  1:1092595/150‑1 (MQ=255)
                                           aTGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGAt                                                                                                            >  1:515209/1‑150 (MQ=255)
                                           aTGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGAt                                                                                                            >  2:346173/1‑150 (MQ=255)
                                                cacaAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAc                                                                                                       >  2:1752815/1‑150 (MQ=255)
                                                cacaAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAc                                                                                                       <  1:1198234/150‑1 (MQ=255)
                                                   aaGAGTAAAAAACTTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCt                                                                                                    <  1:1167620/150‑1 (MQ=255)
                                                   aaGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCt                                                                                                    <  1:1175180/150‑1 (MQ=255)
                                                     gagTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCAGATGCCTGATGCGACGCTTg                                                                                                  >  1:830695/1‑150 (MQ=255)
                                                            aaaCGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTc                                                                                           <  2:340333/150‑1 (MQ=255)
                                                             aaCGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCt                                                                                          <  2:140324/150‑1 (MQ=255)
                                                             aaCGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCt                                                                                          <  2:147486/150‑1 (MQ=255)
                                                                 tAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                      <  2:736109/150‑1 (MQ=255)
                                                                  aGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                     <  2:1182779/150‑1 (MQ=255)
                                                                         cgtcgcATCAGGCATGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGc                                                                              >  1:169803/4‑150 (MQ=255)
                                                                              catcatGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAAt                                                                         >  1:1086319/1‑150 (MQ=255)
                                                                                 catGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAAtgtg                                                                      >  1:1185479/1‑150 (MQ=255)
                                                                                 catGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTACACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAAtgtg                                                                      <  2:1990313/150‑1 (MQ=255)
                                                                                    gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCa                                                                   >  2:268658/1‑150 (MQ=255)
                                                                                          tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACcg                                                             >  1:1945705/1‑150 (MQ=255)
                                                                                             gACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGt                                                          >  1:725846/1‑150 (MQ=255)
                                                                                             gACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGt                                                          <  1:1796211/150‑1 (MQ=255)
                                                                                              aCGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTa                                                         <  2:1431899/150‑1 (MQ=255)
                                                                                                    gCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                                   <  2:515209/150‑1 (MQ=255)
                                                                                                    gCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                                   <  2:169803/150‑1 (MQ=255)
                                                                                                      cgcgTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGAt                                                 <  1:1752815/150‑1 (MQ=255)
                                                                                                       gcgTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATa                                                <  1:270213/150‑1 (MQ=255)
                                                                                                           cTTATCAGGCCTACAAAAGGTGCCAGAACAGTAGGCCGCATAAGTCGATCACTCCGCATACGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGc                                            >  2:774449/1‑150 (MQ=37)
                                                                                                             tATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                          >  1:1507589/1‑150 (MQ=255)
                                                                                                             tATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGAAAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                          >  1:363618/1‑150 (MQ=255)
                                                                                                              aTCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGtt                                         >  2:1185003/1‑150 (MQ=255)
                                                                                                               tCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTc                                        >  1:883568/1‑150 (MQ=255)
                                                                                                                cAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCa                                       <  2:1818646/150‑1 (MQ=255)
                                                                                                                  ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                     >  2:1884641/1‑150 (MQ=255)
                                                                                                                  ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                     >  2:258503/1‑150 (MQ=255)
                                                                                                                  ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                     >  1:1994334/1‑150 (MQ=255)
                                                                                                                    ccTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                          <  2:684619/143‑1 (MQ=255)
                                                                                                                    ccTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                          >  1:684619/1‑143 (MQ=255)
                                                                                                                     cTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccg                                  >  2:173407/1‑150 (MQ=255)
                                                                                                                       aCAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCa                                >  1:206851/1‑150 (MQ=255)
                                                                                                                         aaaaGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATc                              >  2:706514/1‑150 (MQ=255)
                                                                                                                         aaaaGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATc                              >  1:1885546/1‑150 (MQ=255)
                                                                                                                             ggTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGc                          <  2:1086319/150‑1 (MQ=255)
                                                                                                                              gTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCa                         <  1:1469502/150‑1 (MQ=255)
                                                                                                                               tGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCaa                        >  2:1570673/1‑150 (MQ=255)
                                                                                                                               tGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCaa                        >  1:878860/1‑150 (MQ=255)
                                                                                                                                gCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCCGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAAt                       >  1:1439255/1‑150 (MQ=255)
                                                                                                                                 ccAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATa                      >  2:1938188/1‑150 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGc            <  1:1884641/150‑1 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCAGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGc            <  2:264287/150‑1 (MQ=255)
                                                                                                                                           aGGACGTATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGc            <  2:363618/150‑1 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCa           >  1:1167332/1‑150 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAAAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCa           >  1:1006320/1‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAc          >  2:1978695/1‑150 (MQ=21)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAc          >  2:1338680/1‑150 (MQ=21)
                                                                                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                                   >  1:1456734/1‑109 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                                   <  2:1456734/109‑1 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAcc         <  2:1185479/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCga       >  1:1321427/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACta  >  1:681583/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAAATAGTAATGAGCACCGAGACta  <  1:1693955/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTA  >  NC_000913/410994‑411292

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