Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,085 G→T 6.2% intergenic (+94/+313) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,0850GT 6.2% 486.5 / 17.8 145intergenic (+94/+313)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (45/91);  new base T (3/6);  total (48/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 148 total aligned reads displayed.

GGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACG  >  NC_000913/1551936‑1552234
                                                                                                                                                     |                                                                                                                                                     
ggTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                                                                                       <  1:23331/150‑1 (MQ=255)
ggTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                                                                                       <  1:22998/150‑1 (MQ=255)
 gTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                      <  1:400386/150‑1 (MQ=255)
 gTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                      <  1:480279/150‑1 (MQ=255)
 gTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                      <  1:1688243/150‑1 (MQ=255)
   aTCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                                                                                    <  2:353945/150‑1 (MQ=255)
   aTCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                                                                                    <  2:2345354/150‑1 (MQ=255)
   aTCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                                                                                    <  2:473378/150‑1 (MQ=255)
    tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATAAAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                   <  2:828009/150‑1 (MQ=255)
     cGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                                  >  2:804670/1‑150 (MQ=255)
      gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                 <  2:508557/150‑1 (MQ=255)
      gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCATCCGGCATGAAcaaa                                                                                                                                                 <  1:2568686/150‑4 (MQ=37)
      gTGAAATCGAGAACGCAGAAGCGAAAAAAGAGAGTGAAGAAGTGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                 <  2:2886889/150‑1 (MQ=255)
       tGAACTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAAGGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                <  1:89460/150‑1 (MQ=255)
       tGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGccg                                                                                                                                                >  1:2611826/1‑149 (MQ=255)
          aaTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                             >  1:2175322/1‑150 (MQ=255)
            tCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                           <  2:3062903/150‑1 (MQ=37)
            tCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                           <  2:2611826/150‑1 (MQ=37)
             cGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  2:2027597/150‑1 (MQ=255)
              gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  2:2727574/150‑1 (MQ=25)
              gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAAGCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         <  1:1480909/150‑1 (MQ=255)
               agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                        <  1:642469/150‑1 (MQ=255)
               agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                        <  1:642536/150‑1 (MQ=255)
                gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                       <  1:535310/150‑1 (MQ=255)
                gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                       <  1:2603990/150‑1 (MQ=255)
                 aaGGCAGAAGCGAAATAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                                                                                      >  1:827081/1‑150 (MQ=255)
                 aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                      <  1:2918105/150‑1 (MQ=25)
                   ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                                                                                                    >  2:1538104/1‑150 (MQ=255)
                        aaGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                                                                                               >  2:271942/1‑150 (MQ=255)
                        aaGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                                                                                               <  2:1132232/150‑1 (MQ=255)
                                  agagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                                                                                                     >  2:1591796/1‑150 (MQ=255)
                                   gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCa                                                                                                                    <  2:1856123/150‑1 (MQ=255)
                                   gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCa                                                                                                                    <  1:2829758/150‑1 (MQ=255)
                                    agagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                                                                                   <  2:1499164/150‑1 (MQ=255)
                                    agagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                                                                                   <  1:2382794/150‑1 (MQ=255)
                                     gagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                                                                                                  <  2:1666681/150‑1 (MQ=255)
                                      agTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                                                                                                 >  1:335368/1‑150 (MQ=255)
                                           agaagGGATATAATCTCTTTTGAACTTTACGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                                                                                            <  1:1404693/150‑1 (MQ=255)
                                           agaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                                                                                            <  1:2346883/150‑1 (MQ=255)
                                                 gATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGa                                                                                                      <  1:1449051/150‑1 (MQ=255)
                                                   tataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                                                                                                    >  1:1471471/1‑150 (MQ=255)
                                                        tctcTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                                                                                               >  2:689164/1‑150 (MQ=255)
                                                         ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                              >  2:361544/1‑150 (MQ=255)
                                                            ttttGAACTTTTAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                           <  2:2612369/150‑1 (MQ=255)
                                                            ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                           <  2:2612365/150‑1 (MQ=255)
                                                            ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                           >  2:2197524/1‑150 (MQ=255)
                                                                    tttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                   >  1:2879952/1‑150 (MQ=40)
                                                                    tatcaGCTTAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                   <  2:1902557/146‑1 (MQ=0)
                                                                     ttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                                  <  2:686643/150‑1 (MQ=39)
                                                                      tAAGCTTAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt                                                                                 <  1:1658595/150‑1 (MQ=25)
                                                                           tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            >  2:3067486/1‑150 (MQ=34)
                                                                           tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            <  1:2677414/150‑1 (MQ=34)
                                                                           tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            <  2:2347198/150‑1 (MQ=34)
                                                                           tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            <  1:1192694/150‑1 (MQ=34)
                                                                           tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                            <  2:1624847/150‑1 (MQ=34)
                                                                            gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa                                                                           <  1:804670/150‑1 (MQ=33)
                                                                            gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa                                                                           >  1:1792535/1‑150 (MQ=33)
                                                                            gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa                                                                           <  1:2398359/150‑1 (MQ=33)
                                                                             aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGa                                                                                                      >  2:2437321/1‑122 (MQ=33)
                                                                              aaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCaac                                                                         <  1:3067486/150‑3 (MQ=17)
                                                                              aaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAGGAAAACCCTTCAAATCCAATAGATTGCAGTGACCGGGTAGGCCTGATAAGCGTCGCGCATc                                                                         >  1:2663904/1‑150 (MQ=255)
                                                                                 tGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                      <  2:1820577/150‑1 (MQ=32)
                                                                                       tGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGt                                                                >  2:2896742/1‑150 (MQ=32)
                                                                                         tAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTg                                                              <  1:962100/150‑1 (MQ=32)
                                                                                          aaaaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                             <  1:1338924/150‑1 (MQ=32)
                                                                                          aaaaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                             <  1:1043525/150‑1 (MQ=32)
                                                                                            aaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGt                                                           <  2:888039/150‑1 (MQ=32)
                                                                                            aaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGt                                                           <  1:2138906/150‑1 (MQ=32)
                                                                                            aaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGt                                                           <  1:1249510/150‑1 (MQ=32)
                                                                                                cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                                                                                              >  1:2420062/1‑79 (MQ=34)
                                                                                                cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGt                                                       <  2:2445534/150‑1 (MQ=32)
                                                                                                cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAGAAGCGTAGCGCATCAGGCAATGTTGCGTTTGt                                                       <  2:1078316/150‑1 (MQ=17)
                                                                                                 gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAAGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtc                                                      <  1:1001980/150‑1 (MQ=14)
                                                                                                          atatTGTAGACCCCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAGTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                             <  1:918242/150‑1 (MQ=9)
                                                                                                          atatTGTAGACACCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                             <  1:917847/150‑1 (MQ=14)
                                                                                                             ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTGCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATACGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                          >  2:2976178/1‑150 (MQ=18)
                                                                                                                 aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                      <  2:2171936/150‑1 (MQ=34)
                                                                                                                 aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                      >  1:1010534/1‑150 (MQ=21)
                                                                                                                   aCAACGTAGGCTTTGTTCATGCCGGGTGCGGCGTGTACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGTACGTGTAGGCCAGATAAGCGTAGCGCATCAGGCCATGTTGCGTGTGTCATCAGTTTGAAATGgcgc                                    >  1:2916615/1‑150 (MQ=11)
                                                                                                                   aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                    <  2:2992964/150‑1 (MQ=34)
                                                                                                                      aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTTAAATGGCGCTGt                                 >  1:1782818/1‑150 (MQ=14)
                                                                                                                      aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                 <  1:3113642/150‑1 (MQ=35)
                                                                                                                            gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                           <  1:2877973/150‑1 (MQ=40)
                                                                                                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                          <  2:1719861/150‑1 (MQ=255)
                                                                                                                               ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGt                        <  2:85446/150‑1 (MQ=255)
                                                                                                                                    cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                   <  2:2736820/150‑1 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                <  1:2896742/150‑1 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                <  1:2543904/150‑1 (MQ=255)
                                                                                                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGTTTGCAGTGAACGTGTAGGCCTGTTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGGGCCAttttgtt               >  2:2766592/1‑147 (MQ=11)
                                                                                                                                            ctGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg           <  1:1378254/149‑1 (MQ=255)
                                                                                                                                             tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCTATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGGTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt          <  1:2073043/150‑1 (MQ=255)
                                                                                                                                             tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt          <  2:1321466/150‑1 (MQ=255)
                                                                                                                                              gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa         <  1:115989/150‑1 (MQ=255)
                                                                                                                                              gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa         <  1:1894389/150‑1 (MQ=255)
                                                                                                                                              gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa         <  1:1271930/150‑1 (MQ=255)
                                                                                                                                              gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa         <  2:898350/150‑1 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGGGTGCGTGTGTCGTCGGTTACAAGTGGCGCGGTAATAGTCGTCCTCTTCATAGTGTGGa       >  2:1881862/1‑150 (MQ=11)
                                                                                                                                                 gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGGACAAGGCGGCATTGTCCTAGTggcgac      >  2:920796/1‑145 (MQ=14)
                                                                                                                                                   gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACaa    <  2:2002788/150‑1 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCATGAAAACCCTACAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACg  >  1:1992396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACG  >  NC_000913/1551936‑1552234

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