BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,224,239	A→T	6.3%	intergenic (+332/+40)	ycgI → / ← minE	pseudogene/cell division topological specificity factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,224,239	0	A	T	6.3%	943.6 / 37.6	267	intergenic (+332/+40)	ycgI/minE	pseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand): ref base A (125/125); new base T (8/9); total (133/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

Only 100 of 274 total aligned reads displayed.

TGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGACGGCGTTCAGCAACAATAATCTGCAGCCGTTCTTTTGCAATGTTGGCTGTGTTT  >  NC_000913/1224017‑1224510

  ggAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                                                                                                                                                            <  2:637828/275‑1 (MQ=255)

             aataatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGANATAGCAATGAGGAGTATCAGCAAGAATACTAGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAAANCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAANGNGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGt                                                                                                                                                                                                            >  3:83540/1‑279 (MQ=255)

                                                       cTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTAAGCTCAAGAATAGAAATATCGCCATCttttt                                                                                                                                                                      >  2:84257/1‑275 (MQ=255)

                                                                aGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTAGCCGCTTTACGACAACGTGGATGAGAGGGATGAAAGACTCAAGGCAGAGAGAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCGTATTTCCGCTCTTCTGCT‑TCCGGTAAGGTACGGTTCCGCTCAAGAATAGAAATAACGCATTCTTTTTggccaaag                                                                                                                                                              >  4:670705/1‑267 (MQ=255)

                                                                   ttAATTAGTCAACCACGGATAATACTTTTGTTGTGTAACTGGGAAAATCCGTAATGTTTAAATGGCGAGTTAGAAAGAGCAATGAGTAGTATTAGCAAGAATACTCGCCCCTTTACCATAACGTGGATGAGAGGTATGAAAAACGCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCCGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTa                                                                                                                                                          <  4:442047/275‑1 (MQ=255)

                                                                      tttACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                                                                                                                                                                                                                           <  2:374554/176‑1 (MQ=255)

                                                                      tttACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCCGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGc                                                                                                                                                              >  2:442051/1‑268 (MQ=255)

                                                                      tttACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCACGAACAGAAATATCCCCATCTTTTTGCTCAAGCTGtcagt                                                                                                                                                       >  4:369049/1‑271 (MQ=255)

                                                                      tttACTCAACCAGGCATAAACATTTTGCCTTGTGCGTGGGAACAGCCTTAAGGCGTAAAGGTGGTGGTGGAAATAGCAATGAGGATTACAAGCAAGAATACTCGCCGCTTTACCACAACGCGGTTGAGAGGGATGAAAAACTCAAGGCCGAGATAACTCTGCCTTGAAGATAAATCCGATGTTACAGCGGGCTTACTTCAGCGCTTCTGTT‑TACGGTAAGGTCACGTTCAGCTCAAGAACAGCATCATCTCCCTCTTTTTGCTCAAGCTTTACg                                                                                                                                                        >  4:91665/1‑274 (MQ=255)

                                                                           tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCACTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                                                                                                                                                 >  6:896650/1‑149 (MQ=255)

                                                                           tCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                                                                                                                                                 <  8:10968/149‑1 (MQ=255)

                                                                            caacaaGGCATAAACATTTTCTTTTGAGCGTGGGATCAGCATTAAGGTGCAAAGGGGGAGGTGGAAATAGCAAAGAGGAGTATCAGCCAGAATACTCGCCGCTTTACCACACCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTTAAGATAAATGCGTTTCAACAGCGGGCTTATATCAGCTCTTCTGCT‑TCAGGTAAGGTCACGTTCCGCTCAAGAATAGAAATATCGCCCTTTTTTTGCTCAAGCTGTTCGGTTAcc                                                                                                                                                  >  4:42183/1‑274 (MQ=255)

                                                                             aaCCAGACATACACATTTTTTGTTTTGCGTGGGAACAGCCTTAAGGTGTATGGGGCTAGGTGGCAATAGCTTTGAGACGTATCAGCAAAAATACTTGCCCCTGTACCACCGCGGGGCTGTGAGGGTTGAAAAACCCAAGGCAGAGATAACTCTGCCTTGAAGACAATTGCGTTCTTCCAGCGCGTTTATGTCATCTCTT‑TGCTGAACGAAAAAGTTACGTTCACCTCCAGTATAGGAATACCTCCTTCTTTTTGCTTTAGCTTTAGGATTCCCAt                                                                                                                                                >  4:311912/1‑275 (MQ=255)

                                                                                cAGGCATAAACATTTTGTTTTGTGCGTTGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATGGAAATATCGCCATCTTTTTGCTCAAGCTg                                                                                                                                                            <  2:42183/260‑1 (MQ=255)

                                                                                   gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                                                                                                                                                         >  6:485335/1‑149 (MQ=255)

                                                                                   gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAgg                                                                                                                                          <  2:115531/275‑1 (MQ=255)

                                                                                     aTAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                                                                                                                                                                                                       >  8:1453037/1‑149 (MQ=255)

                                                                                       aaaCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                                                                                                                                                                                            >  1:48471/1‑174 (MQ=255)

                                                                                        aaCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                                                                                                                                                                                                    >  6:1582732/1‑149 (MQ=255)

                                                                                        aaCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                                                                                                                                                                                                    <  7:1392502/149‑1 (MQ=255)

                                                                                         aCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGag                                                                                                                                                                                                                                                                   <  5:1085477/149‑2 (MQ=255)

                                                                                         aCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                                                                                                                                                                                                   <  5:1440666/149‑1 (MQ=255)

                                                                                         aCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                                                                                                                                                                                                   <  8:270235/149‑1 (MQ=255)

                                                                                           aTTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                                                                                                                                                       >  5:1658885/1‑143 (MQ=255)

                                                                                           aTTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                                                                                                                                                       <  6:1658885/143‑1 (MQ=255)

                                                                                            ttttgttttgTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                                                                                                                                                                                                >  6:140939/1‑149 (MQ=255)

                                                                                             tttgttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                                                                                                                                                                                                                             <  1:130883/151‑1 (MQ=255)

                                                                                             tttgttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                                                                                                                                                                                                                             >  3:130883/1‑151 (MQ=255)

                                                                                                gttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                                                                                                                                                            <  5:1009747/149‑1 (MQ=255)

                                                                                                  tttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                                                                                                                                                                                                          <  5:140939/149‑1 (MQ=255)

                                                                                                    tgtgCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                                                                                                                                                                                                        <  6:310706/149‑1 (MQ=255)

                                                                                                     gtgCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                                                                                                                                                                                                                                       <  5:773285/149‑1 (MQ=255)

                                                                                                         gTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                                                                                                                                                                                                                                          <  5:404079/126‑1 (MQ=255)

                                                                                                               aCAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                                                                                                                                                                                                      <  6:1121750/140‑1 (MQ=255)

                                                                                                                    cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                                                                                                                                                                                                                           <  2:514091/130‑1 (MQ=255)

                                                                                                                           tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                                                                                                                                                                                                 <  6:379827/149‑1 (MQ=255)

                                                                                                                             tAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAACTATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCCCAGGCTCAATTTGTACGTTTTTAAAAATTACCTCAAGAAAATCctt                                                                                                >  4:618679/1‑275 (MQ=255)

                                                                                                                                aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                                                                                                                                                                                              <  7:657870/99‑1 (MQ=255)

                                                                                                                                aGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                                                                                                                                                                                              >  8:657870/1‑99 (MQ=255)

                                                                                                                                 gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                                                                                                                                                                                                                           <  5:485335/149‑1 (MQ=255)

                                                                                                                                  ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                                                                                                                                                  <  8:788614/141‑1 (MQ=255)

                                                                                                                                  ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                                                                                                                                                  >  7:788614/1‑141 (MQ=255)

                                                                                                                                   gggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGc                                                                                                                                                                                                                         <  8:639877/149‑1 (MQ=255)

                                                                                                                                    ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCt                                                                                                                                                                                                                        >  5:699455/1‑149 (MQ=255)

                                                                                                                                     gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAAt                                                                                                                  >  4:175770/1‑249 (MQ=255)

                                                                                                                                      aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑Tc                                                                                                                                                                                                                     >  5:30598/1‑149 (MQ=255)

                                                                                                                                       ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑Tcc                                                                                                                                                                                                                    <  6:696890/149‑1 (MQ=255)

                                                                                                                                        gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                                                                                                                                                                                                                             <  6:774806/108‑1 (MQ=255)

                                                                                                                                        gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                                                                                                                                                   <  6:745990/134‑1 (MQ=255)

                                                                                                                                        gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                                                                                                                                                   >  5:745990/1‑134 (MQ=255)

                                                                                                                                        gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCg                                                                                                                                                                                                                   >  7:112475/1‑149 (MQ=255)

                                                                                                                                         tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTg                                                                                                                                 <  4:250860/230‑1 (MQ=255)

                                                                                                                                          ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                                                                                                                                                                                                          <  5:593992/109‑1 (MQ=255)

                                                                                                                                           gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                                                                                                                                                            <  8:1204147/138‑1 (MQ=255)

                                                                                                                                            aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTaa                                                                                                                                                                                                               <  7:1316716/149‑1 (MQ=255)

                                                                                                                                            aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTaa                                                                                                                                                                                                               <  8:1607086/149‑1 (MQ=255)

                                                                                                                                            aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTaa                                                                                                                                                                                                               <  7:91425/149‑1 (MQ=255)

                                                                                                                                               tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGt                                                                                                                                                                                                            >  7:1556342/1‑149 (MQ=255)

                                                                                                                                                  cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCAc                                                                                                                                                                                                         >  6:1715378/1‑149 (MQ=255)

                                                                                                                                                      gaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAat                                                                                                      <  2:202370/244‑1 (MQ=255)

                                                                                                                                                             aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCaa                                                                                                                                                                                              <  5:986307/149‑1 (MQ=255)

                                                                                                                                                                aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                                                                                                                                                 >  5:1359357/1‑64 (MQ=255)

                                                                                                                                                                aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                                                                                                                                                                                                 <  6:1359357/64‑1 (MQ=255)

                                                                                                                                                                aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGaa                                                                                                                                                                                           >  5:214575/1‑149 (MQ=255)

                                                                                                                                                                 gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                                                                                                                                                                                                        <  8:812522/88‑1 (MQ=255)

                                                                                                                                                                 gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAAt                                                                                                                                                                                          >  7:407316/1‑149 (MQ=255)

                                                                                                                                                                 gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCa                                                                                                           >  2:677308/1‑228 (MQ=255)

                                                                                                                                                                   aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAg                                                                                                                                                                                        >  8:341882/1‑149 (MQ=255)

                                                                                                                                                                      aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTAATTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGa                                                        <  2:284881/274‑1 (MQ=255)

                                                                                                                                                                        tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAg                                                                                                                                                                                                              <  6:582651/122‑1 (MQ=255)

                                                                                                                                                                         aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                                                                                                                                                   >  7:282073/1‑101 (MQ=255)

                                                                                                                                                                           tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                                                                                                                                                   >  5:1585610/1‑99 (MQ=255)

                                                                                                                                                                            cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCg                                                                                     >  1:167183/1‑239 (MQ=255)

                                                                                                                                                                             gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCAc                                                                                                                                                                                                         >  8:1293786/1‑122 (MQ=255)

                                                                                                                                                                             gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCAc                                                                                                                                                                                                         <  7:1293786/122‑1 (MQ=255)

                                                                                                                                                                                gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc                                                                                                                                                                           >  8:430462/1‑149 (MQ=255)

                                                                                                                                                                                gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGACGGCGTTCAGc                                             >  2:598278/1‑275 (MQ=255)

                                                                                                                                                                                  tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTa                                                                                                                                                                                                                <  8:1515010/110‑1 (MQ=255)

                                                                                                                                                                                  tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTa                                                                                                                                                                                                                >  7:1515010/1‑110 (MQ=255)

                                                                                                                                                                                   ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttt                                                                                                                                                                        <  7:877926/149‑1 (MQ=255)

                                                                                                                                                                                     aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttttt                                                                                                                                                                      <  6:27528/149‑1 (MQ=255)

                                                                                                                                                                                        acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa                                                                                                                                                                                      <  5:1176069/130‑1 (MQ=255)

                                                                                                                                                                                        acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCt                                                                                                                                                                   >  6:737916/1‑149 (MQ=255)

                                                                                                                                                                                              tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGACGGCGTTCAGCAACAATAATCTGCa                               >  2:543285/1‑275 (MQ=255)

                                                                                                                                                                                                gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCCGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTTCGACGGCGTTACACAACAATAATCTGc                                >  4:496282/1‑272 (MQ=255)

                                                                                                                                                                                                     gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAcc                                                                                                                                                  <  4:731036/153‑1 (MQ=255)

                                                                                                                                                                                                      agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGACGGCGTTCAGCAACAATAATCTGCAGCCGTTCttt                     >  1:211008/1‑277 (MQ=255)

                                                                                                                                                                                                       gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTa                                                                                                                                                    <  8:553406/149‑1 (MQ=255)

                                                                                                                                                                                                        ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAc                                                                                                                                                   <  5:1301869/149‑1 (MQ=255)

                                                                                                                                                                                                        ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAc                                                                                                                                                   >  6:1238766/1‑149 (MQ=255)

                                                                                                                                                                                                         gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCAc                                                                                                                                                                                                         <  6:1474549/94‑1 (MQ=255)

                                                                                                                                                                                                         gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGc                                                                                                                                                                                                  >  5:1153193/1‑101 (MQ=255)

                                                                                                                                                                                                         gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGc                                                                                                                                                                                                  <  6:1153193/101‑1 (MQ=255)

                                                                                                                                                                                                             aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtct                                                                                                                                              <  5:737916/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                                                                                                                                                                                                       >  6:1338306/1‑55 (MQ=255)

                                                                                                                                                                                                                   tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGAt                                                                                                                                        >  5:1523470/1‑149 (MQ=255)

                                                                                                                                                                                                                       ggCAGAGTTATCTCTGCCTTGAGGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTa                                                                                                                                                    >  7:487600/1‑133 (MQ=255)

                                                                                                                                                                                                                       ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAt                                                                                                                                    <  7:389316/149‑1 (MQ=255)

                                                                                                                                                                                                                            agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTa                                                                                                                               <  6:1431982/149‑1 (MQ=255)

                                                                                                                                                                                                                            agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTa                                                                                                                               <  6:1523470/149‑1 (MQ=255)

TGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCT‑TCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACGTATTTACAAATGACCTCAAGAATATCTTTACGCAACTGCGGCAGATAATGCGGTTCTGCATCGCTGCGACGGCGTTCAGCAACAATAATCTGCAGCCGTTCTTTTGCAATGTTGGCTGTGTTT  >  NC_000913/1224017‑1224510

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`