Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 864,307 G→C 5.7% intergenic (+3/+73) fsaA → / ← moeB fructose‑6‑phosphate aldolase 1/molybdopterin synthase sulfurylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913864,3070GC 5.7% 439.0 / 13.6 141intergenic (+3/+73)fsaA/moeBfructose‑6‑phosphate aldolase 1/molybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base G (69/64);  new base C (5/3);  total (74/67)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 145 total aligned reads displayed.

CGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTACA  >  NC_000913/864159‑864456
                                                                                                                                                    |                                                                                                                                                     
cgcTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCgg                                                                                                                                                      <  2:1790866/150‑1 (MQ=255)
 gcTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGa                                                                                                                                                     >  1:1335356/1‑150 (MQ=255)
   tGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACg                                                                                                                                                   <  1:795335/150‑1 (MQ=255)
   tGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGAGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACg                                                                                                                                                   >  1:1282253/1‑150 (MQ=255)
        tGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTg                                                                                                                                              <  2:1297989/150‑1 (MQ=255)
         gCTTACTGGCAGGAGGTGAAACATTTACTCTGCCCCTCGAGGTGTCACAACAGCGGAAAAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCCGGTTAACTCGACGTTCTGc                                                                                                                                             <  1:1035103/150‑1 (MQ=255)
          cTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCa                                                                                                                                            <  1:938051/150‑1 (MQ=255)
            tACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATc                                                                                                                                          >  2:1008960/1‑150 (MQ=255)
            tACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATc                                                                                                                                          >  2:1836906/1‑150 (MQ=255)
               tGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTc                                                                                                                                       <  2:1017338/150‑1 (MQ=255)
                        gtgAATCAATTACTCAGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGc                                                                                                                              <  2:1120361/150‑1 (MQ=255)
                          gAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCtg                                                                                                                            >  2:1951822/1‑150 (MQ=255)
                           aaTCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCtga                                                                                                                           <  1:835209/150‑1 (MQ=255)
                           aaTCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCtga                                                                                                                           >  1:1990788/1‑150 (MQ=255)
                               aaTTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAAATGCTGATGAc                                                                                                                       >  1:1048166/1‑150 (MQ=255)
                               aaTTACTCTGCCACTGGAAGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGAATGGCAGGGAGCGTTTGGCAGAACGTCGATTAAACCGGACGTTCTGCATCCTCATAAATTGCTGATGAc                                                                                                                       >  1:1739682/1‑150 (MQ=255)
                                  tACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTg                                                                                                                    <  2:894646/150‑1 (MQ=255)
                                    ctctGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGc                                                                                                                  >  2:1682232/1‑150 (MQ=255)
                                    ctctGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGc                                                                                                                  >  2:137140/1‑150 (MQ=255)
                                    ctctGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGc                                                                                                                  >  2:1171695/1‑150 (MQ=255)
                                      ctGCCACGGGATGGGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGCGCAGGACTGGCAGGGAGCGATTGGCAGAAAGTCGATTTACCCGGACGGTCTGCATCCTCAGTGAATGCTGATGCCTTCGCgg                                                                                                                >  2:1817005/1‑150 (MQ=255)
                                           aCTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGc                                                                                                           <  2:590367/150‑1 (MQ=255)
                                           aCTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTCCATCCTCATAAATTGCTGATGACGTGGCGGAGTGc                                                                                                           <  2:1823129/150‑1 (MQ=255)
                                            cttgATGTGGCACAACTGTTGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGTACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGcc                                                                                                          <  1:1999690/147‑1 (MQ=255)
                                            cTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGcc                                                                                                          <  2:1311646/150‑1 (MQ=255)
                                             tGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCcg                                                                                                         <  1:1625650/150‑1 (MQ=255)
                                                 tgtgGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAAAGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGAAGTGGCGGAGTGCCGCGTc                                                                                                     >  1:61788/1‑150 (MQ=255)
                                                   tgGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGCCGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCta                                                                                                   >  1:478201/1‑149 (MQ=255)
                                                   tgGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCtt                                                                                                   >  1:1118679/1‑150 (MQ=255)
                                                   tgGCACAACAGATGATTAGCTATCCCGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGCGAGCGTTAGGCAGAACGTCGATTTAACCGCACGTTCGGAATCCTCATAAATTGCTGATGACTTGCAGCAGTGCCGCGTCtt                                                                                                   >  2:858492/1‑150 (MQ=255)
                                                      cacaACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGCCGTGGCGGAGTGCCGCGTCTTATc                                                                                                >  2:1047405/1‑150 (MQ=255)
                                                      cacaACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATc                                                                                                >  2:11582/1‑150 (MQ=255)
                                                      cacaACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATc                                                                                                >  2:2028376/1‑150 (MQ=255)
                                                      cacaACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATc                                                                                                >  2:1534624/1‑150 (MQ=255)
                                                         aacaGATGAGTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAgg                                                                                             <  2:113887/150‑1 (MQ=255)
                                                           caGATGCTTAGCTATCCGGCGGTTGATGCCGCTGTGGCTAAGTGTGAGCAGGACTGTCAGGGAGCGTTTGGCAGATCGTCTTTTTAACCGGACGTTCTGCATCCTCATATATTGCTGATGACGGGGCGGAGGGCCGCGTATTATCAGGcc                                                                                           <  1:313059/150‑1 (MQ=255)
                                                            aGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCt                                                                                          <  1:1752212/150‑1 (MQ=255)
                                                            aGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCt                                                                                          >  1:870571/1‑150 (MQ=255)
                                                              atgatTAGCTATCCGGAGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTgg                                                                                        <  2:1282252/150‑1 (MQ=255)
                                                               tgatTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGa                                                                                       >  2:935987/1‑150 (MQ=255)
                                                                 atTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAgg                                                                                     <  2:191205/150‑1 (MQ=255)
                                                                     gCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGc                                                                                 <  2:732957/150‑1 (MQ=255)
                                                                     gCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGGTTGGCAGAACGGCGATTTAACCGGACGTGCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGc                                                                                 <  1:1867193/150‑1 (MQ=255)
                                                                     gCTAGCCTGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGc                                                                                 <  1:979030/150‑1 (MQ=255)
                                                                       tATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCaa                                                                               >  2:602069/1‑150 (MQ=255)
                                                                           cggcggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTAc                                                                           >  1:1666373/1‑150 (MQ=255)
                                                                            ggcggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACt                                                                          >  1:1470980/1‑150 (MQ=255)
                                                                            ggcggTTGATGCCGCTGAGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTAGATTAAACCGGACCTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACt                                                                          >  1:510827/1‑150 (MQ=255)
                                                                               ggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGcc                                                                       >  1:499637/1‑150 (MQ=255)
                                                                               ggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGcc                                                                       >  1:1445209/1‑150 (MQ=255)
                                                                               ggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGcc                                                                       <  1:1680817/150‑1 (MQ=255)
                                                                               ggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGcc                                                                       <  1:201957/150‑1 (MQ=255)
                                                                                  tGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGGTGGGCAGAACGGCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCcac                                                                    <  2:276191/150‑1 (MQ=255)
                                                                                    aTGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCCGGCCTGGAGGTGGCAATTACTGCCcacac                                                                  >  1:1297433/1‑150 (MQ=255)
                                                                                      gCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCacacac                                                                >  1:1526747/1‑150 (MQ=255)
                                                                                      gCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACATCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCacacac                                                                >  1:1919065/1‑150 (MQ=255)
                                                                                       ccGCTGGGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGAGTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACAcc                                                               <  1:223590/150‑1 (MQ=255)
                                                                                                 gAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCg                                                     <  2:1715088/150‑1 (MQ=255)
                                                                                                  aaGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCgg                                                    <  1:411714/150‑1 (MQ=255)
                                                                                                   aGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGa                                                   <  1:1212474/150‑1 (MQ=255)
                                                                                                      ttGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTa                                                >  1:626103/1‑150 (MQ=255)
                                                                                                      ttGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTAAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTa                                                <  2:261691/150‑1 (MQ=255)
                                                                                                        gAGCAGGACTGGAAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGGGGCAATTACTGCCCACACACCTCACACCCCGGATTACg                                              <  2:1779947/150‑1 (MQ=255)
                                                                                                         aGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGc                                             <  2:1384242/150‑1 (MQ=255)
                                                                                                          gCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCa                                            <  2:917150/150‑1 (MQ=255)
                                                                                                          gCAGGACTGGCAGGGAGCGTTTGGCAGAACGGCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCa                                            <  2:1541615/150‑1 (MQ=255)
                                                                                                             ggACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCCCACACCTCACACCCCGGATTACGCATCa                                         >  2:1362682/1‑150 (MQ=255)
                                                                                                             ggACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCa                                         >  2:423626/1‑150 (MQ=255)
                                                                                                                cTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCGGGATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGtt                                      >  1:1865739/1‑150 (MQ=255)
                                                                                                                 tGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGttt                                     <  1:1925088/150‑1 (MQ=255)
                                                                                                                  ggCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTc                                    <  2:1074659/150‑1 (MQ=255)
                                                                                                                  ggCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTc                                    >  1:726587/1‑150 (MQ=255)
                                                                                                                   gCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCa                                   >  2:1855347/1‑150 (MQ=255)
                                                                                                                   gCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCa                                   <  2:2011400/150‑1 (MQ=255)
                                                                                                                    cAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCAt                                  <  2:253120/150‑1 (MQ=255)
                                                                                                                    cAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGAGTACGCATCAGTTTCAt                                  >  1:1032059/1‑150 (MQ=255)
                                                                                                                      gggAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCAttt                                >  2:499113/1‑150 (MQ=255)
                                                                                                                       ggAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTc                               >  2:1708988/1‑150 (MQ=255)
                                                                                                                         aGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcgc                             <  1:11582/150‑1 (MQ=255)
                                                                                                                         aGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcgc                             >  2:292517/1‑150 (MQ=255)
                                                                                                                          gCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACTCCTCACACCCCGGATTACGCATCAGTTTCATTTcgcg                            <  1:1682233/150‑1 (MQ=255)
                                                                                                                           cGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCATCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGACAACACACATCACAACCCGGATTTCGCATCAGTTTGATTTAGCGa                           >  1:1068915/1‑150 (MQ=255)
                                                                                                                            gTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGAGTCTTATCAGGCCTGGAGGTTGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGaa                          <  2:1431071/150‑1 (MQ=255)
                                                                                                                             ggatgCAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGaaa                         >  2:505810/5‑150 (MQ=255)
                                                                                                                                ggCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTg                      >  2:615762/1‑150 (MQ=255)
                                                                                                                                ggCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTg                      >  2:475189/1‑150 (MQ=255)
                                                                                                                                 gCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAACCTGa                     >  1:1523813/1‑150 (MQ=255)
                                                                                                                                 gCAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCCCACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGa                     >  1:584185/1‑150 (MQ=255)
                                                                                                                                 gCAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTgc                     >  1:2049968/1‑149 (MQ=255)
                                                                                                                                 ccAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTAATGCCCACACACCTCACACCCCGGATTTCGCATCAGTTTCATTTCGCGAAACTGa                     >  1:16665/2‑150 (MQ=255)
                                                                                                                                   aGAACGTAGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTCGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCTAAACTGACa                   <  2:1941711/150‑1 (MQ=255)
                                                                                                                                    gAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAg                  <  1:1603278/150‑1 (MQ=255)
                                                                                                                                           gATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatcg           <  1:1225807/150‑1 (MQ=255)
                                                                                                                                               tAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatcgcatc       >  2:2022221/1‑150 (MQ=255)
                                                                                                                                                aaCCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatcgcatcg      <  1:1708989/150‑1 (MQ=255)
                                                                                                                                                aaCCGGACGTGCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGACCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatcgcatcg      <  1:302606/150‑1 (MQ=255)
                                                                                                                                                  ccGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTa    <  1:912858/150‑1 (MQ=255)
                                                                                                                                                    ggACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTACa  <  1:410215/150‑1 (MQ=255)
                                                                                                                                                    ggACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTACa  <  1:103637/150‑1 (MQ=255)
                                                                                                                                                    ggACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTACa  <  1:602069/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATCGCATCGTACA  >  NC_000913/864159‑864456

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