Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,077 T→G 5.7% intergenic (+25/+220) 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,0770TG 5.7% 222.1 / 4.5 70intergenic (+25/+220)mak/araJmanno(fructo)kinase/putative transport protein AraJ
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base G (2/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCC  >  NC_000913/410931‑411226
                                                                                                                                                  |                                                                                                                                                     
cGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGcc                                                                                                                                                    <  1:1197280/150‑1 (MQ=255)
           ttGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGc                                                                                                                                         <  1:928752/150‑1 (MQ=255)
           ttGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGc                                                                                                                                         <  2:15088/150‑1 (MQ=255)
            tGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCt                                                                                                                                        >  1:1248344/1‑150 (MQ=255)
             gATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCtt                                                                                                                                       >  1:443882/1‑150 (MQ=255)
               ttAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGc                                                                                                                                     >  2:1069711/1‑150 (MQ=255)
                tAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCAGGCATGATGCGCCAATTGcc                                                                                                                                    <  1:409389/150‑1 (MQ=255)
                 aaaCAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCcg                                                                                                                                   >  1:499982/1‑150 (MQ=255)
                 aaaCAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCcg                                                                                                                                   >  2:1104322/1‑150 (MQ=255)
                   aCAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCcgcg                                                                                                                                 >  2:345634/1‑150 (MQ=255)
                   aCAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCcgcg                                                                                                                                 <  2:209280/150‑1 (MQ=255)
                    cAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                >  2:1151576/1‑150 (MQ=255)
                    cAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                <  1:171272/150‑1 (MQ=255)
                         tGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCCCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                           >  1:1065391/1‑150 (MQ=255)
                         tGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                           >  2:61165/1‑150 (MQ=255)
                          gTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                          <  2:1248344/150‑1 (MQ=255)
                           tCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                         <  1:388375/150‑1 (MQ=255)
                            cTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                        >  1:719374/1‑150 (MQ=255)
                            cTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                        <  1:1103762/150‑1 (MQ=255)
                                 gcggcgAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtt                                                                                                                   >  1:1214508/1‑149 (MQ=255)
                                     cgAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaaa                                                                                                               >  2:1134993/1‑150 (MQ=255)
                                     cgAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaaa                                                                                                               >  1:901481/1‑150 (MQ=255)
                                     cgAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCAGGc                                                                                                                                                     >  1:170657/1‑112 (MQ=255)
                                     cgAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCAGGc                                                                                                                                                     <  2:170657/112‑1 (MQ=255)
                                         tgtgAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGt                                                                                                           >  1:1420772/1‑150 (MQ=255)
                                         tgtgAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGt                                                                                                           >  1:406916/1‑150 (MQ=255)
                                          gtgAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTg                                                                                                          <  2:1167080/150‑1 (MQ=255)
                                           tgAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGc                                                                                                         >  1:515104/1‑150 (MQ=255)
                                           tgAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGc                                                                                                         >  2:246925/1‑150 (MQ=255)
                                            gAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGcc                                                                                                        >  2:688682/1‑150 (MQ=255)
                                            gAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGcc                                                                                                        <  2:557064/150‑1 (MQ=255)
                                            gAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGcc                                                                                                        >  1:1088669/1‑150 (MQ=255)
                                            gAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGcc                                                                                                        >  1:1411270/1‑150 (MQ=255)
                                              aaCGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAg                                                                                                      <  1:1104322/150‑1 (MQ=255)
                                                   cGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGATCCg                                                                                                 <  2:214771/150‑1 (MQ=255)
                                                         gTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGcc                                                                                           <  1:796698/150‑1 (MQ=255)
                                                           aaGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCgg                                                                                         <  2:1272595/150‑1 (MQ=255)
                                                              gCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATa                                                                                      <  1:684769/150‑1 (MQ=255)
                                                               cGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATaa                                                                                     >  2:724907/1‑150 (MQ=255)
                                                                 aaGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAgg                                                                                   <  1:1065844/150‑1 (MQ=255)
                                                                 aaGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAgg                                                                                   <  1:1061015/150‑1 (MQ=255)
                                                                 aaGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGGGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAgg                                                                                   <  1:1218704/150‑1 (MQ=255)
                                                                     aCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGtt                                                                               >  2:765832/1‑150 (MQ=255)
                                                                           aTTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgcc                                                                         <  1:246925/150‑1 (MQ=255)
                                                                            ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                        <  1:1069711/150‑1 (MQ=255)
                                                                                   ggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                 <  2:1267543/150‑1 (MQ=255)
                                                                                    gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCCACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                <  1:345634/150‑1 (MQ=255)
                                                                                      gTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                              <  2:901481/150‑3 (MQ=37)
                                                                                        aCGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATa                                                            >  2:1349276/1‑150 (MQ=255)
                                                                                         cgcgGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTAGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                                                           >  2:144126/1‑150 (MQ=255)
                                                                                               gcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTc                                                     <  1:1134993/150‑1 (MQ=255)
                                                                                                  gCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGa                                                  >  2:835119/1‑150 (MQ=255)
                                                                                                  gCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGa                                                  >  2:1162900/1‑150 (MQ=255)
                                                                                                  gCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGa                                                  >  2:247693/1‑150 (MQ=255)
                                                                                                    gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATg                                                <  2:499982/150‑1 (MQ=255)
                                                                                                        ttATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTg                                            <  2:515104/150‑1 (MQ=255)
                                                                                                           tGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATg                                         <  1:61165/150‑1 (MQ=255)
                                                                                                                    gagTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTg                                <  1:724907/150‑1 (MQ=255)
                                                                                                                     agTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGc                               >  1:496677/1‑150 (MQ=255)
                                                                                                                                tAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATc                    <  2:1420772/150‑1 (MQ=255)
                                                                                                                                tAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATc                    <  2:1014249/150‑1 (MQ=255)
                                                                                                                                 tggCACTTGGGGCAGCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGGGACGCTTGCCGCGTCTTATCa                   <  1:353271/149‑1 (MQ=255)
                                                                                                                                  ggCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTGTCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                  <  1:133993/150‑1 (MQ=255)
                                                                                                                                  ggCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                  >  2:567655/1‑150 (MQ=255)
                                                                                                                                  ggCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                  <  2:985449/150‑1 (MQ=255)
                                                                                                                                  ggCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                  >  2:1039800/1‑150 (MQ=255)
                                                                                                                                  ggCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                  >  1:985449/1‑150 (MQ=255)
                                                                                                                                   gCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                 <  2:1065391/150‑1 (MQ=255)
                                                                                                                                     aaTTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTTCAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc               <  2:1214508/150‑1 (MQ=255)
                                                                                                                                       gcgtcgcATCAGGCATGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTg             >  2:322479/5‑150 (MQ=255)
                                                                                                                                               tcatGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAAtgt     >  2:980847/1‑150 (MQ=255)
                                                                                                                                               tcatGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAAtgt     >  1:399938/1‑150 (MQ=255)
                                                                                                                                                  tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGcc  <  2:1088669/150‑1 (MQ=255)
                                                                                                                                                  tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGcc  >  2:102945/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCC  >  NC_000913/410931‑411226

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