BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,247,843:1	+G	5.8%	coding (71/663 nt)	yqjA →	DedA family protein YqjA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,247,843	1	.	G	5.8%	424.0 / 21.8	103	coding (71/663 nt)	yqjA	DedA family protein YqjA
Reads supporting (aligned to +/- strand): ref base . (36/61); new base G (4/2); total (40/63)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 109 total aligned reads displayed.

TACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGA  >  NC_000913/3247695‑3247986

tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑g                                                                                                                                                <  2:1672571/150‑1 (MQ=255)

tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑g                                                                                                                                                >  2:946715/1‑150 (MQ=255)

 aCCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑ga                                                                                                                                               >  1:1427960/1‑150 (MQ=255)

  ccAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑gat                                                                                                                                              >  1:852625/1‑150 (MQ=255)

  ccAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑gat                                                                                                                                              >  2:1928044/1‑150 (MQ=255)

    aaaaCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGGAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTg                                                                                                                                            >  1:1314425/1‑150 (MQ=255)

       aCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCa                                                                                                                                         <  2:106002/150‑1 (MQ=255)

        cAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAAT‑CATCGATGGATTGGcca                                                                                                                                        <  1:989561/150‑3 (MQ=255)

          gCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCAtgt                                                                                                                                      >  1:1813207/1‑150 (MQ=255)

          gCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAAT‑CATCGATGGATTGGccagg                                                                                                                                      <  2:1935410/150‑5 (MQ=255)

           cAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCAtgtt                                                                                                                                     >  2:1786194/1‑150 (MQ=255)

           cAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCAtgtt                                                                                                                                     <  1:1204299/150‑1 (MQ=255)

              tAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTa                                                                                                                                  >  2:1182856/1‑150 (MQ=255)

                aCATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTAtt                                                                                                                                >  2:1800772/1‑150 (MQ=255)

                 cATTGCGGTAGTGCATCTTTTAATACCAGCGTGGCGTTAACCGATTCACCAGGAATATTGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTAttt                                                                                                                               <  2:880506/150‑1 (MQ=255)

                   ttGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTg                                                                                                                             <  2:883173/150‑1 (MQ=255)

                     gCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTc                                                                                                                           >  1:964139/1‑150 (MQ=255)

                     gCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTc                                                                                                                           >  1:1185557/1‑150 (MQ=255)

                      cGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCt                                                                                                                          >  2:872942/1‑150 (MQ=255)

                          aGTGCGTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtt                                                                                                                      <  1:636461/150‑1 (MQ=255)

                          aGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtt                                                                                                                      >  2:709150/1‑150 (MQ=255)

                          aGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtt                                                                                                                      <  1:792810/150‑1 (MQ=255)

                          aGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtt                                                                                                                      <  1:202642/150‑1 (MQ=255)

                            tGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttg                                                                                                                    <  2:1917682/150‑1 (MQ=255)

                             gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttgt                                                                                                                   <  2:1596887/150‑1 (MQ=255)

                             gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttgt                                                                                                                   <  2:709134/150‑1 (MQ=255)

                              cATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTa                                                                                                                  <  1:1786194/150‑1 (MQ=255)

                              cATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTa                                                                                                                  <  1:477581/150‑1 (MQ=255)

                                tCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                                                                                                >  1:282173/1‑150 (MQ=255)

                                     tAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGt                                                                                                           <  2:221619/150‑1 (MQ=255)

                                        aaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTc                                                                                                         >  1:1988656/1‑149 (MQ=255)

                                        aaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTc                                                                                                         <  2:1988656/149‑1 (MQ=255)

                                          ccAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgt                                                                                                                       >  1:783770/1‑133 (MQ=255)

                                          ccAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgt                                                                                                                       <  2:783770/133‑1 (MQ=255)

                                           cAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                                                                     <  1:576017/150‑1 (MQ=255)

                                            aGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGa                                                                                                    <  2:1978002/150‑1 (MQ=255)

                                              cGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaaa                                                                                                  <  1:590981/150‑1 (MQ=255)

                                               gTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                                                                                                <  2:1667268/135‑1 (MQ=255)

                                               gTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                                                                                                >  1:1667268/1‑135 (MQ=255)

                                                 ggCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACg                                                                                               >  1:1046997/1‑150 (MQ=255)

                                                   cGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                                                                                             <  2:282173/150‑1 (MQ=255)

                                                     ttAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgctt                                                                                           <  2:852625/150‑1 (MQ=255)

                                                          cGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTc                                                                                      <  1:1359275/150‑1 (MQ=255)

                                                                 ccAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                                                                               <  1:1182856/150‑1 (MQ=255)

                                                                 ccAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                                                                               >  1:178402/1‑150 (MQ=255)

                                                                    ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCtt                                                                            >  1:242774/1‑150 (MQ=255)

                                                                    ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTAGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCtt                                                                            <  2:469164/150‑1 (MQ=255)

                                                                          atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg                                                                                   <  1:1927723/137‑1 (MQ=255)

                                                                          atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg                                                                                   <  1:1535976/137‑1 (MQ=255)

                                                                          atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg                                                                                   >  2:1927723/1‑137 (MQ=255)

                                                                          atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg                                                                                   >  2:1535976/1‑137 (MQ=255)

                                                                            gaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCgg                                                                    >  1:1382228/1‑150 (MQ=255)

                                                                             aatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCAGGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCggg                                                                   <  1:1930280/150‑1 (MQ=255)

                                                                                ggAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGa                                                                >  2:1962153/1‑150 (MQ=255)

                                                                                      tttGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGttt                                                          <  2:1427960/150‑1 (MQ=255)

                                                                                         gACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACt                                                       >  2:37644/1‑150 (MQ=255)

                                                                                         gACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGGGTTGTTTGTAATTTTGTTCGTTGAAAACGGCTTGGTTCCGGCGGCCTTTTTACCGGGCGACAGTTTagt                                                       >  2:374602/1‑148 (MQ=255)

                                                                                           cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGTATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                     >  2:1687332/1‑150 (MQ=255)

                                                                                           cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                     <  2:932427/150‑1 (MQ=255)

                                                                                             cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTa                                                   <  2:835071/150‑1 (MQ=255)

                                                                                             cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTa                                                   >  1:1121510/1‑150 (MQ=255)

                                                                                               aTTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTAtt                                                 <  2:95882/150‑1 (MQ=255)

                                                                                                 tgctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTgg                                               <  1:652985/150‑1 (MQ=255)

                                                                                                    tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGtt                                                           >  2:1764974/1‑135 (MQ=255)

                                                                                                    tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGtt                                                           <  1:1764974/135‑1 (MQ=255)

                                                                                                    tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCg                                            <  2:242774/150‑1 (MQ=255)

                                                                                                     gcAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCgg                                           <  2:1314425/150‑1 (MQ=255)

                                                                                                     gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCgg                                           >  2:1826563/1‑150 (MQ=255)

                                                                                                     gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCgg                                           <  1:37644/150‑1 (MQ=255)

                                                                                                       aaGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCg                                         <  2:528727/150‑1 (MQ=255)

                                                                                                         gCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCgtg                                       <  2:1046997/150‑1 (MQ=255)

                                                                                                         gCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCgtg                                       <  1:1099260/150‑1 (MQ=255)

                                                                                                          cccTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCTTCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCgtgt                                      <  1:1490732/150‑1 (MQ=255)

                                                                                                           ccTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGtt                                     <  1:1826563/150‑1 (MQ=255)

                                                                                                           ccTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGtt                                     <  2:1598664/150‑1 (MQ=255)

                                                                                                            cTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCTTGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTg                                    <  1:772756/150‑1 (MQ=255)

                                                                                                            cTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTg                                    <  2:1935471/150‑1 (MQ=255)

                                                                                                              gtgGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGAt                                  <  1:449750/150‑1 (MQ=255)

                                                                                                                 ggCGCAGGATTTTGAAACCCTGGGCAATCCATCGAT‑GATTGGCATGTTGTATTTTTTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGc                               <  2:725747/150‑1 (MQ=255)

                                                                                                                 ggCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGc                               <  1:466728/150‑1 (MQ=255)

                                                                                                                  gcgcAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCg                              <  2:1593407/150‑1 (MQ=255)

                                                                                                                      aGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAg                          <  1:1928044/150‑1 (MQ=255)

                                                                                                                       ggATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAgg                         >  1:1519111/1‑150 (MQ=255)

                                                                                                                           tttGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                              >  2:1885654/1‑125 (MQ=255)

                                                                                                                           tttGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGcgcg                     <  2:1121510/150‑1 (MQ=255)

                                                                                                                             tGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGAt                   >  2:1152936/1‑150 (MQ=255)

                                                                                                                             tGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGAt                   >  2:182414/1‑150 (MQ=255)

                                                                                                                              gAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGGGTTGATTGCGAAAGGCGCGATg                  <  2:685861/150‑1 (MQ=255)

                                                                                                                               aaaCCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATgg                 <  1:1800772/150‑1 (MQ=255)

                                                                                                                                acaacatgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                >  1:1655338/8‑150 (MQ=255)

                                                                                                                                aaCCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                >  1:536986/1‑150 (MQ=255)

                                                                                                                                 aCCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGc               <  1:1031792/150‑1 (MQ=255)

                                                                                                                                   ccTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTa             >  1:602586/1‑150 (MQ=255)

                                                                                                                                    catgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGATAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTAt            >  1:1885948/4‑150 (MQ=255)

                                                                                                                                    cTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTAt            >  2:1186426/1‑150 (MQ=255)

                                                                                                                                     atgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATc           >  2:1779445/3‑150 (MQ=255)

                                                                                                                                      ggCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATcc          >  2:74684/1‑150 (MQ=255)

                                                                                                                                            tCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAc    <  1:182414/150‑1 (MQ=255)

                                                                                                                                            tCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAc    <  2:765900/150‑1 (MQ=255)

                                                                                                                                              cATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGa  >  2:1187272/1‑150 (MQ=255)

TACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGA  >  NC_000913/3247695‑3247986

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