BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,092,779	Δ1 bp	7.8%	coding (902/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,779	0	T	.	7.8%	264.8 / 19.8	102	coding (902/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand): ref base T (41/53); new base . (4/4); total (45/57)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

Only 100 of 107 total aligned reads displayed.

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926

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)

  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:266709/149‑1 (MQ=255)

          cccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATg                                                                                                                                           >  2:607266/1‑149 (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)

               cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      <  1:212949/149‑1 (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)

                   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: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)

                                                          gTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCACCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCCTCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  2:1177969/1‑149 (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:1242158/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                             ccTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                      <  1:1268299/119‑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)

                                                                                                                                             gtg‑tTGTTCACCGGAATGTTAATGGATGCCATCGAAGCAAGTTTACAGCATCAGTAACCCACCTTAGCGAGAAGGTTGTCGTTGAGTCTCTGAGTGTGATCGCCGTTCTTTCCACGCATACTGGGCGCTGTTTCTTTTTTGAACCAGGa        <  2:18719/149‑1 (MQ=255)

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

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

                                                                                                                                                 tCGATCTCCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACGCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAAc     <  2:1476033/149‑1 (MQ=255)

                                                                                                                                                    aTCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  >  2:1563577/1‑149 (MQ=255)

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926

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