Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,027 A→T 6.0% intergenic (+36/+371) 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,0270AT 6.0% 244.0 / 2.9 81intergenic (+36/+371)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base T (2/3);  total (37/44)
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

GGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCA  >  NC_000913/1551881‑1552168
                                                                                                                                                  |                                                                                                                                             
gggATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTaaaa                                                                                                                                             <  3:210949/149‑1 (MQ=255)
 ggATGATCGGAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAg                                                                                                                                            >  8:319504/1‑149 (MQ=255)
 ggATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGAATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGCATATAATAACATATCAACTATAAGCTGAAAATAGCGCAGTaaaaa                                                                                                                                            >  4:223730/1‑148 (MQ=255)
  gatgatCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAgg                                                                                                                                           >  3:328606/1‑149 (MQ=255)
      atCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGgcgc                                                                                                                                        <  3:304767/148‑1 (MQ=255)
         gaagggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAtt                                                                                                                                    <  6:228920/149‑1 (MQ=255)
         gaagggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAtt                                                                                                                                    >  2:362459/1‑149 (MQ=255)
          aagggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAttt                                                                                                                                   >  4:397457/1‑149 (MQ=255)
          aagggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAttt                                                                                                                                   <  1:324475/149‑1 (MQ=255)
            gggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCTTTTACAGCGCCATTTTc                                                                                                                                 <  3:93200/149‑1 (MQ=255)
            gggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCTTTTACAGCGCCATTTTc                                                                                                                                 <  7:66405/149‑1 (MQ=255)
               aaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCata                                                                                                                              <  3:385004/149‑1 (MQ=255)
                  gTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCa                                                                                                                                <  1:587/144‑1 (MQ=255)
                  gTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCa                                                                                                                                >  2:587/1‑144 (MQ=255)
                   tAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGt                                                                                                                          <  6:239365/149‑1 (MQ=255)
                    aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                           >  4:243257/1‑147 (MQ=255)
                      gtcgtcGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGa                                                                                                                       <  7:455810/149‑1 (MQ=255)
                            gCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGCAGCGAAAAAAGAGAGTGAAGAGGGGCTAGAGTCGCTTTTGAACTGTAAGCTGAAAATGGCGCGGTATAAGGCGGCAGTTTCATATTGTGGagcggat                                                                                                                 >  4:105925/1‑143 (MQ=255)
                            gCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                           <  3:105925/139‑1 (MQ=255)
                                      gaaACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                       >  7:331495/1‑149 (MQ=255)
                                         acacCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCAt                                                                                                    <  7:161336/149‑1 (MQ=255)
                                          cacCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCTTTTACAGCGCCATTTTCagc                                                                                                                              <  2:374625/122‑3 (MQ=38)
                                          cacCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCTTTTACAGCGCCATTTTCagc                                                                                                                              >  1:374625/1‑120 (MQ=38)
                                               tCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                              >  3:358911/1‑149 (MQ=255)
                                                 cgcgCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                             <  8:331495/148‑1 (MQ=255)
                                                   cgcTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                          <  6:17762/149‑1 (MQ=255)
                                                         tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCTGATGCGGCGTg                                                                                     >  6:14641/1‑148 (MQ=255)
                                                         tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGttc                                                                                    <  4:358911/149‑3 (MQ=255)
                                                         tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCATccggca                                                                                       <  2:217257/146‑2 (MQ=255)
                                                         tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCATccggca                                                                                       >  1:217257/1‑145 (MQ=255)
                                                           tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                   <  4:322222/148‑1 (MQ=255)
                                                           tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                   <  3:243257/148‑1 (MQ=255)
                                                              tGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                               <  7:385487/149‑1 (MQ=255)
                                                                  aTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                         <  8:329934/135‑1 (MQ=255)
                                                                  aTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                         >  7:329934/1‑135 (MQ=255)
                                                                     gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                        <  8:251533/149‑1 (MQ=255)
                                                                     gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCg                                                                        >  3:411378/1‑149 (MQ=25)
                                                                        aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                     <  1:362459/149‑1 (MQ=255)
                                                                        aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                     <  4:411378/149‑1 (MQ=25)
                                                                          ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                      >  7:212586/1‑146 (MQ=255)
                                                                          ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                      <  8:212586/146‑1 (MQ=255)
                                                                           gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCata                                                                                                                              <  5:162478/89‑1 (MQ=255)
                                                                           gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCata                                                                                                                              >  6:162478/1‑89 (MQ=255)
                                                                             aGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaa                                                                 >  4:85231/1‑148 (MQ=255)
                                                                              gAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGATAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                               >  1:184724/1‑149 (MQ=255)
                                                                                 gCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                            <  1:424027/149‑1 (MQ=255)
                                                                                    aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAg                                                                                                               <  8:382890/95‑1 (MQ=255)
                                                                                    aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAg                                                                                                               >  7:382890/1‑95 (MQ=255)
                                                                                    aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                          >  3:141451/1‑148 (MQ=255)
                                                                                    aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                          >  3:147611/1‑148 (MQ=255)
                                                                                      aaaaGAGAGTGAAGAAGGTATATAATCTCTTTTGAACTTTAAGCTGATAATGGCGCTGTAAATGTCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                         <  6:423924/131‑1 (MQ=255)
                                                                                      aaaaGAGAGTGAAGAAGGGATATAAACTCTTTTGACCTATAAGCTGAAAATGGCGCAGTAAAAGGCGCCATTTTCATATTGTAGACAACGAAGGATTTGTTCATGCCGGATGAGGCGTGAACGCCTTATcc                                                                         >  5:423924/1‑131 (MQ=21)
                                                                                          gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                                   <  7:319504/149‑1 (MQ=255)
                                                                                            gagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                 <  3:85231/149‑1 (MQ=255)
                                                                                                     agGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        <  3:223730/149‑1 (MQ=255)
                                                                                                     agGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        >  3:222677/1‑149 (MQ=255)
                                                                                                       ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                         >  5:23390/1‑114 (MQ=255)
                                                                                                       ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                         <  6:23390/114‑1 (MQ=255)
                                                                                                       ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                      >  4:146114/1‑149 (MQ=255)
                                                                                                         atatAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                     >  4:174060/1‑148 (MQ=255)
                                                                                                          tataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc                                                                                                  >  2:85878/1‑86 (MQ=255)
                                                                                                          tataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc                                                                                                  <  1:85878/86‑1 (MQ=255)
                                                                                                             aaTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                >  4:111296/1‑149 (MQ=255)
                                                                                                                ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtg                               >  1:358305/1‑147 (MQ=255)
                                                                                                                 tctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                            >  2:200123/1‑149 (MQ=255)
                                                                                                                   ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                             <  7:98527/146‑1 (MQ=255)
                                                                                                                   ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                          >  7:270068/1‑149 (MQ=255)
                                                                                                                    tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                       >  2:209066/1‑119 (MQ=255)
                                                                                                                       gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        >  4:198126/1‑131 (MQ=255)
                                                                                                                       gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        <  3:198126/131‑1 (MQ=255)
                                                                                                                        aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                      <  2:358305/148‑1 (MQ=255)
                                                                                                                        aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                     <  3:397457/149‑1 (MQ=255)
                                                                                                                           tttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                  >  1:204677/1‑149 (MQ=40)
                                                                                                                               aGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa              <  5:14641/149‑1 (MQ=37)
                                                                                                                                 cTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc            <  3:174060/149‑1 (MQ=35)
                                                                                                                                 cTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc            <  1:313148/149‑1 (MQ=35)
                                                                                                                                  tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg           >  7:371510/1‑149 (MQ=34)
                                                                                                                                   gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        >  2:131774/1‑119 (MQ=34)
                                                                                                                                   gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                        <  1:131774/119‑1 (MQ=34)
                                                                                                                                   gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg           >  3:297774/1‑148 (MQ=33)
                                                                                                                                    aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa         <  1:200123/149‑1 (MQ=32)
                                                                                                                                        tGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGATAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg     <  2:184724/149‑1 (MQ=17)
                                                                                                                                           ctcTGTAAAAGGCGCCATTTTCATGTTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGCCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa  <  4:315752/147‑1 (MQ=12)
                                                                                                                                                  |                                                                                                                                             
GGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCA  >  NC_000913/1551881‑1552168

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