Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 157,585 C→T D50N (GAT→AAT)  folK ← 2‑amino‑4‑hydroxy‑6‑hydroxymethyldihyropteridine pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913157,5850CT97.1% 1405.5 / 19.1 419D50N (GAT→AAT) folK2‑amino‑4‑hydroxy‑6‑hydroxymethyldihyropteridine pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base T (181/226);  total (186/233)
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.52e-01

Only 100 of 421 total aligned reads displayed.

ATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCAT  >  NC_000913/157437‑157732
                                                                                                                                                    |                                                                                                                                                   
atgTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGAtttt                                                                                                                                                   <  3:1132128/151‑1 (MQ=255)
  gTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTTCGTTTAAGTAATCGCGTTGATTTTGc                                                                                                                                                 >  6:1577594/1‑151 (MQ=255)
    cGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCg                                                                                                                                                <  2:1278599/150‑1 (MQ=255)
    cGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCg                                                                                                                                                <  2:1574974/150‑1 (MQ=255)
     gagaTCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAACACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCTCGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGc                                                                                                                                              >  3:1513727/1‑151 (MQ=255)
      cgctCCCGCGGGCGTGGTCCCCAGCGTTCAGCTTTGCGGCCGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGGAGCTCTTCAGGTGCAAGAGCGGTGTCCAGCGCCACGGCTGCGGTTAAGTAATCGGGTTGATTTTGCGGcc                                                                                                                                             <  2:1191080/148‑1 (MQ=255)
      agaTCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGcc                                                                                                                                             <  1:399790/151‑1 (MQ=255)
          ccAGCGTGCGTGTTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAg                                                                                                                                         <  4:376643/151‑1 (MQ=255)
            aGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCg                                                                                                                                       >  7:881700/1‑151 (MQ=255)
            aGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCg                                                                                                                                       <  5:211741/151‑1 (MQ=255)
            aGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCg                                                                                                                                       >  8:1211310/1‑151 (MQ=255)
               gtgcgtgGTCGCCAGCGTGCAGCTTTGCGGACGCGAGCTTGCTGCAATTCAATACGGTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGCGTTGATTTTGCGGCGCCAGCGGTg                                                                                                                                    <  1:1432930/151‑1 (MQ=255)
               gtgcgtgGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTg                                                                                                                                    >  3:736481/1‑151 (MQ=255)
                 gcgtgGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTgg                                                                                                                                   <  8:717139/150‑1 (MQ=255)
                     ggTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGc                                                                                                                              <  4:1488214/151‑1 (MQ=255)
                       tCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCgg                                                                                                                            >  8:403580/1‑151 (MQ=255)
                          ccAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTa                                                                                                                          >  2:946610/1‑150 (MQ=255)
                            aGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTaaaa                                                                                                                       <  5:1419876/151‑1 (MQ=255)
                             gCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAAc                                                                                                                      <  6:528160/151‑1 (MQ=255)
                               gTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGGTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACg                                                                                                                     >  5:1370275/1‑150 (MQ=255)
                                 tCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCCGCGGTGGGGTGCGGTAAAACga                                                                                                                    >  1:1627194/1‑149 (MQ=255)
                                 tCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACgaa                                                                                                                   <  7:1184544/150‑1 (MQ=255)
                                    gCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGaaa                                                                                                               <  4:616695/151‑1 (MQ=255)
                                       ttGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACgg                                                                                                            >  2:384891/1‑151 (MQ=255)
                                       ttGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACg                                                                                                             <  4:165589/150‑1 (MQ=255)
                                        tGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGt                                                                                                           >  3:929446/1‑151 (MQ=255)
                                         gCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGt                                                                                                           <  2:1377003/150‑1 (MQ=255)
                                         gCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTa                                                                                                          <  7:132857/151‑1 (MQ=255)
                                             aCGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGa                                                                                                       >  7:1486553/1‑150 (MQ=255)
                                              cgcgACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAAt                                                                                                     <  2:1138891/151‑1 (MQ=255)
                                                 gACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAAGGGTAAGAAtgtg                                                                                                  >  6:855036/1‑151 (MQ=255)
                                                  aCCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTgg                                                                                                 >  6:1226302/1‑151 (MQ=255)
                                                  aCCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTgg                                                                                                 <  3:246875/151‑1 (MQ=255)
                                                      tgctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTGTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCt                                                                                               <  6:1036489/149‑1 (MQ=255)
                                                       gctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCttt                                                                                             <  3:56224/150‑1 (MQ=255)
                                                       gctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                            >  6:177954/1‑151 (MQ=255)
                                                        ctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                            >  3:508846/1‑150 (MQ=255)
                                                        ctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                            >  1:235650/1‑150 (MQ=255)
                                                         tgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                            >  7:285047/1‑149 (MQ=255)
                                                            aaTTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGa                                                                                       <  7:710748/151‑1 (MQ=255)
                                                             aTTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGat                                                                                      <  8:313765/151‑1 (MQ=255)
                                                               tCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGatat                                                                                    >  6:1626452/1‑151 (MQ=255)
                                                               tCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGatat                                                                                    >  6:1162160/1‑151 (MQ=255)
                                                                cAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGatat                                                                                    <  2:339357/150‑1 (MQ=255)
                                                                cAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATc                                                                                   <  6:45898/151‑1 (MQ=255)
                                                                   tACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGcc                                                                                >  7:453083/1‑151 (MQ=255)
                                                                    aCGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCt                                                                               <  5:733513/151‑1 (MQ=255)
                                                                       cTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAAt                                                                            >  8:1511414/1‑151 (MQ=255)
                                                                           gtgtgtTTGAGTAGCTCTTCATGTTGAAGTGCTGTTTCCATCGCCACGGCGGGGTTTAAGTAATCGGGTTGATTTTGCGGCCCCATCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCtt                                                                        <  1:1089400/151‑1 (MQ=255)
                                                                           gtgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCt                                                                         <  2:627995/150‑1 (MQ=255)
                                                                             gtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTc                                                                      <  2:291021/151‑1 (MQ=255)
                                                                              tgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTc                                                                      >  4:24377/1‑150 (MQ=255)
                                                                               gATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGGTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCa                                                                     <  4:1428345/150‑1 (MQ=255)
                                                                                  tGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgc                                                                 >  8:144775/1‑151 (MQ=255)
                                                                                   gAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgca                                                                <  4:527881/151‑1 (MQ=255)
                                                                                       aGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCAt                                                            <  1:1274420/151‑1 (MQ=255)
                                                                                           ctTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGACTGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTg                                                          <  5:432335/149‑1 (MQ=255)
                                                                                           ctTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGAc                                                        >  3:1059682/1‑151 (MQ=255)
                                                                                             tCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGa                                                         <  2:771357/148‑1 (MQ=255)
                                                                                             tCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGa                                                         <  3:137566/148‑1 (MQ=255)
                                                                                             tCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGa                                                         <  3:145159/148‑1 (MQ=255)
                                                                                               aGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGc                                                    >  8:585941/1‑151 (MQ=255)
                                                                                                ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCt                                                   <  3:837924/151‑1 (MQ=255)
                                                                                                 gTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTc                                                  >  5:1166980/1‑151 (MQ=255)
                                                                                                  tGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGCTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTcc                                                 >  7:1587454/1‑151 (MQ=255)
                                                                                                   gCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCa                                                <  1:331383/151‑1 (MQ=255)
                                                                                                    cAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAg                                               <  5:370687/151‑1 (MQ=255)
                                                                                                     aaGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGc                                              >  1:644505/1‑151 (MQ=255)
                                                                                                            gTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAgg                                       <  1:937271/151‑1 (MQ=255)
                                                                                                             tttCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGc                                      >  3:614211/1‑151 (MQ=255)
                                                                                                              ttCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGcc                                     >  1:780680/1‑151 (MQ=255)
                                                                                                               tCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGcc                                     <  7:38412/150‑1 (MQ=255)
                                                                                                                   gcgcCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGAtt                                <  7:1016438/151‑1 (MQ=255)
                                                                                                                    cgcCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATtg                               >  1:1147601/1‑151 (MQ=255)
                                                                                                                     gcCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATtgc                              <  5:642368/151‑1 (MQ=255)
                                                                                                                        aCGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATtgctgc                           >  3:198659/1‑151 (MQ=255)
                                                                                                                         cGGGTGCTTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGcc                          <  5:967748/151‑1 (MQ=255)
                                                                                                                          ggCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCt                         <  3:302553/151‑1 (MQ=255)
                                                                                                                           gCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTa                        <  4:1562216/151‑1 (MQ=255)
                                                                                                                           gCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTACAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTa                        <  5:1513940/151‑1 (MQ=255)
                                                                                                                              gCGTTTAAGTAATCGGGTTGATTTTGCGGCCCNAGCGGTGGGGTGCGGNAAAACGAAGAAACGNTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCANCGGAGAGGCNAGATTGCTGCCTATg                      <  8:1144981/150‑1 (MQ=255)
                                                                                                                              gCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATgg                     >  5:889358/1‑151 (MQ=255)
                                                                                                                               cGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGc                    <  1:1613949/151‑1 (MQ=255)
                                                                                                                                  ttAAGTAATCGTGTTGATTTTGCGGCCGCTGCGGTGTGTTGCGGTATAACGATGAAACGGTTAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTGATTGCAGCAGTGTCCTGATCCAGTGGAGAGGCCAGATTTCTGCCTATGGCaa                  <  1:502046/150‑1 (MQ=255)
                                                                                                                                  ttAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAat                 >  5:561577/1‑151 (MQ=255)
                                                                                                                                  ttAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAat                 >  7:179126/1‑151 (MQ=255)
                                                                                                                                  ttAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAat                 >  4:1279215/1‑151 (MQ=255)
                                                                                                                                     aGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAAtata              <  8:676535/151‑1 (MQ=255)
                                                                                                                                      gTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATAc             <  5:1005619/151‑1 (MQ=255)
                                                                                                                                       tAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg            >  2:207314/1‑151 (MQ=255)
                                                                                                                                       tAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg            >  8:1681414/1‑151 (MQ=255)
                                                                                                                                       gcctCGGGTTGATTTTGCGGCCCCAGCGGGTGGGGGCGGTAAAAGGAAGAAACGGTACGAATGGGGCGTTCAGGGATTTCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg            <  6:31659/148‑1 (MQ=255)
                                                                                                                                        aaTCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg            <  6:94604/150‑1 (MQ=255)
                                                                                                                                         aTCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGcc          <  4:813580/151‑1 (MQ=255)
                                                                                                                                           cGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACTCCAc        >  1:610619/1‑151 (MQ=255)
                                                                                                                                           cGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCa         <  1:1387933/150‑1 (MQ=255)
                                                                                                                                           cGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCa         <  4:1096478/150‑1 (MQ=255)
                                                                                                                                            gggTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACt       <  6:127036/151‑1 (MQ=255)
                                                                                                                                             gggtGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATTGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTg      <  4:870621/148‑1 (MQ=255)
                                                                                                                                                 gATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCa   <  6:676161/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCAT  >  NC_000913/157437‑157732

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