Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,149 T→G 11.1% intergenic (+97/+148) 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,1490TG11.1% 220.1 / 20.6 83intergenic (+97/+148)mak/araJmanno(fructo)kinase/putative transport protein AraJ
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base G (4/5);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

ACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACC  >  NC_000913/411000‑411297
                                                                                                                                                     |                                                                                                                                                    
aCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGtt                                                                                                                                                      <  2:1834025/150‑1 (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:1802537/150‑1 (MQ=38)
            gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                          <  2:1795219/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: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:169803/150‑1 (MQ=255)
                                                                                              gCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCgg                                                        <  2:515209/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                                          >  1:1994334/1‑150 (MQ=255)
                                                                                                            ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                          >  2:258503/1‑150 (MQ=255)
                                                                                                            ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                          >  2:1884641/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:1978695/1‑150 (MQ=21)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCAc               >  2:1338680/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)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACta       >  1:681583/1‑150 (MQ=255)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAAATAGTAATGAGCACCGAGACta       <  1:1693955/150‑1 (MQ=255)
                                                                                                                                                    ttCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAAcc  <  2:1290223/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACC  >  NC_000913/411000‑411297

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: