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%	380.6 / 31.0	102	coding (897/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand): ref base . (42/52); new base A (4/4); total (46/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

Only 100 of 107 total aligned reads displayed.

GTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAA  >  NC_000913/3092626‑3092922

gTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtg                                                                                                                                                       <  1:1167009/149‑1 (MQ=255)

  gACTGGAAAATCGCCCGTCAGAGCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGGTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑Tc                                                                                                                                                    <  2:193307/149‑1 (MQ=255)

     tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAGGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                 <  2:965847/149‑1 (MQ=255)

     tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGTCAGATTTTTGTTGGTCTGGATATCATCGGCGCCCGTCTGCCTGCACTTAACGTCCCCCTCCGAGCCGGTATTCGTGAGAGTGTAGCTGAGTTTCCGGTG‑TCGAt                                                                                                                                                 >  2:317426/1‑149 (MQ=255)

     tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTTGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                 <  2:320474/149‑1 (MQ=255)

     tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTAGTCGTGAGCTTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                 >  1:395721/1‑149 (MQ=255)

      ggAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGACTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATc                                                                                                                                                >  1:259666/1‑149 (MQ=255)

       gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGCCCGGCTGACTGATATTAACGTCACCAGCCCAACCTGGATGCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCa                                                                                                                                               >  1:1468795/1‑149 (MQ=255)

            tCGCCCGTCAGATCGGGCCGTCGCTGATATAAAATGGTCTGATTTTTGTTGGTATGGATATCATCGGCGAACGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGa                                                                                                                                          <  1:905131/149‑1 (MQ=255)

            ccGCCCGGCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTTCATTTTGGTGGGCAGCCTATCATCGGCGACCGTCTGACTGAAATTCACGTCACCAGCCCAACCTGTATTCGTGCGCTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGa                                                                                                                                          <  1:1026699/148‑1 (MQ=255)

             cGCCCGTCAGATCGTGCCGAAGGACAAAGAACAAGGGCTGTTTTTTGTTGGTCTGGATAGCATCGGCGACCGTCTGACTGAAATTCACGTCACCAGCCCAACCTGTTTTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaa                                                                                                                                         <  2:1053353/149‑1 (MQ=255)

              gCCCGTCAGAGCGGGCCGACGCTGAAAGTAACAGGGCTGATGTTTGTTGGTCTGGATCTCATCGGCGACCGTCTGACTGAAATTAACCTCACCCGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                        <  2:1046309/149‑2 (MQ=255)

               cccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATg                                                                                                                                       >  2:607266/1‑149 (MQ=255)

               cccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATg                                                                                                                                       <  2:266709/149‑1 (MQ=255)

                 cGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGATATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaactc                                                                                                                                     <  1:582275/149‑5 (MQ=255)

                 cGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAAtttt                                                                                                                                     >  2:1230685/1‑149 (MQ=255)

                  gTCAGATCGGGCCGACGCTGAAAGATAAAGTGCTGATTTTTGTTGGTCTGGATATGATGGGCGACCGTCTGACTGATATTTACGTCACCAGCACAACCTGTATACGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCAGCGGAATGtaa                                                                                                                                    <  2:1130642/149‑3 (MQ=255)

                  gTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaactct                                                                                                                                    <  1:1351913/149‑6 (MQ=255)

                   tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTaa                                                                                                                                   >  1:1263046/1‑149 (MQ=255)

                    cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                  >  2:271200/1‑149 (MQ=255)

                    cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACAGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                  <  1:597464/149‑1 (MQ=255)

                      caTCGGGCCGACGCTGAAAGAAAAAGGGCTGATATTTGTTGCTCTGGATATCATCGGCGACCGTCTGACTGAACTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGGTTCCGGTG‑TCGATCACCGGAATGTTAATgg                                                                                                                                <  1:65607/148‑1 (MQ=255)

                       aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCACCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGa                                                                                                                               >  2:1368179/1‑149 (MQ=255)

                       aTCGGGCCGACGCTGAAAGAAAAAGGGATGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTACGGTG‑TCGATCACCGGAATGTTAATGGa                                                                                                                               >  1:721212/1‑149 (MQ=255)

                        tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGAt                                                                                                                              <  1:607266/149‑1 (MQ=255)

                            gCCGAAGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATAAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCa                                                                                                                          <  2:1065757/149‑1 (MQ=255)

                             ccGACGCTGAAAGAAAAAGGGCGGATTTTTGTTGGTCTGGATATCCTCGGCGACCGTCTGACTGAACTTCACGTCACCAGCCCAACCTGTATTCGTGAGATTGACGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCAt                                                                                                                         <  2:61045/149‑1 (MQ=255)

                              cGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATc                                                                                                                        >  2:158745/1‑149 (MQ=255)

                               gACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCg                                                                                                                       >  2:70200/1‑149 (MQ=255)

                               gACGCTGAAAGAAAAAGGGCTGACTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCg                                                                                                                       <  2:259666/149‑1 (MQ=255)

                                aCGCTGAACGAAAAAGGGCTGATTTTTGTTGGTGTGGATATCATCGGCGAGCGTCTGACTGTTATTTACGTCTCCATCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTGCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGa                                                                                                                      <  1:997794/149‑1 (MQ=255)

                                 agCTGAAAGAAAATGGGCGGATTTTTGTTGGTCTGTATCTGATCGGCGACCGGCTGACTGAAATTTACGTCTCCAGCCCTACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGGTAATGGATGCCATCGaa                                                                                                                     <  2:772486/148‑1 (MQ=255)

                                  gCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAg                                                                                                                    <  1:838460/149‑1 (MQ=255)

                                    tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                  <  1:1345118/149‑1 (MQ=255)

                                    tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                  >  1:193711/1‑149 (MQ=255)

                                     gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTACCTGCCCCTGCCCCACCTGTATTCGTGAGAGTGAAGCAGAGTTTCCGGGC‑TCGATCACCGGAATGTGAAAGGATGCCATCGAAGCAc                                                                                                                 >  1:1013198/1‑149 (MQ=255)

                                           aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTATTGGATGCCATCGAAGCACGTTTAc                                                                                                           >  2:30925/1‑149 (MQ=255)

                                           aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                           <  1:1209884/149‑1 (MQ=255)

                                            aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTATTGGATGCCATCGATGCACGGTTAca                                                                                                          >  1:1060023/1‑149 (MQ=255)

                                            aaaGGGCGGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAca                                                                                                          <  2:1080469/149‑1 (MQ=255)

                                             aaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                         <  2:196141/149‑1 (MQ=255)

                                               gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                       <  1:1507806/149‑1 (MQ=255)

                                               cggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCATAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                       <  2:228594/148‑1 (MQ=255)

                                                ggCTGATCTTTGTTGGTCTGGATATCATAGGCTACCGTCTGACTGAAATTAACGGCACCAGCCCCACCTGTCTTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGGTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                      <  2:1383020/149‑1 (MQ=255)

                                                 gCTTATTTTTGTTGGGCTGGATATCATCGGGGACCGTCTGACTGAAATTAACTTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTT‑TCGATCACCTTAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagc                                                                                                     <  2:887108/149‑1 (MQ=255)

                                                    gATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGt                                                                                                  >  1:1557866/1‑149 (MQ=255)

                                                    gATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGt                                                                                                  >  1:810504/1‑149 (MQ=255)

                                                      tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTaa                                                                                                <  2:1495070/149‑1 (MQ=255)

                                                      tttttGTTGGTCTGGATATCATCGGCGACCGTATGACTGAAATTAACTTCACCAGCCCAACCTGTATTCGTGAGATTTAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGAAGCAGTaa                                                                                                <  1:286300/149‑1 (MQ=255)

                                                       ttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCGCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCCCCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAc                                                                                               >  1:1476033/1‑149 (MQ=255)

                                                             tgGTCTGGATATCATCGGCGACCGTCTGACTGAAATAAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCt                                                                                         <  2:161690/149‑1 (MQ=255)

                                                                tCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAg                                                                                      >  2:1168045/1‑149 (MQ=255)

                                                                 cgggATTTCATCGGCGACCGTCTGACTGAAATTCACGCCACCAGCCCAACCTCTATTCGTGAGATTGAAGCAGAGTCTCCGGTG‑GCGCTCACCGGAATGTTACTGGAAGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGc                                                                                     <  2:488489/147‑1 (MQ=255)

                                                                 cTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAgg                                                                                     >  1:60837/1‑148 (MQ=255)

                                                                 cTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGc                                                                                     >  1:1001709/1‑149 (MQ=255)

                                                                      tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                <  1:208947/149‑1 (MQ=255)

                                                                              gCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGTATCACGTTTACAGCAGCAGTAACCCACCATAGCGAGAAGGATCTCg                                                                        >  2:407065/1‑149 (MQ=255)

                                                                                gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                      <  2:376508/149‑1 (MQ=255)

                                                                                  ccGTCTGACTGAAATTCACGTCACCAGCCCAACCTGTATTCGTGAGATAGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                    <  2:594552/149‑1 (MQ=255)

                                                                                    gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGACGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                  <  2:1465632/149‑1 (MQ=255)

                                                                                     tCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                                 >  1:834213/1‑149 (MQ=255)

                                                                                          ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTATGa                                                            >  1:511692/1‑149 (MQ=255)

                                                                                           tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCGCCGGAATGTTAATGGATGCAATCGAAGCACGTTTTCAGCAGCAGTACCCCACCTTAGCGAGACGGATCTCGGTGGGAATCTGAg                                                           >  2:742595/1‑149 (MQ=255)

                                                                                           tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                           >  1:991075/1‑149 (MQ=255)

                                                                                           tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTGAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                           >  1:1205961/1‑149 (MQ=255)

                                                                                           tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATCTTAATGGATGCCATCGCAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGat                                                           >  2:1520678/1‑148 (MQ=255)

                                                                                            taaaTTAACGTCTCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                          <  1:206883/148‑1 (MQ=255)

                                                                                            gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                          <  2:936537/149‑1 (MQ=255)

                                                                                            gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                          <  1:158745/149‑1 (MQ=255)

                                                                                            gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                          <  2:1242158/149‑1 (MQ=255)

                                                                                             aaaTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGCTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTg                                                         >  2:735057/1‑149 (MQ=255)

                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGAc                                                       >  2:1568679/1‑149 (MQ=255)

                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGAc                                                       >  2:670119/1‑149 (MQ=255)

                                                                                                      tCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGccc                                                <  1:1357641/149‑1 (MQ=255)

                                                                                                             cccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTc                                         <  1:314/149‑1 (MQ=255)

                                                                                                                aaCCAGTATTCTTTCGATGGAAGCAGAGTTTCCGGTG‑TCGATTACCGGAATGTTAATGGATGCCATCGAAGTACGTTTACAGCAGCAGTAACCCTCCTTAGCGAGATGGATCTCGTTTAGACTCTGAGGGACAGCGCCCTGCTTTCCAc                                      <  2:1067044/149‑1 (MQ=255)

                                                                                                                  ccTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                  >  2:1268299/1‑119 (MQ=255)

                                                                                                                   cTGTATTCGTGAGATTGAAGCATAGTTTCATGTG‑TCGATCACCGGAATGTTAATGGATTCCATCGAATCACGTTTAAAGCAGCAGTAACCCACCTAAGCTAGAAGTATCTAGTTGAGACACTGATTGACCGCGCCCTTCTTTCCACGCa                                   >  2:847260/1‑149 (MQ=255)

                                                                                                                        ttCGTGAGATTGACGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTg                              <  1:326539/149‑1 (MQ=255)

                                                                                                                            tGAGATTGAAGCAGACTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCAATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgca                          >  1:1172907/1‑148 (MQ=255)

                                                                                                                             gagaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcgc                         >  1:1239158/1‑149 (MQ=255)

                                                                                                                                    aaGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                  >  1:1056692/1‑149 (MQ=255)

                                                                                                                                    aaGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                  >  1:654539/1‑149 (MQ=255)

                                                                                                                                       ttaacattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCGTTCCACGCATACTGGGCGCTGTTGCtttt               >  1:1247933/7‑149 (MQ=255)

                                                                                                                                       cAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttt               >  2:152300/1‑149 (MQ=255)

                                                                                                                                        agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGATTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt              <  2:411526/149‑1 (MQ=255)

                                                                                                                                          acattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGGTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg            >  2:891281/4‑149 (MQ=255)

                                                                                                                                           gTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGa           <  1:1479463/149‑1 (MQ=255)

                                                                                                                                           gTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGa           <  1:1168045/149‑1 (MQ=255)

                                                                                                                                            tttCCGGTG‑TCGATCACCGTAATGTTAATGGATGCCATCGAAGCACGTTTCCAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa          <  1:151256/149‑1 (MQ=255)

                                                                                                                                            tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTACCCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCCCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTGTTTTGaa          >  2:742443/1‑149 (MQ=255)

                                                                                                                                            attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTCACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa          >  1:93600/2‑149 (MQ=255)

                                                                                                                                             ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGACGGACCTTGACCGCAGAAGTAACCCACCATAGCGAGAATGAACTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACTCATACTTGGACCTGTTGCTTTTTTGatc         >  1:1215279/1‑147 (MQ=255)

                                                                                                                                             ttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc         <  1:1549689/149‑1 (MQ=255)

                                                                                                                                              tCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAcc        <  2:282239/149‑1 (MQ=255)

                                                                                                                                                 ggTG‑TCGATCACCGGAATGTTAATGGATGCGATCGATGCAAGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGGATACTGGGCGCTGTTGCTTTTTTGAACCAgg     <  2:1529332/149‑1 (MQ=255)

                                                                                                                                                 ggTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAgg     <  1:1230685/149‑1 (MQ=255)

                                                                                                                                                    g‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGaaa  <  1:271200/149‑1 (MQ=255)

                                                                                                                                                    g‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGaaa  >  1:454991/1‑149 (MQ=255)

                                                                                                                                                    g‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGAAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCATTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGaaa  <  2:1187496/149‑1 (MQ=255)

GTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAA  >  NC_000913/3092626‑3092922

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