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

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  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)
                                                                                     tgaaGTGAAAGTGTAGGTATTATAAGAGTATCTCATCAGGCAATTTTTCTTTTTTCATCAGTTTCAAATGGCGCTGTTAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCTTGCCGGATGCGGCGTGAACGCCTTATC‑Cg                                                                  <  2:1349039/147‑1 (MQ=17)
                                                                                       cAGTGAACGTGTAGTCCTGATAATCGGATCTCATCAGGCAATGTTGCTTTTGTAATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGc                                                                <  2:254920/150‑1 (MQ=18)
                                                                                       cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATC‑CGGc                                                                >  1:1698997/1‑150 (MQ=255)
                                                                                         ttGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCAt                                                              <  1:783860/149‑1 (MQ=255)
                                                                                         gcgAACGTGTAGGCCTGATAACCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCAt                                                              <  2:16304/148‑1 (MQ=255)
                                                                                         gTGAACGTGTAGGCATGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTACGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCAt                                                              <  1:1295421/150‑1 (MQ=255)
                                                                                         ctGAACGAGAAGGCATTATAAGCGTAGCGCATCAGGCAATGTTTCGGTAGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCAt                                                              <  2:1310647/149‑1 (MQ=18)
                                                                                             aCGTCTAGGCATGATATTCGTATCGAATCAGTTAATGTTGATTTTTTAATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGaaa                                                          <  2:227638/150‑1 (MQ=17)
                                                                                               gtgtCGGCCTGATAAGAGTATCTAATCATTCAATGTTTAGTGAGTCATCATATTCAACTGTCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGTTGCGGCGTGAACGCCTTATC‑CGGCATGAAAAc                                                        <  1:1259295/150‑1 (MQ=14)
                                                                                                 gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAAccc                                                      <  2:692421/150‑1 (MQ=255)
                                                                                                 gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAAccc                                                      >  2:1707664/1‑150 (MQ=255)
                                                                                                  tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCt                                                     <  2:1683115/150‑1 (MQ=255)
                                                                                                  tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAAccat                                                     >  2:1903878/1‑148 (MQ=255)
                                                                                                    ggcattaTAAGCGTAGCGCTTCAGGCAATGTTGCGTTTGTCATCAGGTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTc                                                   <  2:395553/144‑1 (MQ=18)
                                                                                                    ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTTTACATTCAGGATGCGTCTTTAACGTCTAAAA‑CGGCATGAAAAGCCtta                                                   >  2:411496/1‑149 (MQ=21)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGGCATTTTAATATTGAATACAACGTAGGCTTTGTTCATGCCGGAAGCGTCGTTAAAGGCGGAAC‑CCGCATGAAAACCCTTCAAATc                                              >  2:781464/1‑150 (MQ=18)
                                                                                                            aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCaa                                           <  2:1131950/150‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGCCTTATC‑CGGCATGAAAACCCTTCAAATACAAt                                          >  1:2018184/1‑150 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATA‑CGGCATGAAAACCCTTCAACTCCAAt                                          >  2:449970/1‑150 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCCGAAGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAAACCAAt                                          >  2:2087646/1‑150 (MQ=18)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACACCGAAGGCTTAGTTCATACCGGATCCGGCGTGAAAGAC‑AATCACGGCCTCACAACCCTTCAAATCCAATa                                         >  2:275422/1‑150 (MQ=18)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATC‑CGGCATGAAAACCCTTCAAATCCAAtg                                         >  2:1174465/1‑149 (MQ=255)
                                                                                                              ccGTATCGAATCAGTAAATGTTGATTTAGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATa                                         <  2:1896033/149‑1 (MQ=14)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATA‑CGGCATGAACACCCTTCAAATCCAATAGa                                       >  1:732173/1‑150 (MQ=255)
                                                                                                                 tAGCGCATCAGGAAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTACAAATCCAATAGAt                                      <  2:933583/150‑1 (MQ=255)
                                                                                                                   tcgcATCAGTAAATGTTGCTTTTGTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTg                                    <  1:1877360/149‑1 (MQ=255)
                                                                                                                    ggcATCAGTCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:1144906/149‑1 (MQ=255)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGGGTCATTTTATTATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCCGCTTTAACGACTTATA‑CGGAAAGAAAACCCTTAAAATCAAATAGATAGCatttaaa                            >  1:196817/1‑144 (MQ=17)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGAAGGCTTTGTTCATTCCGGATGCGGCGTGAACGCCGGATC‑CGCCATCAACACCCTTCAAAAACAATAGATTCCAGTTAAc                            >  2:1296843/1‑150 (MQ=17)
                                                                                                                               aTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTATGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                        <  2:1657332/150‑1 (MQ=255)
                                                                                                                                tGTTTCTTTTGGCTTCAGGTTCAAATGGCGTTTAAATAGGGGTCTTTTTTATATTGTTGACAAGTTAGGATTTGTTCATGCCGTAGGCGGCGTTAACGCCTTATC‑CGGCATGAATACCCTTCAAATCCAATAGATTGGAGTGAACGTTTa                       <  2:659144/150‑1 (MQ=2)
                                                                                                                                  ttGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACGTATA‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTTAACGTGTatg                     >  1:150366/1‑148 (MQ=17)
                                                                                                                                    gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                   >  2:680764/1‑150 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAAGTGTAGGcc                   >  2:651437/1‑150 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAACGACTTATC‑ACGCATGAAAACCCTTCAAATCCAAAAGATTTCATTGAACGTGTAGGcc                   >  1:243591/1‑150 (MQ=14)
                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                  >  2:747310/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  1:898668/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  2:1685670/1‑150 (MQ=255)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCTGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAACTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  1:1263743/1‑150 (MQ=255)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAACCCCTTCAAATCCAATACATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  1:155583/1‑150 (MQ=255)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:1333524/1‑150 (MQ=39)
                                                                                                                                                 gTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      <  2:788237/150‑1 (MQ=17)
                                                                                                                                                  tttCAAATGGCGCTGTAAAATGCGTTATTTTCATATTGTTGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg     <  1:1612272/150‑1 (MQ=255)
                                                                                                                                                  tttCAAATGGCGCTGTAAAAGGCGTCATTTTCATTTTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg     <  1:1295805/150‑1 (MQ=25)
                                                                                                                                                   ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgg    >  1:317938/1‑149 (MQ=21)
                                                                                                                                                     cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:189624/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATC‑CGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  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: ‑