Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,904,212 G→A intergenic (+838/+535) ygcE → / ← queE putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,904,2120GA95.4% 255.1 / NA 109intergenic (+838/+535)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (52/52);  minor base C (4/1);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCT  >  NC_000913/2904068‑2904355
                                                                                                                                                |                                                                                                                                               
aGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCGATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCCGGGAACACGGGAGAACGTAATAAAAAGCCGTGGAGCTGCGGGGTTTATCCCCGCGGAcgcg                                                                                                                                             >  2:741519/1‑149 (MQ=255)
 ggCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                            <  1:244731/149‑1 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                           >  1:283218/1‑149 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                           >  1:935415/1‑149 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                           >  1:989105/1‑149 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                           <  1:242404/149‑1 (MQ=255)
    tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCggggg                                                                                                                                         >  1:769345/1‑148 (MQ=255)
    tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGa                                                                                                                                         >  1:915442/1‑149 (MQ=255)
    tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGGATCTGCTCGGTTTATCCCCGCTGACGCGGGGa                                                                                                                                         >  1:109104/1‑149 (MQ=255)
     ccGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTAATGCTATTCGGTGTATCCCCGCTGGCGCGGGGAACACGGGCGCAAGGAATACAAAGCCGTGTATCTGCTCTGTTTATCCCCGCTGACGCGGGGaa                                                                                                                                        <  2:1061668/149‑1 (MQ=255)
       ggTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAaca                                                                                                                                      <  1:394121/149‑1 (MQ=255)
           tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATCGACTATTGCTATTCGGATTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTgg                                                                                                                                  <  1:42228/149‑1 (MQ=255)
            aTCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGc                                                                                                                                 <  1:770350/149‑1 (MQ=255)
             tCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGct                                                                                                                                >  2:1017505/1‑149 (MQ=255)
              ccccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctc                                                                                                                               >  2:281958/1‑149 (MQ=255)
               cccGCTGGCGCGTGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctct                                                                                                                              <  2:526614/149‑1 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGTCGCACGGATTACTAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                             >  1:85128/1‑149 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                             <  2:201730/149‑1 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAAGGAATACAAAGCCGTGTATATGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                             >  1:331543/1‑149 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAAGGAATACAAAGCAGTGTATCTGCTAGGTTTATCACCGCAGACGCGGGGACCACTGGctctt                                                                                                                             >  1:498743/1‑148 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAAACTATGGCCTATAGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                             <  2:832898/149‑1 (MQ=255)
                 cGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGc                                                                                                                            <  1:316926/149‑1 (MQ=255)
                 cGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGc                                                                                                                            <  2:1192815/149‑1 (MQ=255)
                 cGCTGGCCCGTGGAAAACCTGTTTACGCAAATCTATGTACTATTGCTATTCGGAGTATCCCCGCGGGCGCGGGGAACACGGGCGCTCGGAATACACAGCCCTGTATCTGCTCGGTTTATCCCCGCTGACGCGTGGAACACTGGCTCTGc                                                                                                                            <  1:939165/149‑1 (MQ=255)
                   cTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCaa                                                                                                                          <  1:886206/149‑1 (MQ=255)
                      gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcag                                                                                                                       <  1:281958/149‑1 (MQ=255)
                        gcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagca                                                                                                                     <  2:283218/149‑1 (MQ=255)
                        gcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagca                                                                                                                     <  2:935415/149‑1 (MQ=255)
                          ggggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACACGGGGAACACTGGCTCTGCAACagcagc                                                                                                                   <  2:848668/149‑1 (MQ=255)
                             gAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAcc                                                                                                                >  1:434237/1‑149 (MQ=255)
                                cacCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCAt                                                                                                             <  1:347970/149‑1 (MQ=255)
                                 acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATg                                                                                                            >  1:349698/1‑149 (MQ=255)
                                 acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATg                                                                                                            >  1:339343/1‑149 (MQ=255)
                                   cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTTCTCTGTTTATCCCCGCTGACGCGGGGAAAAATGGCGATGCAACAGCAGCACACAGGAc                                                                                                          >  2:812128/1‑149 (MQ=255)
                                   cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATgcc                                                                                                          >  1:94138/1‑147 (MQ=255)
                                   cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCAGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAc                                                                                                          >  1:891653/1‑149 (MQ=255)
                                    tGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGTGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAcc                                                                                                         <  2:907281/149‑1 (MQ=255)
                                    tGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCACCAGCAGCACCCATGAcc                                                                                                         >  1:766405/1‑149 (MQ=255)
                                         tCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCAGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAAAACGTc                                                                                                    >  2:589278/1‑149 (MQ=255)
                                           gCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACCCTGGCGCTGCAACCGCAGCACCCATGCCCACGTCgg                                                                                                  >  1:571209/1‑149 (MQ=255)
                                           gCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTAAAACAGCATCACACATGACAACGTCgg                                                                                                  >  1:976250/1‑149 (MQ=255)
                                            ctaaTCTATGTACTATTGCTATTCGGTTTATCCCCGCTTGCGCGGGGAACAGGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTTCAACAGCAGCACCCATGACCACGTCGGt                                                                                                 <  1:93484/147‑1 (MQ=255)
                                             aaaTCTATGGACTATTGCTATTAGGTTTATCCCCGCTGCCGCGGGGAACACGTGCGCACGGAATAGAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGtt                                                                                                <  1:1187632/149‑1 (MQ=255)
                                               aTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACAACGTCGGTTTa                                                                                              >  1:536974/1‑149 (MQ=255)
                                                 cTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATc                                                                                            <  1:975011/149‑1 (MQ=255)
                                                    tGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATcccc                                                                                         <  1:714592/149‑1 (MQ=255)
                                                     ggACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACGCATGACCACGTCGGTTTATCCCCg                                                                                        >  2:399294/1‑149 (MQ=255)
                                                        cTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCTACAGCATCACCCATGTCCACGTCGGTTTATCCCCGCTg                                                                                     <  1:120080/149‑1 (MQ=255)
                                                           ttGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGAACCCATGACCACGTCGGTTTATCACCGCTGgcg                                                                                  >  1:165976/1‑149 (MQ=255)
                                                             gCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGTGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGcgcg                                                                                >  1:682041/1‑149 (MQ=255)
                                                              ctcttCGGTTTATCCCCGGTGGCGCGGGGAACACGGGCTCACGGAATAGAAAGCCGTGTATCTGCTCGGTTTTTACCCGCTGACTCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCgg                                                                               <  2:668810/146‑1 (MQ=255)
                                                               tATGCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                              <  1:588653/149‑1 (MQ=255)
                                                                 ttCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggcgc                                                                            >  1:189451/1‑146 (MQ=255)
                                                                    ggTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACGACGTCGGTTTATCCCCGCTGGCGCGGGGAaca                                                                         >  1:374691/1‑149 (MQ=255)
                                                                    ggTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACACATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAaca                                                                         >  1:278638/1‑149 (MQ=255)
                                                                     gTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAacac                                                                        >  2:338386/1‑149 (MQ=255)
                                                                     gTTTATCCCCGCAGGCGCGGGGAACACGGGCGCACGGGATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGTGAACACTCGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAacac                                                                        <  1:389217/149‑1 (MQ=255)
                                                                      tttATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCGGCCGGGGGGAACACTGGCGCTGCAACAGCAGCACCGATGACCACGTCGGTGTATCCCCGCTGGCGCGGGGAacagg                                                                       >  1:788739/1‑147 (MQ=255)
                                                                            cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACCCGAAATGc                                                                 >  2:862809/1‑149 (MQ=255)
                                                                            cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGc                                                                 <  1:268646/149‑1 (MQ=255)
                                                                             ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCGGGCGCGGGGAACACGAAATGcg                                                                >  1:1083399/1‑148 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                               <  1:161339/149‑1 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                               <  2:154866/149‑1 (MQ=255)
                                                                                     gcgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGTCTCTGCATCAGCAGCACCGATGACCTCTTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCg                                                        <  1:37487/149‑1 (MQ=255)
                                                                                     gcgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGATCACGAAATGCTGGTGAGCg                                                        >  2:161104/1‑149 (MQ=255)
                                                                                      cgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGCACACGAAATGCGGGTGCGCGt                                                       >  1:1173234/1‑149 (MQ=255)
                                                                                         ggAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTGTCTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTaa                                                    <  2:915442/149‑1 (MQ=255)
                                                                                           aaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATg                                                  <  1:86500/149‑1 (MQ=255)
                                                                                           aaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATg                                                  <  1:334708/149‑1 (MQ=255)
                                                                                           aaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATAACCACGTCTGTTTATCCCCGCTGGCGCGGGGAACACGAAATTCTGTTGAGGGTTAATg                                                  >  2:809034/1‑149 (MQ=255)
                                                                                            ccccGGGCCCACCGAAGCCAAAGCCGAGTATCGGCTCGGAGTATCCCCGCTCACGCGGCGAACACTGGCGCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGCCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGc                                                 <  2:806612/146‑1 (MQ=255)
                                                                                             cacGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGcc                                                >  2:640063/1‑149 (MQ=255)
                                                                                              acGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCg                                               <  2:1083399/149‑1 (MQ=255)
                                                                                                  gcgcACGGAATACAAAGCCGTGTATCTGCTCGGTATATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                           <  1:157036/149‑1 (MQ=255)
                                                                                                  gcgcACGGAATACAAAGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                           <  1:315962/149‑1 (MQ=255)
                                                                                                     cACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAacac                                        >  1:394609/1‑149 (MQ=255)
                                                                                                        ggAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACACATGACCACGTCGGTTTATCCCCGCTGGCGCGGCGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAgg                                     >  1:75188/1‑149 (MQ=255)
                                                                                                              cAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACCGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATc                               >  2:869031/1‑149 (MQ=255)
                                                                                                                     gtgtATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTg                        >  2:528381/1‑149 (MQ=255)
                                                                                                                        tATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATCACCACGTCGGTTTATCCCCGCTGGCGCGGCGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGgcc                     >  2:519257/1‑148 (MQ=255)
                                                                                                                          tCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGcgcg                   <  1:1184592/149‑1 (MQ=255)
                                                                                                                            tGCTCGGTTTATCCCCGCTGACGCGGGGAACACAGGCTCTGCAACAGCAGCACCCATGACCACTTCTGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAATCACAGGTTTATCCCCGCTGGCGCggg                 <  2:729231/149‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCCCCGCTGACGCGGGAAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgggg                >  1:898959/1‑149 (MQ=255)
                                                                                                                              cTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGCCCACGTCGGTTTATCCCCGCTGGGGCGGGGAACACGATATGCTGGGGAGCGTTAATGCAGAAAAAACAGGTTTATCCCCGCGTGCGCggggg               >  2:969170/1‑148 (MQ=255)
                                                                                                                              cTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGATTATCCCCGCTGGCGCGGGGa               >  1:251156/1‑149 (MQ=255)
                                                                                                                               tCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGaa              <  2:339343/149‑1 (MQ=255)
                                                                                                                               tCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGaa              >  2:201377/1‑149 (MQ=255)
                                                                                                                               gcGGTTTATCCCCGCTGACGCGGGGAACACTGTCTCTGCAACAGCATCACCGATGACAACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTTTTGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGaa              <  2:533511/148‑1 (MQ=255)
                                                                                                                                  gattATCCCCCGTGACGCGGGGAACAGTGGCTCTGCTTCAGCAGCAACAATGACCTCTGCGGTTTATCCCCGCTGGCTCGGGGAACACGAATTGCTGGTGAGCTTTTATGGCGCAATCACAGGTTTATCCCCGCTGGCGCGGGGAacac           <  2:101004/147‑1 (MQ=255)
                                                                                                                                  gTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCTCGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAacac           >  1:706920/1‑149 (MQ=255)
                                                                                                                                   tttATCCCCGCTGACGCGGGGAACACTTGCTCTGCAACAGCATCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAACTGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAAcaca          <  1:1084107/149‑1 (MQ=255)
                                                                                                                                   tttATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACGCGAAATGCTGGTGCGCGTTAATGCCGCATACCCAGGGTTGTGCCCGCTGGCggggggggccca          >  1:200740/1‑143 (MQ=255)
                                                                                                                                     tATCCCCGCTGCCGCGGGGAACACTGGCTCTCCAACAGCAGCACCCATGACCACGTCGGGTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGACAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAtt        <  2:680659/149‑1 (MQ=255)
                                                                                                                                      aTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTa       >  1:536490/1‑149 (MQ=255)
                                                                                                                                      aTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTAATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGATTATCCACGCTGGCGCGGGGACCACATTa       >  1:5742/1‑149 (MQ=255)
                                                                                                                                       tCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAGTAc      >  1:609430/1‑149 (MQ=255)
                                                                                                                                        ccccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACg     <  2:377371/149‑1 (MQ=255)
                                                                                                                                         cccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGc    <  1:528381/149‑1 (MQ=255)
                                                                                                                                          ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCGGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGcc   >  1:373906/1‑149 (MQ=255)
                                                                                                                                           cGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCt  >  2:633561/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCT  >  NC_000913/2904068‑2904355

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