Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,904,212 G→A 100% 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,2120GA94.9% 198.7 / ‑2.5 91intergenic (+838/+535)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (40/43);  minor base C (3/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.

TCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGC  >  NC_000913/2904066‑2904353
                                                                                                                                                  |                                                                                                                                             
tCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTATCGCACATCTATGGACTATTTCTATTCGGTTTATCCACTCTGGCGCGGGGAACAAGGGAGCACGGAATACAATGAAGTGTATATGCTAGGGTTATGCCCGCTGacg                                                                                                                                             >  2:477105/1‑149 (MQ=255)
   ggCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                          <  2:694557/149‑1 (MQ=255)
   ggCTCCGGTGTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTCTGGACTATTGCTATTCGGTTTATCCCCGCGGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                          <  1:226246/149‑1 (MQ=255)
     cTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAAGACGGGCGCACGGAATACAAAGCCGAGTATCGGCGCGGTTTATCCCCGCTGACGCgggg                                                                                                                                        >  2:949253/1‑149 (MQ=255)
        cGGTTTATCCCCGCTGGCGCGGGGAACACGTGTTATCGCAAATCTATGGACTGTTGGTATTCGGTTTATCCCCGCTGGGGCGGGGAACACGGGAGAAGGGAATAGAAAGCAGTGGTTCTGGTGGGGTTTTCGCGGCTGAGGCggggggc                                                                                                                                     >  2:817518/1‑146 (MQ=255)
         ggTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGataa                                                                                                                                    <  2:145946/149‑4 (MQ=255)
          ggtgttCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCACGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGCacac                                                                                                                                   >  2:35722/6‑149 (MQ=255)
          gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACGATTGCTATGCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATGCTAAGCCGTGTATGTGGTCGGTTTATCGCCGCTGGCGCGGGGTacac                                                                                                                                   >  2:684630/1‑149 (MQ=255)
          gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATATATGGACTATTGCTATTCGGTTTAACCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAACGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAacac                                                                                                                                   >  2:788979/1‑149 (MQ=255)
           tttATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACt                                                                                                                                  <  2:455538/149‑1 (MQ=255)
                 cccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctct                                                                                                                            <  1:140535/149‑1 (MQ=255)
                  ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                           <  1:543380/149‑1 (MQ=255)
                     cTGGCGCGGGGAACACCTGTTTTCGCAAATCTATCGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCaa                                                                                                                        <  1:197131/149‑1 (MQ=255)
                      tGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAc                                                                                                                       <  2:794276/149‑1 (MQ=255)
                       ggCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAca                                                                                                                      <  1:554847/149‑1 (MQ=255)
                          gcgGGGAACACCTGTTTTCGCAAATCTATGGACTTTTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagcc                                                                                                                   >  2:321095/1‑148 (MQ=255)
                             gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAgcagca                                                                                                                <  2:390807/149‑1 (MQ=255)
                             gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAgcagca                                                                                                                >  2:873597/1‑149 (MQ=255)
                                aaCACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAccc                                                                                                             >  2:424510/1‑149 (MQ=255)
                                aaCACCTCTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCAGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAccc                                                                                                             <  1:1017994/149‑1 (MQ=255)
                                  tacCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCAt                                                                                                           <  2:251289/148‑1 (MQ=255)
                                   acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGACCCCATg                                                                                                          >  1:995991/1‑149 (MQ=255)
                                   acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGGATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCAGCACCAGCACCACCCaag                                                                                                          >  1:992063/1‑147 (MQ=255)
                                     cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAc                                                                                                        <  1:939100/149‑1 (MQ=255)
                                     cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAc                                                                                                        <  2:356485/149‑1 (MQ=255)
                                     cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAc                                                                                                        <  2:271554/149‑1 (MQ=255)
                                      tGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACTCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAcc                                                                                                       >  1:781457/1‑149 (MQ=255)
                                       gTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAACCATGAccc                                                                                                      >  2:38468/1‑148 (MQ=255)
                                          ttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGt                                                                                                   >  1:264962/1‑149 (MQ=255)
                                            cGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCg                                                                                                 <  1:910037/149‑1 (MQ=255)
                                             gCAAATCTTTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCgg                                                                                                >  2:811611/1‑149 (MQ=255)
                                                aaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGttt                                                                                             <  1:532434/149‑1 (MQ=255)
                                                 aTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACTGAAAACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAGCACTGGCTCTGCAACAGCAGCACCGATGACCAAGTCGGTTTa                                                                                            >  1:348924/1‑149 (MQ=255)
                                                  tCTTTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTAt                                                                                           <  2:14452/149‑1 (MQ=255)
                                                     aTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATccc                                                                                        <  2:956447/149‑1 (MQ=255)
                                                      tGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTCTCTGCTCGGTTTATCCCCGCTCACGAGGGGAACACTGGCACTGCATCAGCAGCGCCCATGCCGACGTCGGTTTAGcccc                                                                                       >  2:868479/1‑149 (MQ=255)
                                                         aCTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCt                                                                                    <  1:255371/149‑1 (MQ=255)
                                                          cTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAAAACAAAGCCGTGTATCTGCTCGGGTTAACCCCGGTGCCGCGGGGAACACTGGCTCTGCAAAAACATCACCCATGAGGACGTCGGTTTATCCCCGCTg                                                                                   >  1:608463/1‑149 (MQ=255)
                                                                   ttCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggggg                                                                          >  1:537507/1‑148 (MQ=255)
                                                                   ttCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCAGCAACAGCAGCTCACATGACCACGTCGGTTTATCCCCGCTGGCGCggggg                                                                          >  1:54690/1‑148 (MQ=255)
                                                                    tCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGaa                                                                         >  1:168231/1‑149 (MQ=255)
                                                                      ggTTTATCCCCGCTGGCGCGGGGAACCCGGGCGCACGTAATACAAAGCCGTGTATCTGCTCGGTATCTCCCCGCAGACGCGGGGAACACAGGCTCTGCAACACAAGCACACATGAACAAGACGGTTTATAAACGCTGTcgcgtgtacca                                                                       >  2:210964/1‑142 (MQ=255)
                                                                       gTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAgcac                                                                      >  1:306581/1‑149 (MQ=255)
                                                                           aTCCCCGCTGGCGCGGGGAACACGGGCGCACGCAATACAAAGGCGTGGATCTGCTAGGTTTATCCCCGCTGACACGGCGAACACTGGCTCTGCAACAGCAGCACCCATCACCACGTCGCTATATCCCCCATGGCGCGGGGCACACGaaa                                                                  >  2:67924/1‑149 (MQ=255)
                                                                            tCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAAt                                                                 <  1:308943/149‑1 (MQ=255)
                                                                              cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGc                                                               >  2:1017438/1‑149 (MQ=255)
                                                                              cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACAGGCTCTGCAACAGCAGCACCCAAGACAACGTCGGTTTATCCCCGCAGGCGCGGGGAACACCAAATGc                                                               >  1:989110/1‑149 (MQ=255)
                                                                              cccGCTGGCGCGGGGAACACCGGCGCACGGAATACATAGCCGTGTATCTGCTCGGGTTATCCCCGCTGACGCGGGGAACACTGGCTCTGAAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGc                                                               <  1:681825/149‑1 (MQ=255)
                                                                               ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGACCCCGAAATGcg                                                              >  1:178863/1‑148 (MQ=255)
                                                                               ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATACCCGCTGGCGAGGGGACCACGAAATgat                                                              >  1:166040/1‑147 (MQ=255)
                                                                                 gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTTTATCTGCTCGGTTTATCCCCGCTGACGCTGGGAACACTGGCTCTGCAACAGCAGCACACATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTgg                                                            <  1:648292/149‑1 (MQ=255)
                                                                                  cTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                           <  2:475177/149‑1 (MQ=255)
                                                                                   tGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGCACACTGGCTCTGCAACAGCAGCGCCCATGACCACGTCGGTTTCTCCCCGCTGGCGCGGGGAACACGAAATGCTGGTg                                                          >  1:371596/1‑149 (MQ=255)
                                                                                    ggCGCTGGGAACTCGGGCGCCCGGTATACATAGACGTGTATCTGCTCGGTTTATCCCCGCTGTCGCGGGGAACACTGGCTCTGCAACAGCACCAACCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGa                                                         <  2:549067/149‑1 (MQ=255)
                                                                                        cgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCACCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGt                                                     <  2:983358/149‑1 (MQ=255)
                                                                                          gggAACACGGGCGCACGGACTACAAAGCCGGGTATCTGCGCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACTCGAAATGCTGGTGAGCGTTa                                                   <  1:340007/149‑1 (MQ=255)
                                                                                             aaCACCGGCGCACGGAATACCTCGCCGTGTATCTGCTCGGTTAATCGCCGCTGACGCGGGGAACACTTGCTCTGCTACAGCAGCACCCTTGACCACGTGGGTTTATCCCCGCTGGCGCGGGGAACACGGAATGCTGGTGAGCGTTAATg                                                <  1:88815/149‑1 (MQ=255)
                                                                                               cacGGTCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGcc                                              <  2:11829/149‑1 (MQ=255)
                                                                                                acGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAGCACGAAATTCTGGTGAGCGTTAATGCCg                                             >  1:928745/1‑149 (MQ=255)
                                                                                                    gcgcACGGAATACAAAGCCGTGTAGCTGCTCGGGTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                         <  1:655699/149‑1 (MQ=255)
                                                                                                        aCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAAcaca                                     >  2:599253/1‑149 (MQ=255)
                                                                                                        aCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTCACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCCCAGGCGCGTGGAACACGAAAAGCTGTTGAGCGATAATGCCGCataagca                                     >  1:565545/1‑143 (MQ=255)
                                                                                                         cGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                    <  1:873597/149‑1 (MQ=255)
                                                                                                         cGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                    <  1:120832/149‑1 (MQ=255)
                                                                                                         cGGAATACAAAGCAGTGGTTCGGCTCTGTGTATCCACGCTGACGCGGGGAACTCTGGCGCTGCAAAAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                    <  2:195397/149‑1 (MQ=255)
                                                                                                          ggAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGCGGAACACGAAATGCTGGTGAGCGTTAATGCAGCAAACACAgg                                   >  1:834362/1‑149 (MQ=255)
                                                                                                           taatcgaaaGCCGGTTATCTGCGCTGTTTATCCCCGCTGACTCGGGGAACACTGGCTCTGCTAAAGCAGAACCCATGATCACGGCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTCGTGAGCGTTAATGCCGCAAACACAGGt                                  <  2:93985/143‑1 (MQ=255)
                                                                                                           gAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGTCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGt                                  <  2:644715/149‑1 (MQ=255)
                                                                                                           gAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGt                                  <  1:321095/149‑1 (MQ=255)
                                                                                                             aTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCTGGGCACACGAAATGCATGAGAGAGTTAATCCCGCAAACACAGGttt                                >  1:680949/1‑149 (MQ=255)
                                                                                                              tACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGGGCGTTTATGACGCAAACACAGGTTTa                               >  1:193582/1‑149 (MQ=255)
                                                                                                              tACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGAAGCGGGGAACACTGGCTCTGCAACAGCCGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAAGGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTa                               >  2:834196/1‑149 (MQ=255)
                                                                                                               aCAAATCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCTCCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTAt                              <  1:994896/149‑1 (MQ=255)
                                                                                                                cAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATc                             <  1:69677/149‑1 (MQ=255)
                                                                                                                 aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGTACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                            >  1:348849/1‑149 (MQ=255)
                                                                                                                 aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                            <  2:428217/149‑1 (MQ=255)
                                                                                                                 aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                            >  2:711344/1‑149 (MQ=255)
                                                                                                                 aaaGCCGGGTATCTGCTCGGTATATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                            <  1:480124/149‑1 (MQ=255)
                                                                                                                   aGCCGTGTATCTGCTCGGTTTATCCCCGATGAAGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACTTCGTTTTATCCCCGCTGGAGCGGGGAACAAGAAATGATGGTGATCGTTAATGAATCAAACACAGTTTTATcccc                          >  2:296679/1‑149 (MQ=255)
                                                                                                                       gtgtATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGACCACGAAATGCTGGTGAGCGGTAATGCCGCAAACACAGGGTTATCCCCGCTg                      >  1:210574/1‑149 (MQ=255)
                                                                                                                              tGCTCGGCTTAACCCCGCTGCCGCGGGGATCACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCggg               <  1:684630/149‑1 (MQ=255)
                                                                                                                              tGCTCGGATTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCggg               <  1:788979/149‑1 (MQ=255)
                                                                                                                                cTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgggca             >  2:127064/1‑147 (MQ=255)
                                                                                                                                  cGGCTTATCCCCGCTGACGCGGGGAACACTGGCGCTGCAACAGCCCCACCCCTGACCCCGGCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAac           <  1:948062/149‑1 (MQ=255)
                                                                                                                                   ggTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGCCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCGGGTGAGCGTTTATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAaca          >  2:677360/1‑149 (MQ=255)
                                                                                                                                    gTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAACCACAGGTTTATCCCCGCTGGCGCGGGGAacac         >  1:438897/1‑149 (MQ=255)
                                                                                                                                      ttATCCCCGCTGACGCGGGGAACACTTGCTCTGCAACAGCATCACCCATGACCACGGCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAt       <  2:245150/149‑1 (MQ=255)
                                                                                                                                      ttATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCGCCCATGACCACGTCGGTTTATCCCCGCGGGCGGGGGGAACACGAAATGCTGGTGAGCGTTACTGCCGCACACACAGGTTTATCCCCGCGGGCGCggggcccacct       >  1:919536/1‑142 (MQ=255)
                                                                                                                                      ttATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAt       >  2:97590/1‑149 (MQ=255)
                                                                                                                                       tATCCCCGCAGAAGCTTGGAACACTTGCTCTTCAACAGGAGCACCCAGGACCACGTCGGTTTATCCCCTCTGTCGCGTGTAACACGAAATGCTGGTGAGAGTTAATGCCTGAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAtt      <  2:644259/149‑1 (MQ=255)
                                                                                                                                        aTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTa     <  1:424510/149‑1 (MQ=255)
                                                                                                                                           cccGCGGCCGCGGGGAACACTGGCTCTGCAACATCCGAACCCTTGTCCACGTCTGTTTATCCCCGCTTGCTCGGGGAACAAGAAATGCTGGTTAGCGTTAATGTCTCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGc  <  2:608463/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGC  >  NC_000913/2904066‑2904353

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: