BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,042	(G)_6→7	pseudogene (156/351 nt)	yrhC →	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,036	1	.	G	93.5%	369.6 / 18.0	108	pseudogene (150/351 nt)	yrhC	pseudogene fragment
Reads supporting (aligned to +/- strand): ref base . (3/4); new base G (44/57); total (47/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 112 total aligned reads displayed.

AAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATG  >  NC_000913/3623892‑3624180

aaaTACATACCAGACAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATATGCCAGAGTATTACAGTTAAGA‑gggggg                                                                                                                                            >  1:958838/1‑151 (MQ=25)

  atacatacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑gggggggc                                                                                                                                          >  2:1424044/1‑149 (MQ=25)

    acatacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑gggggggca                                                                                                                                         <  2:1208296/150‑4 (MQ=25)

      atacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑gggggggca                                                                                                                                         >  2:1298950/1‑145 (MQ=25)

       tacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGAc                                                                                                                                      >  2:279181/1‑151 (MQ=25)

        acCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGAc                                                                                                                                      <  2:1016889/150‑1 (MQ=25)

         ccAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGc                                                                                                                                    >  2:368035/1‑151 (MQ=25)

          cAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGAc                                                                                                                                      <  1:851004/148‑1 (MQ=25)

             aaaagaaaGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGtt                                                                                                                                >  2:1118046/1‑151 (MQ=25)

              aaagaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGtt                                                                                                                                >  1:252514/1‑150 (MQ=25)

                agaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGa                                                                                                                             <  1:1297833/151‑1 (MQ=25)

                 gaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGtt                                                                                                                                >  2:358326/1‑147 (MQ=25)

                   aagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTg                                                                                                                          <  2:1101603/151‑1 (MQ=25)

                   aagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTg                                                                                                                          >  1:1243891/1‑151 (MQ=25)

                        aaGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAAt                                                                                                                     <  1:359641/151‑1 (MQ=25)

                        aaGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAAtt                                                                                                                    >  2:365447/1‑151 (MQ=35)

                           tAAAGTGTATGGTGTATATAATAAATGTAAGAAATTAAATGCTGATGAATTATTTCGCTTTCTTGATGATCCCAATTCCGTGAAAATGTTTTCATCTGCCAGAGTAGTTCAGTTAATT‑GGGGTGGAAGTCGCGGTTAGATTGGCAATTGa                                                                                                                  <  1:644757/150‑1 (MQ=255)

                           tAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTg                                                                                                                   >  2:22462/1‑150 (MQ=25)

                                agTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGATAGATTGGCAATTGAAttat                                                                                                             >  1:958133/1‑149 (MQ=255)

                                  tATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGc                                                                                                           >  1:212130/1‑151 (MQ=255)

                                   aTGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGc                                                                                                           <  1:548271/150‑1 (MQ=255)

                                    tGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCAc                                                                                                         >  1:617168/1‑151 (MQ=255)

                                     ggTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACt                                                                                                        >  1:1142133/1‑151 (MQ=255)

                                       tttATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTg                                                                                                       >  2:815455/1‑150 (MQ=255)

                                        ttATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGAc                                                                                                     >  2:273040/1‑150 (MQ=255)

                                        ttATATAATAAATGTAAGAAACTAAATGATGATGAATTATTACGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGa                                                                                                      <  1:1260010/150‑1 (MQ=255)

                                              aataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTACGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAAt                                                                                               >  1:1212228/1‑151 (MQ=255)

                                              aataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAAt                                                                                               <  1:550279/151‑1 (MQ=255)

                                              aataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAAt                                                                                               >  1:346478/1‑151 (MQ=255)

                                               ataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAAtt                                                                                              <  1:342963/151‑1 (MQ=255)

                                                 aaaTGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATtat                                                                                            <  1:302363/151‑1 (MQ=255)

                                                  aaTGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATtata                                                                                           <  2:1167010/151‑1 (MQ=255)

                                                     gTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATcc                                                                                        >  1:858536/1‑151 (MQ=255)

                                                       aaGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGt                                                                                      >  1:1166141/1‑150 (MQ=40)

                                                        aGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGt                                                                                      >  1:1353323/1‑150 (MQ=255)

                                                        aGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGt                                                                                      >  2:831384/1‑150 (MQ=255)

                                                         gAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGt                                                                                      >  2:1023352/1‑149 (MQ=255)

                                                            aCTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTag                                                                                    >  2:845832/1‑148 (MQ=255)

                                                               aaaTGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGat                                                                                <  2:212140/149‑1 (MQ=255)

                                                               aaaTGATGATGAATTATTTCGCTTACTTCATGATCTCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACACTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAAat                                                                                <  1:436697/149‑1 (MQ=14)

                                                                aaTGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATc                                                                             <  1:422121/151‑1 (MQ=255)

                                                                 atgatgatgAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCg                                                                            <  2:814092/151‑1 (MQ=255)

                                                                   gatgatgaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCgg                                                                           <  1:16226/150‑1 (MQ=255)

                                                                      gatgaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGc                                                                        <  2:1197921/150‑1 (MQ=255)

                                                                       atgaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCAt                                                                      <  1:1034653/151‑1 (MQ=255)

                                                                         gaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCAttt                                                                    >  2:173800/1‑151 (MQ=255)

                                                                         gaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCAtt                                                                     <  1:1334988/150‑1 (MQ=255)

                                                                           attattTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTtat                                                                  >  2:129699/1‑151 (MQ=255)

                                                                             tattTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGCTTGGCAATTGAATTCTGCACTGACAAAAATCATATCCGTAGAGATATCGGAGCATTTATAc                                                                >  2:780627/1‑151 (MQ=255)

                                                                             tattTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATGATATCCGTAGAGATATCGGAGCATTTATAc                                                                >  2:792981/1‑151 (MQ=255)

                                                                               tttCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTc                                                              >  1:535037/1‑151 (MQ=255)

                                                                                 tCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCgg                                                            >  1:993715/1‑151 (MQ=255)

                                                                                  cGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCgg                                                            <  1:69903/150‑1 (MQ=255)

                                                                                   gCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCggg                                                           <  1:926674/150‑1 (MQ=255)

                                                                                      tACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCa                                                         <  2:530800/149‑1 (MQ=255)

                                                                                      tACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCa                                                         <  1:632586/149‑1 (MQ=255)

                                                                                      tACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCa                                                         <  2:1149125/149‑1 (MQ=255)

                                                                                      tACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGTCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCaaa                                                       <  1:1393809/150‑1 (MQ=255)

                                                                                          tgatgaTCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAatat                                                    <  2:861083/150‑1 (MQ=255)

                                                                                          tgatgaTCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAatat                                                    <  2:360462/150‑1 (MQ=255)

                                                                                           gatgatCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATaa                                                  <  1:610045/151‑1 (MQ=255)

                                                                                           gatgatCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATaaa                                                 <  2:880813/151‑1 (MQ=255)

                                                                                            atgatCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAatat                                                    <  1:224448/148‑1 (MQ=255)

                                                                                               atCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAAttt                                              <  1:104528/151‑1 (MQ=255)

                                                                                               atCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAAttt                                              <  1:132997/151‑1 (MQ=255)

                                                                                               atCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTg                                             <  1:579896/151‑1 (MQ=255)

                                                                                                 cGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGc                                            <  2:1117328/151‑1 (MQ=255)

                                                                                                 cGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGc                                            >  2:1047380/1‑151 (MQ=255)

                                                                                                    aaTTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaa                                         <  1:902973/151‑1 (MQ=255)

                                                                                                    aaTTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaa                                         <  2:26212/151‑1 (MQ=255)

                                                                                                     aTTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaaa                                        <  2:959753/151‑1 (MQ=255)

                                                                                                     aTTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaaa                                        <  2:827285/151‑1 (MQ=255)

                                                                                                      ttCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaaa                                        >  2:427749/1‑150 (MQ=255)

                                                                                                        ccTTGAAAAGGAGTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAAGTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCaaaaaa                                      <  1:755956/150‑1 (MQ=255)

                                                                                                          ttGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAATTTTGCAAAAAATGc                                   >  1:84912/1‑151 (MQ=255)

                                                                                                            gAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCg                                  <  1:1059899/150‑1 (MQ=255)

                                                                                                            gAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCg                                  <  2:1258090/150‑1 (MQ=255)

                                                                                                                   aTTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAAt                          <  2:1206269/151‑1 (MQ=255)

                                                                                                                     ttCATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATg                         <  1:183013/150‑1 (MQ=255)

                                                                                                                      tCATCTGCCAGAGGACTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGCTATCGGAGCATTTATACTCGGGCCAATATAAATTTGCAAAAAATGCGCAGCTAATgct                       >  1:1343315/1‑149 (MQ=14)

                                                                                                                       cATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGttt                      <  2:797092/151‑1 (MQ=255)

                                                                                                                       cATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGttt                      >  2:50219/1‑151 (MQ=255)

                                                                                                                       cATCTGCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGtt                       <  2:541147/150‑1 (MQ=255)

                                                                                                                            gCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTaa                  <  1:839443/150‑1 (MQ=255)

                                                                                                                            gCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTaa                  <  2:886961/150‑1 (MQ=255)

                                                                                                                            gCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTa                   <  1:681551/149‑1 (MQ=255)

                                                                                                                            gCCAGAGTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAat                 <  2:887142/151‑1 (MQ=255)

                                                                                                                               agagTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGAC‑AAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATAttt             >  2:933435/1‑151 (MQ=255)

                                                                                                                                gagTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATAttt             <  1:205905/151‑1 (MQ=255)

                                                                                                                                gagTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATAttt             <  1:144850/151‑1 (MQ=255)

                                                                                                                                 agTATTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATAtttt            <  1:231544/151‑1 (MQ=255)

                                                                                                                                    aTTACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGaa         <  2:1212859/151‑1 (MQ=255)

                                                                                                                                     ttACAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGaa         >  2:95358/1‑150 (MQ=255)

                                                                                                                                       aCAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACaa      <  1:1304465/151‑1 (MQ=255)

                                                                                                                                       aCAGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACaa      <  2:1335487/151‑1 (MQ=255)

                                                                                                                                         aGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAat     >  2:204320/1‑150 (MQ=255)

                                                                                                                                         aGTTAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAat     >  1:1059265/1‑150 (MQ=255)

                                                                                                                                           ttAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATg  >  2:611567/1‑151 (MQ=255)

                                                                                                                                           ttAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATg  >  1:1370492/1‑151 (MQ=255)

                                                                                                                                           ttAAGAGGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATg  >  2:732410/1‑151 (MQ=255)

AAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGA‑GGGGGGCAAGACGCTGTTAGATTGGCAATTGAATTCTGCACTGACAAAAATTATATCCGTAGAGATATCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATG  >  NC_000913/3623892‑3624180

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`