BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,238	T→A	5.3%	intergenic (+16/‑131)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,238	0	T	A	5.3%	670.7 / 10.1	212	intergenic (+16/‑131)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base T (103/94); new base A (10/4); total (114/98)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

GCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781091‑781387

tcaacTCCCTTACCATGTTGGGGTTGGGGGTTCGTTCCCGTCATCAACCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCa                                                                                                                                                      <  1:2770419/146‑1 (MQ=14)

 cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAg                                                                                                                                                     >  2:42494/1‑150 (MQ=40)

  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                    <  1:1904692/150‑1 (MQ=255)

   ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAt                                                                                                                                                   <  2:1954401/150‑1 (MQ=255)

      cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGtt                                                                                                                                                >  1:2819705/1‑150 (MQ=255)

       ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTa                                                                                                                                               >  2:2676156/1‑150 (MQ=255)

         ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTAAACCTCCAACACAGaagagaac                                                                                                                                             >  2:3595675/1‑143 (MQ=11)

             aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTg                                                                                                                                         >  1:3204836/1‑150 (MQ=255)

             aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACTAAGATGTTAAGCGTg                                                                                                                                         >  1:3292634/1‑150 (MQ=255)

                gaggGGGTCGTCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCATTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                      <  1:409025/150‑1 (MQ=255)

                  gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGa                                                                                                                                    <  2:1685088/150‑1 (MQ=255)

                    gggCCGGCGCCTCGATCCTGTCATCAGCCTCCCCTCGGGTCGTTCGACCATTTGCTATAGCATATGTCTCTACATGAATGGGTGGCCTTCTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                  <  2:2156744/150‑1 (MQ=9)

                      gTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                                <  1:1234008/150‑1 (MQ=255)

                         ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTAGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGTGAATGATAACGtt                                                                                                                             >  1:1701171/1‑150 (MQ=255)

                           cggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgc                                                                                                                           <  1:42575/150‑1 (MQ=255)

                            ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCGGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAACCTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                          <  2:1967592/150‑1 (MQ=255)

                             gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcgc                                                                                                                         >  1:3587989/1‑150 (MQ=255)

                             gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcgc                                                                                                                         <  1:1267769/150‑1 (MQ=255)

                              ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCGGATGTTAAGCGGGAAGGATAACGTGGCGcc                                                                                                                        >  1:1656800/1‑150 (MQ=255)

                                 gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                     >  1:3027117/1‑150 (MQ=255)

                                  aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCAAGCa                                                                                                                    >  2:239613/1‑150 (MQ=255)

                                  aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAAACAGATGTTAAGCGTGAAGGATAGCGTAGCGCCAGCa                                                                                                                    >  2:3927968/1‑150 (MQ=255)

                                     ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAACGGGGCGCCAGCAACg                                                                                                                 >  1:3258221/1‑150 (MQ=255)

                                      cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgg                                                                                                                >  1:3782533/1‑150 (MQ=255)

                                      cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTATTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgg                                                                                                                <  2:1118373/150‑1 (MQ=255)

                                        tcttcaCCCACCACTCTGGTCTTTAGCTCTGTTGGTTGAGCTTTTGACTTTTAATCAATTTGTCTATGGTTCTAATCCTTCACTACCAACCTGTTTTAAATTCAAACTCTGATGTTAAGCGTGAAGGATAAGGTTGCGCAAGCAACGGcc                                                                                                              <  2:3808181/147‑1 (MQ=255)

                                           tcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCAACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGATGTATAACGTTGCGCCAGCAACGGCCCGt                                                                                                           <  2:698138/150‑1 (MQ=255)

                                               ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTTAAGCGTGACGGGTAACGTTGCGCCAGCACCGGCCCGTAggg                                                                                                       >  2:2816124/1‑150 (MQ=255)

                                               ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACAAAGAAGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAggg                                                                                                       >  1:2336865/1‑150 (MQ=255)

                                                caccacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCTGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgc                                                                                                      >  2:1184100/1‑150 (MQ=255)

                                                            cGTTAGCTCAGTTGGTTGAGAACAACACTCTTAAACAATTGGCCGCAGGATCGCATCCTGCAATACCCACCA‑TTTTAAACGCCTACACTGGTTTATATCGTGAAGGATAACGTTGCGATAGCAACGGCCCGTAGGGCGAgacgaagtcga                                                                                         <  2:1467146/150‑5 (MQ=2)

                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                            <  2:2876616/98‑1 (MQ=40)

                                                                 gCTCAGTTGGTAGAGCAGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATc                                                                                    >  2:3307899/1‑150 (MQ=255)

                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTTAAGCTTGAAGGATAACTTTGCTAAAGCAACGCCCCGAAGGCCTAG‑CGAAGACAGTCaaaca                                                                                  >  2:3883705/1‑146 (MQ=14)

                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTAAACAACAAACTCAGATGTTAAGCGTGACGGATAACGTAGAGCAAGCACCGGCACTAAGTGCGAG‑CGAACCGAGTCaaaat                                                                                  >  2:3251381/1‑146 (MQ=14)

                                                                     aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGc                                                                                >  2:3558563/1‑150 (MQ=255)

                                                                       ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                            <  1:539665/89‑1 (MQ=255)

                                                                       ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                            >  2:539665/1‑89 (MQ=255)

                                                                         ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGa                                                                            <  2:3727436/150‑1 (MQ=255)

                                                                           tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGAcc                                                                          <  2:102906/150‑1 (MQ=255)

                                                                              tgCAGTTGACTTTTATTCCATTGGTCGCAGGTTCGAAGCCTGCACGACCCACCAGTTTTAACACCAAACTCAGTTTTTATGGGTGTTGGTTAACTTTGCGCCTCCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACccac                                                                       <  1:3239684/149‑1 (MQ=11)

                                                                              agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACccac                                                                       >  1:3135118/1‑150 (MQ=255)

                                                                                      aCTTTTAATCAAGTGGTCGCAGGTTCGAATCCGGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTa                                                               <  2:732762/150‑1 (MQ=255)

                                                                                        ttttAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  1:3307899/150‑1 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCTGTTTCTAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGATGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  1:1357833/150‑1 (MQ=17)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAg                                                                                                                                                     <  1:126058/63‑1 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  1:1116984/1‑150 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  1:1116888/1‑150 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCATaaa                                                             >  2:2522849/1‑150 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAAACTGaaa                                                             >  1:166961/1‑150 (MQ=255)

                                                                                             aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                        >  1:3645139/1‑150 (MQ=255)

                                                                                             aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                        <  2:3301262/150‑1 (MQ=255)

                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                       >  2:1658036/1‑150 (MQ=255)

                                                                                               cAATTGGCCGCAGTTTCGAATCCTGCACGACCCACCAGTTTTAACAACAAACTCAGTTGTTAAGCGTGATGGATACCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                      <  1:3431953/150‑1 (MQ=17)

                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGccc                                                    >  1:580311/1‑150 (MQ=255)

                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGccc                                                    <  2:183478/150‑1 (MQ=255)

                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAAGATTGGccc                                                    >  1:2653878/1‑150 (MQ=255)

                                                                                                     gTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAg                                                >  2:1211898/1‑150 (MQ=11)

                                                                                                         ctggTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                            <  1:658206/148‑1 (MQ=255)

                                                                                                          aGGCTCGCATCCCGCACGACCCACCAGTTTTAACATAACCGTCAGCAGCTATTCGTGAAGTATTACGTTGCGACTGCAACGGCCCGTAGGTCGTG‑CGATGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                           <  2:1666305/150‑1 (MQ=17)

                                                                                                            gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         <  2:3408298/150‑1 (MQ=255)

                                                                                                            gTTCGAATCCTGCACGACCCACCAGTTTTAACAGCAAACTCAGAGGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         >  2:143719/1‑150 (MQ=255)

                                                                                                             ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        >  2:1446669/1‑150 (MQ=255)

                                                                                                             ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        <  1:3390822/150‑1 (MQ=255)

                                                                                                                gAATCCTGCACTACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATATCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGTGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  2:3560317/150‑1 (MQ=255)

                                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                    >  1:1816272/1‑150 (MQ=255)

                                                                                                                 aaTCCTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATACTGAACGACCCACCAATCTTAAAGATTGGAGCCGAGTAAAAATCTATc                                    >  2:3094782/1‑150 (MQ=18)

                                                                                                                  aTCCTGCACGACCCACCTGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAATAATCTTTCa                                   <  1:2966671/150‑1 (MQ=255)

                                                                                                                  aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:2054640/150‑1 (MQ=255)

                                                                                                                     cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                >  1:4067090/1‑150 (MQ=255)

                                                                                                                      tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                               >  1:2221778/1‑150 (MQ=255)

                                                                                                                      tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                               >  1:901552/1‑150 (MQ=255)

                                                                                                                          cGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                           <  1:2005397/150‑1 (MQ=255)

                                                                                                                             cccACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCg                        >  2:1327204/1‑150 (MQ=255)

                                                                                                                               caccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                      <  2:3027117/150‑1 (MQ=255)

                                                                                                                                accaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     <  2:3797303/150‑1 (MQ=255)

                                                                                                                                accaGTTTTAACATCAAACTCAGTTGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTATAAATCTTTCTGGTAACACCCGTAt                     <  1:2897043/150‑1 (MQ=255)

                                                                                                                                accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     <  1:2676156/150‑1 (MQ=255)

                                                                                                                                accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCAGCACGACCCACCAAACTAAAAGATTTGCCCCGATAAAAAATCTTTAAGGTAACACCCGTAt                     >  1:1467146/1‑150 (MQ=18)

                                                                                                                                 ccaGTTTTAACATCTAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                    <  2:3761878/150‑1 (MQ=255)

                                                                                                                                 ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                    <  1:1435307/150‑1 (MQ=255)

                                                                                                                                  caGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGGAAAAATATTTAAGTAACCAccccaatat                   >  2:3176842/1‑144 (MQ=11)

                                                                                                                                  caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                   <  2:2225817/150‑1 (MQ=255)

                                                                                                                                  caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                   <  1:3558563/150‑1 (MQ=255)

                                                                                                                                   aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                  >  1:2525943/1‑150 (MQ=255)

                                                                                                                                   aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                  >  1:256285/1‑150 (MQ=255)

                                                                                                                                   aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                  >  1:1219793/1‑150 (MQ=255)

                                                                                                                                     ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTGTCAGGTAACACCCGTATGGGTc                >  1:3908/1‑150 (MQ=255)

                                                                                                                                       ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  2:3292570/150‑1 (MQ=255)

                                                                                                                                            aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAgtt         >  1:1138375/1‑148 (MQ=255)

                                                                                                                                            aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt         <  1:3459781/150‑1 (MQ=255)

                                                                                                                                            aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAGCTTTCAGGTAACACCCGGATGGGTCGTTAGCt         >  1:1484696/1‑150 (MQ=255)

                                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATTTTTCAGGTAACACCCGTATGGGTCGTTAGCTc        >  1:1845998/1‑150 (MQ=255)

                                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc        >  1:2071075/1‑150 (MQ=255)

                                                                                                                                                aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTACAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  2:1064222/150‑1 (MQ=255)

                                                                                                                                                aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  2:1116984/150‑1 (MQ=255)

                                                                                                                                                aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  2:1116888/150‑1 (MQ=255)

                                                                                                                                                 aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGTCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt    >  2:585332/1‑150 (MQ=255)

                                                                                                                                                   tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg  >  1:2909001/1‑150 (MQ=255)

                                                                                                                                                   tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGGTgg  >  2:499716/1‑150 (MQ=255)

GCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781091‑781387

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