BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,092,774:1	+A	7.8%	coding (897/951 nt)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,092,774	1	.	A	7.8%	478.4 / 34.7	129	coding (897/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand): ref base . (73/46); new base A (4/6); total (77/52)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 132 total aligned reads displayed.

AGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092625‑3092920

aGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtg                                                                                                                                                     <  2:12653/150‑1 (MQ=255)

  tGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGATAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑Tc                                                                                                                                                  >  1:1392458/1‑150 (MQ=255)

  tGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑Tc                                                                                                                                                  >  1:677104/1‑150 (MQ=255)

   tccTGGATAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAATAGGGATGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTTTGCTTGAGATTGAATCAGAGTTTCCGGTG‑TCg                                                                                                                                                 <  2:1506500/148‑1 (MQ=255)

   gACTGGAAAATCGCCCGTCAGATCGGGGCGACGCTGATAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCATCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCg                                                                                                                                                 <  2:673524/150‑1 (MQ=255)

     cTGGAACATCGCCCGGCAGATCGGGCCGACGCTGAAAGAAACAGGGCTGATTTTTGTTGTTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtgatcga                                                                                                                                                <  1:804275/150‑6 (MQ=255)

     cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                               >  1:1120575/1‑150 (MQ=255)

     cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                               >  2:911614/1‑150 (MQ=255)

     cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                               >  1:1661476/1‑150 (MQ=255)

     cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                               >  2:465307/1‑150 (MQ=255)

     cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCAACGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                               >  2:1518053/1‑150 (MQ=255)

      tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtgatcgac                                                                                                                                               <  2:1612117/150‑7 (MQ=255)

        gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAATTTTGTTGGTCTGGATATCATCGGCGACCGACTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGAATGAAGCAGAGTATCCGGAG‑TCGATCAc                                                                                                                                            >  2:627839/1‑150 (MQ=255)

            aTCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGa                                                                                                                                        >  2:114975/1‑150 (MQ=255)

               gCCCGTCAGATCGGGCCGACGCTTAAAGAAAAAGTGCTGATTTTTGTTTGTATGGATATAATCTGCGAACGTCTGACTTAAATTAAATTCACCAGCCCAACCTGTATTCGTTAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                     <  1:1321763/150‑3 (MQ=255)

               gCCCGTCAGATCGGGCCGACGCTGAAAGAACAAGGCCTGACTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACTTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                     <  1:1020256/150‑3 (MQ=255)

               gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAAATAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATg                                                                                                                                     >  2:790498/1‑150 (MQ=255)

               gCCCGGCATATCGGGCAGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                     <  1:333898/150‑3 (MQ=255)

                  cGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑GCGATCACCGGAATGTTa                                                                                                                                  <  1:1574377/150‑1 (MQ=255)

                     cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATg                                                                                                                               <  2:53104/150‑1 (MQ=255)

                     cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATg                                                                                                                               >  2:977970/1‑150 (MQ=255)

                       gATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGAAATCATCGGCGACCGTCTGACTGAAATTAAGGTCACCCGCGCAACCTGTATTCGTGTGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAAAGTTAATGGa                                                                                                                             >  2:319715/1‑150 (MQ=255)

                         tCGGTCCGTCGCTGAAATAAAAAGGGCTGATGTATGTTGGTCTGGATTTCATCGGCGACCGTCTGACTGAAATTCACGGCCCCACCCCAACCTGTATACGTGAGAATCAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATg                                                                                                                           <  2:345369/150‑1 (MQ=255)

                         tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATg                                                                                                                           >  2:1705633/1‑150 (MQ=255)

                            ggCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCa                                                                                                                        >  1:134190/1‑150 (MQ=255)

                                    cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGTATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGTATGCCATCTAAGCa                                                                                                                >  1:1634661/1‑150 (MQ=255)

                                    cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTCACTGAAATTAACGTCACCAGCCCAACCTGAATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                >  2:968537/1‑150 (MQ=255)

                                      gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCAGAATGTTAATGGATGCCATCGAAGCACg                                                                                                              >  2:1122037/1‑150 (MQ=255)

                                        aagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                            >  2:5622/1‑150 (MQ=255)

                                        aagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGAATGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                            >  2:461198/1‑150 (MQ=255)

                                         agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGttt                                                                                                           >  1:1708692/1‑150 (MQ=255)

                                         agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGttt                                                                                                           <  1:454421/150‑1 (MQ=255)

                                          gaaaGAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATGGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTACTGGATGCCATCGAAGCACGTTTa                                                                                                          >  2:1287243/1‑150 (MQ=255)

                                          gaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTa                                                                                                          >  2:73861/1‑150 (MQ=255)

                                           aaaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCTCCTTTAc                                                                                                         >  2:1056272/1‑150 (MQ=255)

                                             aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                       <  2:1226624/150‑1 (MQ=255)

                                                gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTTTTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                    <  1:871905/150‑1 (MQ=255)

                                                gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                    >  1:1042963/1‑150 (MQ=255)

                                                   cTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                 >  1:1614512/1‑150 (MQ=255)

                                                         tttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAccc                                                                                           <  1:109336/150‑1 (MQ=255)

                                                              tgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCtt                                                                                      >  2:1159987/1‑150 (MQ=255)

                                                               ggTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATTGATGCCATCGAAGAACGTTTACAGCAGCAGTAACCCACCTTa                                                                                     >  2:910470/1‑150 (MQ=255)

                                                               ggTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                     >  2:75964/1‑150 (MQ=255)

                                                                gTCTGTATTTCATAGGCGATCGTCTGACTGAAATTAACGTCAACAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAg                                                                                    <  1:1480590/150‑1 (MQ=255)

                                                                    ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGAAGCAGTAACCCACCTTAGCgag                                                                                >  2:337683/1‑150 (MQ=255)

                                                                    ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGACCGACCTGTATTAGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTatacat                                                                                >  1:1474235/1‑145 (MQ=255)

                                                                        atcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAgg                                                                            <  1:1565014/150‑1 (MQ=255)

                                                                          catcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTAGTGAGAGTGAAGAAGAGTTTCCGTAG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTATACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                          >  2:1113332/1‑150 (MQ=255)

                                                                           atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtc                                                                         <  2:1662936/150‑1 (MQ=255)

                                                                             cGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                       >  2:396698/1‑150 (MQ=255)

                                                                             cGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                       <  2:27628/150‑1 (MQ=255)

                                                                              ggCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCg                                                                      <  1:1122037/150‑1 (MQ=255)

                                                                                cGACCGTCTGACTGAAATTAACGTCACCAG‑ACAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                   <  2:500295/150‑1 (MQ=255)

                                                                                 taCCGTCTTACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCATCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                   <  2:444019/149‑1 (MQ=255)

                                                                                  aCCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                  >  1:681054/1‑150 (MQ=255)

                                                                                     gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                               <  2:466191/150‑1 (MQ=255)

                                                                                     gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                               <  1:1433026/150‑1 (MQ=255)

                                                                                     gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                               <  1:5622/150‑1 (MQ=255)

                                                                                       ctgactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGActc                                                             >  2:1379073/1‑150 (MQ=255)

                                                                                         gactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTTATGGATGCCATGGAAGCACGTTTACAGCAGCAGTAACCCACCGTAGCGAGAAGGTTCTCGTTGAGACGCTg                                                           >  2:285435/1‑150 (MQ=255)

                                                                                          actgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                          <  1:1513487/150‑1 (MQ=255)

                                                                                           ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                         <  1:525677/150‑1 (MQ=255)

                                                                                            tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                        <  1:461198/150‑1 (MQ=255)

                                                                                              aaaTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGa                                                      <  1:1123089/150‑1 (MQ=255)

                                                                                                aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCAACGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACa                                                    >  2:286929/1‑150 (MQ=255)

                                                                                                  tAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAgc                                                  >  2:1626751/1‑150 (MQ=255)

                                                                                                    gcGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACCGCAGCAGTAACCCACCTTAGCGAGAAGGATGACGTTGAGACTCTGAGTGACAgcgc                                                >  1:1065722/2‑150 (MQ=255)

                                                                                                      gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGccc                                              >  1:1218132/1‑150 (MQ=255)

                                                                                                       tCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCct                                             <  1:306667/150‑1 (MQ=255)

                                                                                                          ccAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttct                                          >  1:737227/1‑150 (MQ=255)

                                                                                                           cAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttctt                                         <  1:1705633/150‑1 (MQ=255)

                                                                                                               ccAACCTGTATTTGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCa                                     <  2:1538874/150‑1 (MQ=255)

                                                                                                                    cTGTTTTCGTGATATTGAAGGAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGT‑ACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATa                               <  1:1622499/150‑1 (MQ=255)

                                                                                                                      gTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGAGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGAATAc                              >  2:1523026/1‑150 (MQ=255)

                                                                                                                       tATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACt                             <  1:1255156/150‑1 (MQ=255)

                                                                                                                        aTTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCCCCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTg                            >  1:1287640/1‑150 (MQ=255)

                                                                                                                           cGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                         >  1:483347/1‑150 (MQ=255)

                                                                                                                           cGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                         >  2:1495587/1‑150 (MQ=255)

                                                                                                                            gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcg                        >  1:14491/1‑150 (MQ=255)

                                                                                                                             tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcgc                       <  1:566761/150‑1 (MQ=255)

                                                                                                                              gagaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCCCt                      >  2:179386/1‑150 (MQ=255)

                                                                                                                                    gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCt                      <  1:1539841/144‑1 (MQ=255)

                                                                                                                                    gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                <  1:790498/150‑1 (MQ=255)

                                                                                                                                    gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                >  2:1573222/1‑150 (MQ=255)

                                                                                                                                    gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                >  2:1304301/1‑150 (MQ=255)

                                                                                                                                      aGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCGACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttt              >  1:880518/1‑150 (MQ=255)

                                                                                                                                       gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttt             >  1:180617/1‑150 (MQ=255)

                                                                                                                                        cAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            >  2:914362/1‑150 (MQ=255)

                                                                                                                                         agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTATAAAGCACCAGTAACCCACATTAGCGAGAAGGATATCGTTGAGACTCTGAGTGACCGCTCCCTTCTATACACGCATACTGGGCGCAGTTGCttattt           >  1:339597/1‑150 (MQ=255)

                                                                                                                                          gagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg          <  2:681054/150‑1 (MQ=255)

                                                                                                                                          gagTTTCCGGCG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg          <  2:1106146/150‑1 (MQ=255)

                                                                                                                                          gagTCTCCGGTG‑TCGATCACCGGAATGGTAATGGTTGCCCTCGCAGGACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCGCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACTCATACTGGGCGCTGTTGCTTTTTTg          <  1:319715/150‑1 (MQ=255)

                                                                                                                                          cagTTTCCGGTG‑GCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAGCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg          <  2:472935/149‑1 (MQ=255)

                                                                                                                                           agTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGa         >  1:603290/1‑150 (MQ=255)

                                                                                                                                           acattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTGTCCCCGCATACTGGGCGCTGGTGCTTTTTTGa         >  1:347502/4‑150 (MQ=255)

                                                                                                                                            gTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  1:1718140/1‑150 (MQ=255)

                                                                                                                                            cattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  2:979165/3‑150 (MQ=255)

                                                                                                                                             attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc       >  1:768515/2‑150 (MQ=255)

                                                                                                                                              ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAcc      >  2:363400/1‑150 (MQ=255)

                                                                                                                                               tCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa     >  1:306061/1‑150 (MQ=255)

AGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092625‑3092920

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: `‑`