BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,552,191	0	C	A	27.6%	170.0 / 20.0	98	intergenic (+200/+207)	fdnI/yddM	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand): ref base C (19/52); new base A (21/6); total (40/58)
Fisher's exact test for biased strand distribution p-value = 6.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 101 total aligned reads displayed.

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552042‑1552340

atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                       <  1:448804/150‑1 (MQ=31)

  atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                     >  2:1477439/1‑150 (MQ=32)

     gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                                                                                                                                  <  2:1959502/150‑1 (MQ=34)

        gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                               <  1:1485306/150‑1 (MQ=34)

         aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                              <  1:510197/150‑1 (MQ=34)

          cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCt                                                                                                                                             <  2:1393265/150‑1 (MQ=34)

           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  1:1660338/150‑1 (MQ=35)

           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  1:399444/150‑1 (MQ=35)

           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  2:1362118/150‑1 (MQ=35)

           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  1:1494814/150‑1 (MQ=35)

               tAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTAGCGTTTGTCATCACTTAAAAATGGAGCTGTaaa                                                                                                                                        >  2:1520651/1‑150 (MQ=11)

                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                                                                                                                                       >  1:532529/1‑150 (MQ=39)

                    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                   >  1:213556/1‑150 (MQ=255)

                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                          <  2:1841321/150‑1 (MQ=255)

                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                          <  1:930463/150‑1 (MQ=255)

                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTCATCAGTTAAAAATGGCTCTGTAAAATGGGGCAttttatt                                                                                                                         >  1:357667/1‑147 (MQ=11)

                               cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCata                                                                                                                        <  2:788915/150‑1 (MQ=255)

                                ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                       >  2:1158952/1‑150 (MQ=255)

                                   tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGATCGTGGAGGCCTGATAAGCGTAGCGCATCAGGCAATGTAGCGTTAGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                    <  2:1339356/150‑1 (MQ=255)

                                   tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                    <  2:466684/150‑1 (MQ=255)

                                   tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                    <  1:3693/150‑1 (MQ=255)

                                    gcggctTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCATTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                   <  2:836115/150‑1 (MQ=255)

                                    gcggcgTGAACGCATTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                   <  1:1039012/150‑1 (MQ=255)

                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                  <  2:1005719/150‑1 (MQ=255)

                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                  <  2:567323/150‑1 (MQ=255)

                                     cggcgTGAACGCCTTATCCGGAATGAAAACCCTTCAAATCCATTATATTGATTTGTACGTTTAGTCCTTATAAGCGTAGCGATTCATTTAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                  <  2:250156/150‑1 (MQ=17)

                                      ggcgTGAACGCCTTATCCTTCAAGAAAACCCTTCAAATCCAATAGATTTAAGTGAACGAGTAGTCCAGATAATCGTATCGCATCAGTCAATGTTTCGTTAGTCATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                                 <  1:2024163/150‑1 (MQ=17)

                                      ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                                 <  2:424185/150‑1 (MQ=255)

                                           gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACg                                                                                                            <  1:246524/150‑1 (MQ=255)

                                           gAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACg                                                                                                            <  1:1262569/150‑1 (MQ=255)

                                             aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGTCGTCATTTTAATATTGAAGACAACGtc                                                                                                          >  2:1045098/1‑149 (MQ=14)

                                             aCGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTCCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAAGTTGCGATTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGGCATAATAATATTGAATACacaata                                                                                                          >  2:744257/1‑145 (MQ=11)

                                                ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGAGACAGTTACATATTGAATACaaaatagaa                                                                                                       >  2:597769/1‑143 (MQ=11)

                                                  ttATCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCtt                                                                                                     <  2:155167/150‑1 (MQ=255)

                                                   tATCCGGCATGATAACCCTTCTAATCCAATAGATTGCAGTGAACGTGTAGTCCTTATAAGCGTAGCTCATCAGTCAATGTTGCGTTTGTCTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCttt                                                                                                    <  2:1069237/150‑1 (MQ=18)

                                                   tATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCttt                                                                                                    >  1:1441786/1‑150 (MQ=255)

                                                     tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                  <  2:213556/150‑1 (MQ=255)

                                                      ccGGCATGAAAACCCTTCAAATACAATAGATTGAAGTGAAAGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTAGCGTGTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                 <  2:172505/150‑1 (MQ=21)

                                                        tgaaTGAAAACCCTTAAAATCCAATAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                               <  2:1983304/147‑1 (MQ=255)

                                                        ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCCTTTGTCATCAGTTAAAAATGGCGCTGTAAACGGAGTCATTTTAATATTGAATACAACGAAGGCTTTGTTCa                                                                                               >  1:1567141/1‑150 (MQ=17)

                                                            tGAAAACCCTTCAAATCCTATAGATTGCAGTGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                           <  2:502731/150‑1 (MQ=255)

                                                              aaaaCCCTTCAAATCCAATAGTTTGCATTGAACGTGTAGGCCTTATATGCGTAGCGCATCATTCAATGTTGCGTTTGTAATCAGGTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                         <  1:1010747/150‑1 (MQ=25)

                                                              aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTATGCCTGATAAGCGTAGCGCATGAGGCAATGTTTCGTTTGTAATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                         <  2:1904871/150‑1 (MQ=21)

                                                              aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATACAAAGTAGGCTTTATTCAttcacg                                                                                         >  1:1063303/1‑145 (MQ=17)

                                                               aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGTCATTTTAATATTGTATAAAAAGAAGGCTATTTTCATTCCGGa                                                                                        >  2:1602674/1‑150 (MQ=11)

                                                                     ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGGCATTTTAATATTGTAGAAAAAGTATGCTTTTTTCAAGCAGGAttatgc                                                                                  >  2:1751113/1‑145 (MQ=14)

                                                                       cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCAGTATgaggcga                                                                                >  2:1763785/1‑149 (MQ=18)

                                                                        aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAATAGGCGGCTTTTTATTATTGTATAAAACGAAGGCTTTTTTCATTCAGGATCCGGGGtt                                                                               >  1:977242/1‑149 (MQ=17)

                                                                             cattAGATTGCAGTGAACTTGTAGGACTTATAATCGTAGAGAATCAGGCAATGTTGCGTTTGTAATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                          <  2:756/147‑1 (MQ=21)

                                                                             cAATAGATTGCAGCGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTGGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                          <  1:762864/150‑1 (MQ=255)

                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                         <  2:529299/150‑1 (MQ=255)

                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGAAGGCTTTGTTCATGCCGGATGCGGCGTGAACgac                                                                         >  2:1909364/1‑148 (MQ=255)

                                                                                tAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                       <  2:1149313/150‑1 (MQ=255)

                                                                                tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                       <  2:532529/150‑1 (MQ=255)

                                                                                    ttGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                   <  1:420993/150‑1 (MQ=255)

                                                                                         gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGTATGCGGCGTGAACGCCTTATCCGGCAt                                                              >  2:1943808/1‑150 (MQ=255)

                                                                                         gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGATAACGTAGGCTTTTTTCATGCCGGATACGGAGTGAACGCCTTATCCGGCAt                                                              >  2:1179063/1‑150 (MQ=25)

                                                                                         gTGAACGAGTCGGCCTCATAACCGGAGCGCATCAGGCAATCTGGCTGTTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                              <  1:118259/150‑1 (MQ=17)

                                                                                          tgtaagggtAGGCCTGATAATCGTAGCGCATCATGCAATGTTGCGTTTGTCATGAGTTTCAAAAGGCGGTTTAAAATGCGTGATTTTGATATTGTTGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                             <  2:1441786/143‑1 (MQ=21)

                                                                                          tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                             <  2:1199825/150‑1 (MQ=255)

                                                                                          tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                             <  1:1475359/150‑1 (MQ=255)

                                                                                             aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAATGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaa                                                          <  1:1158952/150‑1 (MQ=255)

                                                                                              cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                         >  1:874798/1‑150 (MQ=255)

                                                                                                    ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                   <  2:1593790/150‑1 (MQ=255)

                                                                                                    ggCCTCATAAGCGTACCGCATCAGGCAATGTTTCGTTTGCCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                   <  1:127277/150‑1 (MQ=255)

                                                                                                    ggCCAGATAAGCGTAGGTCATCAGGCAATGTTGCGTTACTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                   <  1:414470/150‑1 (MQ=18)

                                                                                                     tccTGATAATCGTAGCGAATCAGGCAATGTATCGTTTGTCATCAGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                  <  1:458701/149‑1 (MQ=18)

                                                                                                     gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                  >  1:1653869/1‑150 (MQ=255)

                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                 >  1:946423/1‑150 (MQ=255)

                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAAGAGTCATTTTAATATTGAATAAAAAGTAGGCTTTGTACATGCCGGATGCGGAGTGAACGGCTTATACTGCATGAAAACACTTaaa                                                 >  2:1338582/1‑150 (MQ=18)

                                                                                                       cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                                <  1:1477439/150‑1 (MQ=255)

                                                                                                        ttaTACGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                               <  2:238558/148‑1 (MQ=255)

                                                                                                          aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                             <  2:1290569/150‑1 (MQ=255)

                                                                                                            aaGTGTAGCTCATCACGGAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                           <  2:755817/150‑1 (MQ=255)

                                                                                                             ggagTATCGCATAAGGCATTGTTGCGTTTGTCATAAGTTTCAAGTGGCGCTGTAAATCGCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                          <  2:1839356/147‑1 (MQ=17)

                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                          >  2:1694335/1‑150 (MQ=255)

                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTAATATTGAAGACAACGAAGGCTTTGTTCATGCCGGATGCGTCGTGATAGACAATTACGGCATCAACACGCTTAAAATCCAATa                                         >  2:1056122/1‑150 (MQ=18)

                                                                                                              gCGTAGCACATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTATAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGTCGTGAAAGACTAATACGGCATGAAAACCCTTAAAATCCAATa                                         >  1:1092213/1‑150 (MQ=18)

                                                                                                               cGTAGCGCATCAGTCAATGTTGCGTTTGTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCATGAAAACCCTTCAAATCCAATAg                                        <  2:853206/150‑1 (MQ=255)

                                                                                                               cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTTTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCCGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTAAAATCCAATAg                                        >  1:998996/1‑150 (MQ=18)

                                                                                                                  aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCATGAAAACCCTTCAAATCCAATATAtt                                     >  1:797479/1‑150 (MQ=255)

                                                                                                                    cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTCATATTGTAGAAAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTAAAAACCAATATAATGc                                   >  1:1056306/1‑150 (MQ=17)

                                                                                                                        tCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGtt                               >  1:1762591/1‑149 (MQ=255)

                                                                                                                          aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTATACAACGTAGGCTTTGTTCATTCAGAATGAGTCGTTAACTCCTTATCCGACATGAAATCACTTCAAATCCAATAGATTGAAGTGaa                             >  1:152635/1‑150 (MQ=17)

                                                                                                                            tcAATGTTGCTTTAGTCTTCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                           <  2:501137/149‑1 (MQ=255)

                                                                                                                                 gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                      >  1:1542328/1‑150 (MQ=255)

                                                                                                                                 gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGACTTATACGGCATGAAAACCCTTCAAATACAATATATTGCAGTGAACGTGTAg                      >  2:255789/1‑150 (MQ=17)

                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGCTTGCAGTGAACGTGTAGGCCt                  >  2:1210157/1‑150 (MQ=255)

                                                                                                                                      gTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                 <  1:449672/150‑1 (MQ=255)

                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCCGTGAACGTGTAGGCCTGATAAGc           >  2:314682/1‑150 (MQ=255)

                                                                                                                                            catcaGTTTCAAATGGCGCTATAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  2:1770139/1‑150 (MQ=255)

                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:627277/1‑150 (MQ=39)

                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          <  2:643407/150‑1 (MQ=39)

                                                                                                                                              tcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCTTGATAAGCGt         <  2:1354887/150‑1 (MQ=255)

                                                                                                                                               caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa        >  2:140860/1‑150 (MQ=38)

                                                                                                                                                aGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg       <  2:1600173/150‑1 (MQ=18)

                                                                                                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      >  1:1116685/1‑150 (MQ=35)

                                                                                                                                                   ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc    >  2:1618961/1‑150 (MQ=34)

                                                                                                                                                     cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGGCTTTTCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:1748752/1‑150 (MQ=12)

                                                                                                                                                     cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:1898849/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: `‑`