BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,238	T→A	6.1%	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	6.1%	630.1 / 15.6	217	intergenic (+16/‑131)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base T (104/92); new base A (14/3); total (121/96)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781091‑781387

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

 cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAACCTGCACGACCCACCAGTTTTAACATCAAACTCAg                                                                                                                                                    >  1:243017/1‑150 (MQ=255)

  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                   >  1:262473/1‑150 (MQ=255)

  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                   >  1:4672990/1‑150 (MQ=255)

   ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAt                                                                                                                                                  >  1:5483155/1‑150 (MQ=255)

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

    cTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTATCATCAAACTCAGATg                                                                                                                                                 >  1:5354731/1‑150 (MQ=255)

    cTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATg                                                                                                                                                 >  1:3297402/1‑150 (MQ=255)

    cTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACTTCAAAATCAGATg                                                                                                                                                 >  2:1077605/1‑150 (MQ=255)

        cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTaa                                                                                                                                             <  2:3363514/150‑1 (MQ=255)

        cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTaa                                                                                                                                             <  2:4848878/150‑1 (MQ=255)

         ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTGAAAGTCAGATTTTTAg                                                                                                                                            >  2:2581591/1‑150 (MQ=11)

         ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTATAACATCAAACTCAGATGTTAAg                                                                                                                                            >  1:4376168/1‑150 (MQ=255)

          tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                           <  2:1628089/150‑1 (MQ=255)

             aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTatgatgc                                                                                                                                        >  1:3029269/1‑144 (MQ=255)

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

             aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTAAACATCAAACTCAGATGTTAAGCGTg                                                                                                                                        >  2:86228/1‑150 (MQ=255)

             aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAg                                                                                                                                            >  2:1275757/1‑145 (MQ=17)

               tgtgggggTCGGGGGTTTGTTACCGTCAACAACCAACAACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                      <  1:4952559/147‑1 (MQ=11)

                gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCTGTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGTCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                     <  1:3239470/150‑1 (MQ=255)

                gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGATAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATAAAACAAACAAGATAAGCGTGAAg                                                                                                                                     >  1:3722780/1‑150 (MQ=14)

                  gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGa                                                                                                                                   <  1:1033833/150‑1 (MQ=255)

                   ggggTCGGTGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGTCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                  >  2:869002/1‑150 (MQ=255)

                   ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATAAAACTCAGAATTTAACCGTGAAGGAt                                                                                                                                  >  1:2587993/1‑150 (MQ=14)

                   ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                  <  1:1271336/150‑1 (MQ=255)

                     ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                                <  1:3634545/150‑1 (MQ=255)

                             gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATAAAACACAGATGTTAAGCGTGAAGGATAAAGTTgcgc                                                                                                                        >  1:4473123/1‑150 (MQ=255)

                              ttCGCTCCCGTAATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                       <  1:3382235/150‑1 (MQ=255)

                               tCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCCGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGTGAAGGATAACGTTGCGCCa                                                                                                                      >  1:224614/1‑150 (MQ=255)

                                cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                     <  1:1573491/150‑1 (MQ=255)

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

                                  aTCCCGTCATCACCCACCACTCGGGTCGTTAGCACAGTTGGTAGATCAGTTTACTTTTAATAAATTTGTCGCAGGTTAGAATCCTGCACGAACCACCAATTTTAACAACAAAATAAGATGTTAATCGGGAAGTATAACGTTGCGACAgaa                                                                                                                   >  1:1170411/1‑148 (MQ=14)

                                   tCCCGTCAGCA‑CCACCCCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAAc                                                                                                                 <  2:3715072/150‑1 (MQ=255)

                                     ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGAAGTTGAATTTTAATCAATTTGTCGCAGGTTCGCATCCTGCACGACCCACCAGTTTTAACATAAAACTAAGTAGATAAGAGTGATGGAAAACGTGTCGCcaacaaac                                                                                                                >  1:3349580/1‑148 (MQ=14)

                                     ccGTCATCACCCACCACTCGGGTCGTTAGCGCAGGTGGTAGAGCCGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATACTGCACGACCCCCCCCCTTTAACATCCAACGCAGGGGGTAAGCGTGAAGGATAACGGGGCGCCAGccccg                                                                                                                >  1:1939074/1‑146 (MQ=255)

                                      cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGTTGTTAAGCGGGAAGGATAACGTTGCGCCAGCAACgg                                                                                                               >  1:1770138/1‑150 (MQ=255)

                                        tcatcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAACCACCAGTTTTAACATCAAAATAAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGcc                                                                                                             >  1:3056582/1‑150 (MQ=255)

                                            cccccaccACTCGGGTCGTTAGGTCAGTTGGTAGAGCTGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                         <  2:4694444/148‑1 (MQ=255)

                                             acccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTACCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAg                                                                                                        <  2:2448430/150‑1 (MQ=255)

                                             acccaccA‑TCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                               <  1:1829300/126‑1 (MQ=255)

                                              cccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACAGCACGACCCACCAGTTTTAACATCAAACTCAGTTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCgaggg                                                                                                       >  1:3534830/1‑146 (MQ=255)

                                               ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTAGCAGGTTAGAATACTGCACGACCCACCAGTTTTAAAATCAAAATAAGATGTTAAGCGTGAAGGATAAAGTTGCGCCAGCAACCGCCCCTAggg                                                                                                      >  1:604679/1‑150 (MQ=14)

                                                caccacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgc                                                                                                     <  2:4092473/150‑1 (MQ=255)

                                                      tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAACTAAGAAGTTAAGCGTCAAGCAAAACGGTGCGCCAGCAACGGCCCCGAGCgcgagcga                                                                                               >  2:2694042/1‑150 (MQ=14)

                                                           tCGTTAGCACAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTATTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                          <  2:3620958/150‑1 (MQ=255)

                                                           gcGTTCGCTCAGTTGGTAGAGCAGTTGACTATTAATCCATCGGACGCAGGATCGAATCCTGCACGACCCACCTGATATATCATCAAACTCAGATGTTTAGCGTGTAGGATAACGTTGCGCCAGCAACGTCCCGTAGGGCGagcgaagcga                                                                                          <  2:3977862/149‑1 (MQ=11)

                                                               tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAACACAGATGTTAAGCGGGACGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGCGTCa                                                                                      >  1:477841/1‑150 (MQ=11)

                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                             <  2:1916964/139‑1 (MQ=255)

                                                                      gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCATTTTAAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTCCGCCAGCAACGCCCCGTAGGCCTAGCGAAGAGAGTCaaacagcc                                                                               >  2:399511/1‑143 (MQ=14)

                                                                      gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                               >  1:2696488/1‑150 (MQ=255)

                                                                         ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                            >  2:512609/1‑150 (MQ=255)

                                                                          gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAc                                                                           >  2:4565217/1‑150 (MQ=255)

                                                                             gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCTAAGCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                                                                        <  1:3118850/150‑1 (MQ=255)

                                                                              agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                         <  2:2745890/116‑1 (MQ=255)

                                                                              agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCAGCCCGAccccc                                                                       >  2:3909004/1‑148 (MQ=255)

                                                                              agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                         >  1:2745890/1‑116 (MQ=255)

                                                                               gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                                      >  1:1367853/1‑150 (MQ=255)

                                                                               gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                                      >  1:1367834/1‑150 (MQ=255)

                                                                                   tttaattttAATCATTGGGTCGCTGTTTCGAATCCGGCACGACCCACCATTTTTAACAACAAACTCTTATTTTAAGCGTGATGGATAACGTTGCGCCACCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATc                                                                  <  1:5082784/145‑1 (MQ=14)

                                                                                    tGACTTTTAATCTTTTGGTCGCTGGTTCGAATCCTGCACGACCCACCTGTTTTAACATCAAACTCAGTTGTTAAGCGTTAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                                 <  1:3406398/150‑1 (MQ=14)

                                                                                    tGAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                                 <  1:5202531/150‑1 (MQ=255)

                                                                                        tttttttCATTTGGTCGCAGGTTCGATTCCTGCTGGTCCCACCTGTTTTAACATCAAACTCTGTTGTTATGGGTGTAGGATAAGGTTGCGCCTGCTACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  2:477841/144‑1 (MQ=14)

                                                                                        ttttAATCATTTGGTCGCTGGTTCGAATCCTGCACGACCCACCTGTTTTAACATCAAACTAAGATGTTATGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  1:2043380/150‑1 (MQ=17)

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

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  2:3786772/150‑1 (MQ=255)

                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAAATCAGAAGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATACAGCACCACCCACCAATCATaaa                                                             >  2:3019327/1‑150 (MQ=14)

                                                                                              tCAATTGGTCGCGGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTACAGATTgg                                                       >  1:5072270/1‑150 (MQ=255)

                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                      >  1:3735746/1‑71 (MQ=255)

                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                      <  2:3735746/71‑1 (MQ=255)

                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                       >  1:1888361/1‑150 (MQ=255)

                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                       >  1:3623800/1‑150 (MQ=255)

                                                                                               cAATTGGGCGCTGTTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCATATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggt                                                      <  1:1287053/150‑2 (MQ=14)

                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggtgg                                                    <  2:3269119/150‑4 (MQ=255)

                                                                                                  ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                   <  2:1341443/150‑1 (MQ=255)

                                                                                                  ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                   >  1:1196253/1‑150 (MQ=255)

                                                                                                    ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                 <  2:4121897/150‑1 (MQ=255)

                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGTAGCGAGTCAGACTGCACGACCCACCAAGCTTAGAGTGTGGCCCCGAGt                                               >  2:2631744/1‑150 (MQ=14)

                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  1:146132/1‑150 (MQ=255)

                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGAGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  2:2964245/1‑150 (MQ=255)

                                                                                                             ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        >  1:4527223/1‑150 (MQ=255)

                                                                                                             ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGGCGAGCGGAGCGAGTCATCCTGCACGACCCACCAAGCTTAGAGTTTGGCCCCGAGTAAAAATc                                        >  2:4909583/1‑150 (MQ=18)

                                                                                                              tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                       >  1:804214/1‑150 (MQ=255)

                                                                                                              tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                       >  1:3804495/1‑150 (MQ=255)

                                                                                                              tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTGAAGCGGGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCCAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGACTAAAAATCt                                       >  1:420039/1‑150 (MQ=255)

                                                                                                                gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     >  1:2150587/1‑150 (MQ=255)

                                                                                                                gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  2:4736470/150‑1 (MQ=255)

                                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                    <  1:2763116/150‑1 (MQ=255)

                                                                                                                  aTCCTGCTCTACCCACCAGTTTTTACATCTAACTCAGTTTTTAATCTTTTTGGATAACGTTGCGACTGCATCGCCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:5422428/150‑1 (MQ=18)

                                                                                                                  aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  1:4986441/150‑1 (MQ=255)

                                                                                                                   tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  2:5483155/150‑1 (MQ=255)

                                                                                                                    ccTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTTAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                 <  1:1962627/150‑1 (MQ=255)

                                                                                                                       gCACGACCAACCCGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATATCGTTGCGCCTGCAACGGCCCGTAGGGCGAGCGAAGCGAGTGTTCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  2:5322567/150‑1 (MQ=21)

                                                                                                                        cACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                             >  2:5216048/1‑150 (MQ=255)

                                                                                                                              ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGt                       >  1:473240/1‑150 (MQ=255)

                                                                                                                                accaGTTTTATCATCGAAGTCAGATGTTAAGCGTGAAGGATATCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     <  2:2905439/150‑1 (MQ=255)

                                                                                                                                 ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         >  1:1311776/1‑129 (MQ=14)

                                                                                                                                 ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                    <  1:4650165/150‑1 (MQ=255)

                                                                                                                                        tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   >  1:232582/1‑128 (MQ=14)

                                                                                                                                        tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:232582/128‑1 (MQ=255)

                                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAACCTGAACGACCCACAAATCAAACAGAATGGCCCCGACTAAAAATCTATCAGCAAACACACGAATGGGTCGTTAGCTc        >  2:4379346/1‑150 (MQ=18)

GCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781091‑781387

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