Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 21.7% intergenic (+200/+207) 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,1910CA21.7% 163.0 / 15.3 89intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (20/47);  new base A (19/3);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 3.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552042‑1552340
                                                                                                                                                     |                                                                                                                                                      
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                        <  1:772128/150‑1 (MQ=31)
  atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                      <  2:742942/150‑1 (MQ=32)
  atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                      <  2:2075898/150‑1 (MQ=32)
      tataCAACGTAGGCTTTGTTCATTCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAATCCCTTCAAATCCAATAGATTGAAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                                                                                                                                  <  1:117598/147‑1 (MQ=255)
       aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                                                                                                                 <  2:763363/150‑1 (MQ=34)
        gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                                <  2:1982757/150‑1 (MQ=34)
         aCAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                               <  2:516470/150‑1 (MQ=21)
         aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                               <  1:2046166/150‑1 (MQ=34)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                             <  1:1187574/150‑1 (MQ=35)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                             <  2:903927/150‑1 (MQ=35)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                                                                                                                                        <  1:1954752/150‑1 (MQ=39)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTaaaa                                                                                                                                        >  2:1958206/1‑150 (MQ=21)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATTCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                       <  2:715844/150‑1 (MQ=255)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAgg                                                                                                                                      >  1:1208267/1‑150 (MQ=18)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCATTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                      <  2:2080900/150‑1 (MQ=255)
                    tttGTTCATGCCGGATGCTGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGAATCAGGCAATGTTTAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                    <  2:277818/150‑1 (MQ=255)
                        ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                                <  2:1906734/150‑1 (MQ=255)
                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                              >  1:1657332/1‑150 (MQ=255)
                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                              <  1:731630/150‑1 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGGGTCAtttta                                                                                                                            >  2:1908082/1‑149 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGTTTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                            <  2:778579/150‑1 (MQ=255)
                            tGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGAGACAtttta                                                                                                                            >  2:1122557/1‑149 (MQ=11)
                                ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATAAGGCAATGATGCGTTTTTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                        >  2:1612272/1‑150 (MQ=255)
                                   ttcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                     <  1:195830/148‑1 (MQ=255)
                                    gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                    <  2:1189440/150‑1 (MQ=255)
                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                   >  1:1775720/1‑150 (MQ=255)
                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTAATATTGAAg                                                                                                                   >  1:10604/1‑150 (MQ=255)
                                          tGAACGCCTTATCCGGAATGAAAACACTTCAAATCCAATAGATTGCTGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                              <  1:454534/150‑1 (MQ=18)
                                           taaCGGCTTATCCGCAATGAAAACCCTTCAAATCCAATATATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCCGTTGAGATCAGTATCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACg                                                                                                             <  1:73627/149‑1 (MQ=17)
                                              cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                          <  1:1312052/150‑1 (MQ=255)
                                               gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                         <  1:443474/150‑1 (MQ=255)
                                                tcTTATCCTGCATGATAACCCTTAAATTCCAATAGATTGCAGTGATCGTGTAGGCCTGATAAGCGTAGCGAATGAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                        <  2:1760370/149‑1 (MQ=18)
                                                ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAAAGTAgga                                                                                                        >  2:1804929/1‑149 (MQ=255)
                                                    aTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg                                                                                                    <  2:1166264/150‑1 (MQ=255)
                                                         gCATGAAAACGCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATGAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCAt                                                                                               <  2:101026/150‑1 (MQ=255)
                                                              aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTAATGcatt                                                                                          >  2:655883/1‑147 (MQ=255)
                                                                 aCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCAtc                                                                                       <  2:1775720/150‑2 (MQ=255)
                                                                  cccTTCTAATTCTATAGATTGCATTGAAAGTGTCGGACTGATAAGCGTAGCGCATGAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                      <  2:1025735/150‑1 (MQ=25)
                                                                        aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATTTTGATATTGTAGAAAAAGACGGCTTAGTACATGCACGATgcggcttt                                                                                >  2:1080761/1‑147 (MQ=17)
                                                                         aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTTTAGACAACGTAGGCTTTATTCATGCCGGATGCAGCGTGa                                                                               >  2:1950453/1‑150 (MQ=21)
                                                                            ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGCGGCGTGAAcc                                                                            >  2:731329/1‑149 (MQ=255)
                                                                               aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                         <  1:1720101/150‑1 (MQ=255)
                                                                                tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                        <  2:526853/150‑1 (MQ=255)
                                                                                   aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCGTCATTTTCATATTGTTGACCACGTAGGCTTTGTTCATGCCGGATCCGGCGTGAACGCCTTATc                                                                     >  1:862282/1‑150 (MQ=255)
                                                                                   aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATTTTAATATTGTAGACAAAGTAGGCTTTGTTCATGCCCGAAGCGTCGTGATCGCCGTaac                                                                     >  2:98319/1‑148 (MQ=18)
                                                                                     tgaaGTGAAAGTGTAGGTATTATAAGAGTATCTCATCAGGCAATTTTTCTTTTTTCATCAGTTTCAAATGGCGCTGTTAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCTTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                  <  2:1349039/147‑1 (MQ=17)
                                                                                       cAGTGAACGTGTAGTCCTGATAATCGGATCTCATCAGGCAATGTTGCTTTTGTAATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                <  2:254920/150‑1 (MQ=18)
                                                                                       cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGc                                                                >  1:1698997/1‑150 (MQ=255)
                                                                                         ttGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                              <  1:783860/149‑1 (MQ=255)
                                                                                         gcgAACGTGTAGGCCTGATAACCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                              <  2:16304/148‑1 (MQ=255)
                                                                                         gTGAACGTGTAGGCATGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTACGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                              <  1:1295421/150‑1 (MQ=255)
                                                                                         ctGAACGAGAAGGCATTATAAGCGTAGCGCATCAGGCAATGTTTCGGTAGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                              <  2:1310647/149‑1 (MQ=18)
                                                                                             aCGTCTAGGCATGATATTCGTATCGAATCAGTTAATGTTGATTTTTTAATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaa                                                          <  2:227638/150‑1 (MQ=17)
                                                                                               gtgtCGGCCTGATAAGAGTATCTAATCATTCAATGTTTAGTGAGTCATCATATTCAACTGTCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGTTGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                        <  1:1259295/150‑1 (MQ=14)
                                                                                                 gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                      <  2:692421/150‑1 (MQ=255)
                                                                                                 gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                      >  2:1707664/1‑150 (MQ=255)
                                                                                                  tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                     <  2:1683115/150‑1 (MQ=255)
                                                                                                  tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccat                                                     >  2:1903878/1‑148 (MQ=255)
                                                                                                    ggcattaTAAGCGTAGCGCTTCAGGCAATGTTGCGTTTGTCATCAGGTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                   <  2:395553/144‑1 (MQ=18)
                                                                                                    ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTTTACATTCAGGATGCGTCTTTAACGTCTAAAACGGCATGAAAAGCCtta                                                   >  2:411496/1‑149 (MQ=21)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGGCATTTTAATATTGAATACAACGTAGGCTTTGTTCATGCCGGAAGCGTCGTTAAAGGCGGAACCCGCATGAAAACCCTTCAAATc                                              >  2:781464/1‑150 (MQ=18)
                                                                                                            aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                           <  2:1131950/150‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGCCTTATCCGGCATGAAAACCCTTCAAATACAAt                                          >  1:2018184/1‑150 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAACTCCAAt                                          >  2:449970/1‑150 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCCGAAGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAACCAAt                                          >  2:2087646/1‑150 (MQ=18)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACACCGAAGGCTTAGTTCATACCGGATCCGGCGTGAAAGACAATCACGGCCTCACAACCCTTCAAATCCAATa                                         >  2:275422/1‑150 (MQ=18)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCAAtg                                         >  2:1174465/1‑149 (MQ=255)
                                                                                                              ccGTATCGAATCAGTAAATGTTGATTTAGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         <  2:1896033/149‑1 (MQ=14)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAACACCCTTCAAATCCAATAGa                                       >  1:732173/1‑150 (MQ=255)
                                                                                                                 tAGCGCATCAGGAAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTACAAATCCAATAGAt                                      <  2:933583/150‑1 (MQ=255)
                                                                                                                   tcgcATCAGTAAATGTTGCTTTTGTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                    <  1:1877360/149‑1 (MQ=255)
                                                                                                                    ggcATCAGTCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:1144906/149‑1 (MQ=255)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGGGTCATTTTATTATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCCGCTTTAACGACTTATACGGAAAGAAAACCCTTAAAATCAAATAGATAGCatttaaa                            >  1:196817/1‑144 (MQ=17)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGAAGGCTTTGTTCATTCCGGATGCGGCGTGAACGCCGGATCCGCCATCAACACCCTTCAAAAACAATAGATTCCAGTTAAc                            >  2:1296843/1‑150 (MQ=17)
                                                                                                                               aTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTATGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                        <  2:1657332/150‑1 (MQ=255)
                                                                                                                                tGTTTCTTTTGGCTTCAGGTTCAAATGGCGTTTAAATAGGGGTCTTTTTTATATTGTTGACAAGTTAGGATTTGTTCATGCCGTAGGCGGCGTTAACGCCTTATCCGGCATGAATACCCTTCAAATCCAATAGATTGGAGTGAACGTTTa                       <  2:659144/150‑1 (MQ=2)
                                                                                                                                  ttGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACGTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTTAACGTGTatg                     >  1:150366/1‑148 (MQ=17)
                                                                                                                                    gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                   >  2:680764/1‑150 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAAGTGTAGGcc                   >  2:651437/1‑150 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAACGACTTATCACGCATGAAAACCCTTCAAATCCAAAAGATTTCATTGAACGTGTAGGcc                   >  1:243591/1‑150 (MQ=14)
                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                  >  2:747310/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  1:898668/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  2:1685670/1‑150 (MQ=255)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCTGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAACTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  1:1263743/1‑150 (MQ=255)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAACCCCTTCAAATCCAATACATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  1:155583/1‑150 (MQ=255)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:1333524/1‑150 (MQ=39)
                                                                                                                                                 gTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      <  2:788237/150‑1 (MQ=17)
                                                                                                                                                  tttCAAATGGCGCTGTAAAATGCGTTATTTTCATATTGTTGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg     <  1:1612272/150‑1 (MQ=255)
                                                                                                                                                  tttCAAATGGCGCTGTAAAAGGCGTCATTTTCATTTTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg     <  1:1295805/150‑1 (MQ=25)
                                                                                                                                                   ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgg    >  1:317938/1‑149 (MQ=21)
                                                                                                                                                     cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:189624/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552042‑1552340

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: