Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,142 A→G 10.4% intergenic (+90/+155) 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,1420AG10.4% 240.2 / 15.1 86intergenic (+90/+155)mak/araJmanno(fructo)kinase/putative transport protein AraJ
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base G (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 = 1.00e+00

CGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGA  >  NC_000913/410994‑411287
                                                                                                                                                    |                                                                                                                                                 
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:1008789/150‑1 (MQ=38)
                  gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                <  2:1802537/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                                                                                                       >  2:346173/1‑150 (MQ=255)
                                           aTGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGAt                                                                                                       >  1:515209/1‑150 (MQ=255)
                                                cacaAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAc                                                                                                  <  1:1198234/150‑1 (MQ=255)
                                                cacaAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAc                                                                                                  >  2:1752815/1‑150 (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:147486/150‑1 (MQ=255)
                                                             aaCGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCt                                                                                     <  2:140324/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:1796211/150‑1 (MQ=255)
                                                                                             gACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGt                                                     >  1:725846/1‑150 (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                                     >  1:684619/1‑143 (MQ=255)
                                                                                                                    ccTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                     <  2:684619/143‑1 (MQ=255)
                                                                                                                     cTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccg                             >  2:173407/1‑150 (MQ=255)
                                                                                                                       aCAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCa                           >  1:206851/1‑150 (MQ=255)
                                                                                                                         aaaaGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATc                         >  1:1885546/1‑150 (MQ=255)
                                                                                                                         aaaaGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATc                         >  2:706514/1‑150 (MQ=255)
                                                                                                                             ggTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGc                     <  2:1086319/150‑1 (MQ=255)
                                                                                                                              gTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCa                    <  1:1469502/150‑1 (MQ=255)
                                                                                                                               tGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCaa                   >  1:878860/1‑150 (MQ=255)
                                                                                                                               tGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCaa                   >  2:1570673/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:1338680/1‑150 (MQ=21)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAc     >  2:1978695/1‑150 (MQ=21)
                                                                                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                              <  2:1456734/109‑1 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                              >  1:1456734/1‑109 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAcc    <  2:1185479/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCga  >  1:1321427/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGA  >  NC_000913/410994‑411287

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: