BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,493,563	T→A	8.2%	intergenic (‑309/+204)	frc ← / ← yfdX	formyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,493,563	0	T	A	8.2%	371.9 / 19.7	110	intergenic (‑309/+204)	frc/yfdX	formyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (44/56); new base A (4/5); total (48/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 120 total aligned reads displayed.

TTTCTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAA  >  NC_000913/2493415‑2493711

tttCTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGCGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCttttt                                                                                                                                                       >  1:806295/1‑149 (MQ=255)

 ttCTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCttttt                                                                                                                                                       >  4:64815/1‑148 (MQ=255)

   cTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTaa                                                                                                                                                     <  3:1048694/148‑1 (MQ=255)

    tGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAg                                                                                                                                                   <  2:172250/149‑1 (MQ=255)

    tGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAg                                                                                                                                                   >  2:756243/1‑149 (MQ=255)

     gCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAg                                                                                                                                                   <  4:48705/148‑1 (MQ=255)

      cATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGc                                                                                                                                                  <  3:114085/148‑1 (MQ=255)

      cATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCg                                                                                                                                                 <  1:756243/149‑1 (MQ=255)

        tGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGt                                                                                                                                                <  4:591772/148‑1 (MQ=255)

        tGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGt                                                                                                                                                <  4:749598/148‑1 (MQ=255)

         ggTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑a                                                                                                                                              <  4:161712/148‑1 (MQ=255)

         ggTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGT‑aa                                                                                                                                             <  2:939105/149‑2 (MQ=255)

           tAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑Acc                                                                                                                                            <  3:635479/148‑1 (MQ=255)

            aaGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCAt                                                                                                                                          <  2:333162/149‑1 (MQ=255)

                   gAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTCAAAGCGT‑ACCATAAaaca                                                                                                                                    <  3:192251/148‑1 (MQ=255)

                     ggACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataa                                                                                                                                 <  1:324815/149‑1 (MQ=255)

                     ggACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacata                                                                                                                                  <  3:739249/148‑1 (MQ=255)

                      gACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataac                                                                                                                                >  1:580443/1‑149 (MQ=255)

                      gACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataa                                                                                                                                 >  4:305434/1‑148 (MQ=255)

                         gAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAacataacata                                                                                                                             >  2:387691/1‑149 (MQ=255)

                           aaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataa                                                                                                                                 >  2:873169/1‑143 (MQ=255)

                           aaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataa                                                                                                                                 <  1:873169/143‑1 (MQ=255)

                              gaagaaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACAtaataa                                                                                                                         >  3:664169/1‑148 (MQ=255)

                              gaagaaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATaaa                                                                                                                        >  2:728201/1‑149 (MQ=255)

                               aagaaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAAt                                                                                                                       <  2:83504/149‑1 (MQ=255)

                               aagaaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGt                                                                                                                                                >  3:801306/1‑125 (MQ=255)

                               aagaaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGt                                                                                                                                                <  4:801306/125‑1 (MQ=255)

                                  aaTATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGa                                                                                                                     <  4:490257/148‑1 (MQ=255)

                                     atGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGCGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGAttt                                                                                                                  <  4:303028/148‑1 (MQ=255)

                                        aTAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTtgtg                                                                                                               <  3:233698/148‑1 (MQ=255)

                                         tAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGtga                                                                                                              <  3:305434/148‑1 (MQ=255)

                                         tAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGtga                                                                                                              >  4:943725/1‑148 (MQ=255)

                                            atatTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGtgatga                                                                                                           <  3:253382/148‑1 (MQ=255)

                                               ttGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTaa                                                                                                                                                     <  3:309405/104‑1 (MQ=255)

                                               ttGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTaa                                                                                                                                                     >  4:309405/1‑104 (MQ=255)

                                               ttGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACg                                                                                                       >  2:317953/1‑149 (MQ=255)

                                               ttGATTCGTGGTAATCAGAGAGAGTCGCGACAATTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGacac                                                                                                        >  3:670466/1‑148 (MQ=255)

                                                      gTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGatta                                                                                                 <  4:206616/148‑1 (MQ=255)

                                                        ggTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTg                                                                                              >  2:57340/1‑149 (MQ=255)

                                                           aaTCAGTGAGAGTGGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCgg                                                                                           <  1:329372/149‑1 (MQ=255)

                                                           aaTCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCg                                                                                            <  3:114005/148‑1 (MQ=255)

                                                            aTCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGt                                                                                          >  2:309036/1‑149 (MQ=255)

                                                               agagagagTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGttt                                                                                        >  4:769049/1‑148 (MQ=255)

                                                                   agagTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACtt                                                                                    <  3:943725/148‑1 (MQ=255)

                                                                     agTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTg                                                                                  <  3:61627/148‑1 (MQ=255)

                                                                     agTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGa                                                                                 <  1:728201/149‑1 (MQ=255)

                                                                     agTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGa                                                                                 <  2:323991/149‑1 (MQ=255)

                                                                      gTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGa                                                                                 <  3:328990/148‑1 (MQ=255)

                                                                      gTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGa                                                                                 <  3:715631/148‑1 (MQ=255)

                                                                         gcgACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGaa                                                                             >  2:821371/1‑149 (MQ=255)

                                                                               tttGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATaa                                                                        >  3:500062/1‑148 (MQ=255)

                                                                                   tgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                   >  2:1007953/1‑149 (MQ=255)

                                                                                   tgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                   >  1:210911/1‑149 (MQ=255)

                                                                                     aTGTTATTCTTACGACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCg                                                                  <  4:670466/148‑1 (MQ=255)

                                                                                     aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGt                                                                                          >  4:347386/1‑124 (MQ=255)

                                                                                     aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGt                                                                                          <  3:347386/124‑1 (MQ=255)

                                                                                     aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCg                                                                  <  4:241835/148‑1 (MQ=255)

                                                                                     aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCg                                                                  <  3:309507/148‑1 (MQ=255)

                                                                                          aTTCTTACAACTGTAAGAAGCGCTGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGttt                                                             <  4:677883/148‑1 (MQ=255)

                                                                                          aTTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGttt                                                             <  4:664169/148‑1 (MQ=255)

                                                                                           ttCTTACAACTGTAATAAGCGCCGGTGATTTACAATCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtat                                                           <  1:567259/149‑1 (MQ=255)

                                                                                           ttCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAaacataa                                                                                                                                 >  2:753060/1‑79 (MQ=255)

                                                                                            tCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtat                                                           <  4:500062/148‑1 (MQ=255)

                                                                                            tCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtat                                                           <  4:133053/148‑1 (MQ=255)

                                                                                              ttACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGc                                                        >  2:486608/1‑149 (MQ=255)

                                                                                               tACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGc                                                        >  4:174841/1‑148 (MQ=255)

                                                                                               tACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGc                                                        >  4:542829/1‑148 (MQ=255)

                                                                                                aCAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCAc                                                      >  1:955360/1‑149 (MQ=255)

                                                                                                  aaCTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAg                                                                                <  1:5866/121‑1 (MQ=255)

                                                                                                     tGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACcgcg                                                  >  4:589435/1‑148 (MQ=255)

                                                                                                      gTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGa                                                 >  3:25428/1‑148 (MQ=255)

                                                                                                        aagaagCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAgc                                               >  3:384952/1‑148 (MQ=255)

                                                                                                             gcgcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCaa                                          <  3:500356/148‑1 (MQ=255)

                                                                                                             gcgcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCaa                                          <  3:718977/148‑1 (MQ=255)

                                                                                                                  ggTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAaataat                                     <  4:384952/148‑1 (MQ=255)

                                                                                                                     gATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAAtttt                                  <  4:25428/148‑1 (MQ=255)

                                                                                                                         tACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGa                              <  3:769049/148‑1 (MQ=255)

                                                                                                                          aCAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGAcc                            <  2:477274/149‑1 (MQ=255)

                                                                                                                               caagGTTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACaa                        <  4:521596/148‑1 (MQ=255)

                                                                                                                               caagGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAc                       <  2:904660/149‑1 (MQ=255)

                                                                                                                               caagGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAc                       >  2:513879/1‑149 (MQ=255)

                                                                                                                               caagGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAc                       <  2:955360/149‑1 (MQ=255)

                                                                                                                                    gTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTc                  >  1:846543/1‑149 (MQ=255)

                                                                                                                                     tATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCg                 >  1:155532/1‑149 (MQ=255)

                                                                                                                                         ggtACGCTTTAAAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAata                                                                           <  2:1060957/85‑1 (MQ=255)

                                                                                                                                         ggtACGCTTTAAAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAata                                                                           >  1:1060957/3‑87 (MQ=255)

                                                                                                                                          gtACGCTTTAAAAAGCGTAATCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTAt              >  3:230481/2‑148 (MQ=255)

                                                                                                                                            aCGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACTTGTCGTATtgt           >  4:12687/1‑148 (MQ=255)

                                                                                                                                            aCGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg          <  2:13898/149‑1 (MQ=255)

                                                                                                                                            aCGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg          <  2:846543/149‑1 (MQ=255)

                                                                                                                                              gCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt        >  1:617700/1‑149 (MQ=255)

                                                                                                                                              gCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt        >  1:806398/1‑149 (MQ=255)

                                                                                                                                                tttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa      >  2:504468/1‑149 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCCTGAAACATAACATACTAAAAGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGATGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTATACGACCACAACCTGTCGTATTGTGGTCa      >  4:655042/1‑148 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGGGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg          >  4:619645/1‑144 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTTTCGTATtgtg          <  3:619645/144‑1 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGt                   >  4:1018873/1‑135 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGt                   <  3:1018873/135‑1 (MQ=255)

                                                                                                                                                 ttttAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa      >  4:793110/1‑148 (MQ=255)

                                                                                                                                                    tAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaaa  <  2:806295/149‑1 (MQ=255)

TTTCTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGT‑ACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAA  >  NC_000913/2493415‑2493711

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`