BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,765	T→C	5.3%	D338G (GAC→GGC)	ypfG ←	DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,765	0	T	C	5.3%	365.2 / 7.7	114	D338G (GAC→GGC)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand): ref base T (59/49); new base C (4/2); total (63/51)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 114 total aligned reads displayed.

CTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  NC_000913/2581619‑2581910

cTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGACCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCt                                                                                                                                                 <  2:607441/149‑1 (MQ=255)

 tGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCt                                                                                                                                                 <  3:577128/148‑1 (MQ=255)

  gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                >  4:340128/1‑148 (MQ=255)

  gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                >  3:780396/1‑148 (MQ=255)

   ggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  4:146008/1‑148 (MQ=255)

    gAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGc                                                                                                                                              <  4:559672/148‑1 (MQ=255)

     aGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGcc                                                                                                                                             <  4:868472/148‑1 (MQ=255)

         aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATg                                                                                                                                         <  4:455760/148‑1 (MQ=255)

          aaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgt                                                                                                                                        <  4:272516/148‑4 (MQ=255)

                 cTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGc                                                                                                                                              >  2:986783/1‑135 (MQ=255)

                 cTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                <  2:389755/149‑1 (MQ=255)

                  tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                >  3:297633/1‑148 (MQ=255)

                   gATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                               <  4:918011/148‑1 (MQ=255)

                     tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                             <  3:463622/148‑1 (MQ=255)

                     tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                            <  2:115881/149‑1 (MQ=255)

                      ggTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                            <  4:500348/148‑1 (MQ=255)

                       gTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAgg                                                                                                                           >  4:352798/1‑148 (MQ=255)

                        tgattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGt                                                                                                                          >  3:215325/1‑148 (MQ=255)

                        tgattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGt                                                                                                                          <  4:402153/148‑1 (MQ=255)

                         gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        <  1:514396/149‑1 (MQ=255)

                         gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        <  2:751393/149‑1 (MQ=255)

                         gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        <  1:314706/149‑1 (MQ=255)

                         gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTg                                                                                                                         <  4:525581/148‑1 (MQ=255)

                                ggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTg                                                                                                                  <  4:123375/148‑1 (MQ=255)

                                      tgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                                                                                                           >  1:553028/1‑149 (MQ=255)

                                      tgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGAAGGTGGGTTCTGCGGCGTa                                                                                                           >  1:1047296/1‑149 (MQ=255)

                                       gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                             <  2:95360/130‑1 (MQ=255)

                                       gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                             >  1:95360/1‑130 (MQ=255)

                                       gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                          >  1:157043/1‑149 (MQ=255)

                                            ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCAc                                                                                                     <  1:292006/149‑1 (MQ=255)

                                            ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGGTCTGCGGCGTAACGCa                                                                                                      >  3:938526/1‑148 (MQ=255)

                                              ggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                   >  2:185904/1‑149 (MQ=255)

                                               gatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCGAGGCGTCTgg                                                                                                                                               <  3:601256/104‑1 (MQ=255)

                                               gatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  4:601256/1‑104 (MQ=255)

                                                   cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCg                                                                                               >  3:634225/1‑148 (MQ=255)

                                                    gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                             >  1:850753/1‑149 (MQ=255)

                                                     aCGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                             >  3:99218/1‑148 (MQ=255)

                                                      cGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaaaa                                                                                           <  1:493554/149‑1 (MQ=255)

                                                         tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACg                                                                                         <  3:801440/148‑1 (MQ=255)

                                                         tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                        <  2:1038009/149‑1 (MQ=255)

                                                         tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                        >  1:419739/1‑149 (MQ=255)

                                                            cgcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                      <  4:900406/148‑1 (MQ=255)

                                                              cgTCTTATCCGGCCTACAGCGTCTTATCCGGCCCACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAcc                                                                                   <  1:883025/149‑1 (MQ=255)

                                                                tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                 >  1:534542/1‑149 (MQ=255)

                                                                 cTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                 >  3:895412/1‑148 (MQ=255)

                                                                    aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                             >  2:28213/1‑149 (MQ=255)

                                                                    aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGACCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                             >  1:597412/1‑149 (MQ=255)

                                                                      ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                            >  4:188466/1‑148 (MQ=255)

                                                                        ggCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcgc                                                                         >  1:435893/1‑149 (MQ=255)

                                                                         gCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGcc                                                                        >  1:870410/1‑149 (MQ=255)

                                                                            tACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATg                                                                                                                                         >  3:328473/1‑81 (MQ=255)

                                                                            tACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATg                                                                                                                                         <  4:328473/81‑1 (MQ=255)

                                                                            tACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgc                                                                     >  2:426035/1‑149 (MQ=255)

                                                                                gCGTCTTATCCGGGCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcgcg                                                                  <  4:830738/148‑1 (MQ=255)

                                                                                 cGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcgcgc                                                                 <  3:973569/148‑1 (MQ=255)

                                                                                  gTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                                                                <  4:780396/148‑1 (MQ=255)

                                                                                   tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  2:913166/1‑149 (MQ=255)

                                                                                   tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  1:960324/1‑149 (MQ=255)

                                                                                      tATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGTCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                            <  4:598352/148‑1 (MQ=255)

                                                                                      tATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                            >  3:896728/1‑148 (MQ=255)

                                                                                      tATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                            <  3:340128/148‑1 (MQ=255)

                                                                                       aTCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                           <  4:1031204/148‑1 (MQ=255)

                                                                                       aTCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCCGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATg                                                          >  1:187934/1‑149 (MQ=255)

                                                                                        tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc                                                         >  2:982550/1‑149 (MQ=255)

                                                                                              cTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACaa                                                    <  4:777456/148‑1 (MQ=255)

                                                                                                aCATGTCCCGCCGTTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGGCAGTAATCGCTGGAGGGTTCTGCAGCGTAACGCACCAGGCGAAAACGCTGACC‑TCAAAGCGCCAGCGCCCCTTAATGCCACAATCGc                                                >  2:401562/1‑149 (MQ=255)

                                                                                                aCATGTCCCGCCATTTTGTTTAACGGGTGATCCACAAAGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                 >  2:76986/1‑149 (MQ=255)

                                                                                                      cccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTaa                                            >  4:59664/1‑148 (MQ=255)

                                                                                                      cccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGATGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTaa                                            >  3:904471/1‑148 (MQ=255)

                                                                                                       ccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                          >  2:506172/1‑149 (MQ=255)

                                                                                                        cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  2:776013/1‑149 (MQ=255)

                                                                                                            aTTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCgcgc                                     <  2:850753/149‑1 (MQ=255)

                                                                                                             ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                    >  2:265111/1‑149 (MQ=255)

                                                                                                                 gTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCttt                                >  1:938284/1‑149 (MQ=255)

                                                                                                                  tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTACTCCGCGCCCttt                                >  3:489414/1‑148 (MQ=255)

                                                                                                                       cGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                          <  2:68131/149‑1 (MQ=255)

                                                                                                                       cGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                          <  2:246353/149‑1 (MQ=255)

                                                                                                                          gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAgg                        <  3:146008/148‑1 (MQ=255)

                                                                                                                           tGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTc                      >  1:607582/1‑149 (MQ=255)

                                                                                                                             aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                    >  2:818310/1‑149 (MQ=255)

                                                                                                                               ccACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCa                  <  2:477232/149‑1 (MQ=255)

                                                                                                                                cacaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAg                 >  2:86818/1‑149 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                >  4:265958/1‑53 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                <  3:265958/53‑1 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                           <  2:234131/106‑1 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAAt                                                           >  1:234131/1‑106 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                <  4:942474/117‑1 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                >  3:942474/1‑117 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAg                 >  4:173184/1‑148 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAg                 >  3:989101/1‑148 (MQ=255)

                                                                                                                                  caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                >  3:303943/1‑148 (MQ=255)

                                                                                                                                  caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGGCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                >  3:933059/1‑148 (MQ=255)

                                                                                                                                   aaCGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTc              >  1:134674/1‑149 (MQ=255)

                                                                                                                                     cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc            <  2:419739/149‑1 (MQ=255)

                                                                                                                                      atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           >  1:199094/2‑149 (MQ=255)

                                                                                                                                       tGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           <  3:848354/148‑1 (MQ=255)

                                                                                                                                       tGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTTATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           >  3:438154/1‑148 (MQ=255)

                                                                                                                                       tGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           >  3:768652/1‑148 (MQ=255)

                                                                                                                                        gggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcgcg         >  2:887528/1‑149 (MQ=255)

                                                                                                                                          gCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGAt       >  2:902849/1‑149 (MQ=255)

CTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  NC_000913/2581619‑2581910

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