Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,156 G→T 7.9% intergenic (+104/+141) 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,1560GT 7.9% 243.2 / 5.5 74intergenic (+104/+141)mak/araJmanno(fructo)kinase/putative transport protein AraJ
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base T (2/4);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACCTACCCCA  >  NC_000913/411007‑411304
                                                                                                                                                     |                                                                                                                                                    
ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                      <  2:859112/150‑1 (MQ=255)
ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                      <  1:1242670/150‑1 (MQ=255)
  ccAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                                    >  1:167490/1‑150 (MQ=255)
     gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  2:1802537/150‑1 (MQ=38)
     gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  2:1795219/150‑1 (MQ=38)
     gcggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  2:1008789/150‑1 (MQ=38)
        gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                              <  1:1079152/150‑1 (MQ=255)
        gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                              <  2:659356/150‑1 (MQ=255)
        gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                              <  2:221164/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                                                                                                                   <  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: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: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                                                 >  2:1884641/1‑150 (MQ=255)
                                                                                                     ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                                 >  1:1994334/1‑150 (MQ=255)
                                                                                                     ggCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                                 >  2:258503/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                                          >  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                                                               <  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)
                                                                                                                                                 cgccgcATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTAAAACCTAcc     >  2:216503/1‑150 (MQ=255)
                                                                                                                                                  gccgcATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACCTAccc    >  1:121743/1‑150 (MQ=255)
                                                                                                                                                    cgcATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACCTACCCCa  >  1:1590994/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTAATGAGCACCGAGACTATAACCTACCCCA  >  NC_000913/411007‑411304

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: