BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,765	T→C	6.5%	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	6.5%	535.2 / 19.4	170	D338G (GAC→GGC)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand): ref base T (87/72); new base C (6/5); total (93/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 170 total aligned reads displayed.

TGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCAT  >  NC_000913/2581620‑2581912

 gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                                 <  3:88916/149‑1 (MQ=255)

     gTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCa                                                                                                                                             <  1:103436/149‑1 (MQ=255)

        aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                                                                          >  4:406255/1‑149 (MQ=255)

         aaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                                                                         >  4:98075/1‑149 (MQ=255)

           ggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATgttgt                                                                                                                                       <  4:185544/149‑5 (MQ=255)

             tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                     <  4:1009262/149‑1 (MQ=255)

             tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                     >  1:108695/1‑149 (MQ=255)

                cTGATCGTGATTGCCGGATGCTGATTGTCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                  <  4:1489498/149‑1 (MQ=255)

                   aTGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                               >  1:1672246/1‑149 (MQ=255)

                    tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                              >  1:789687/1‑149 (MQ=255)

                        gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                          <  4:1463512/149‑1 (MQ=255)

                         attgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTggg                                                                                                                         >  1:919562/1‑149 (MQ=255)

                          ttgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGt                                                                                                                        >  3:1240936/1‑149 (MQ=255)

                             ccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAATGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                  <  2:1492049/136‑1 (MQ=255)

                              cggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                  <  2:1566441/119‑1 (MQ=255)

                              cggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                  >  1:1566441/1‑119 (MQ=255)

                               ggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgc                                                                                                                   >  2:376012/1‑149 (MQ=255)

                                gatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                  <  4:1240936/149‑1 (MQ=255)

                                  tgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcggg                                                                                                                <  2:804637/149‑2 (MQ=255)

                                      gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                                                                                                             >  4:245121/1‑149 (MQ=255)

                                      gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                            >  2:1494908/1‑149 (MQ=255)

                                          gccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCa                                                                                                        <  2:1023062/149‑1 (MQ=255)

                                                   gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                    >  3:1512924/1‑144 (MQ=255)

                                                   gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                               >  2:301180/1‑149 (MQ=255)

                                                      gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                                 >  1:669104/1‑96 (MQ=255)

                                                      gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                                 <  2:669104/96‑1 (MQ=255)

                                                         gccgagTCTTATCCGGAGTAAAGCCTCTTAGCCGGCCTACAGGTCTCGCCA‑TTTTGTTTTACGGATGATGCACAACGTGGACCAGGCGTCTGGCCCATCCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGTACCTGGCGAAAACGCt                                                                                         >  2:1231144/6‑149 (MQ=255)

                                                         gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                         >  2:1613027/1‑149 (MQ=255)

                                                          aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                        >  3:697712/1‑149 (MQ=255)

                                                           cgcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCAAAAACGCTGa                                                                                       <  1:1494908/149‑1 (MQ=255)

                                                           cgcgTCTCATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                       <  3:764100/149‑1 (MQ=255)

                                                            gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                                                                 <  4:1143947/122‑1 (MQ=255)

                                                            gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                          <  2:1065517/145‑1 (MQ=255)

                                                            gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAc                                                                                      >  2:1186869/1‑149 (MQ=255)

                                                               tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCTCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                   <  4:841235/149‑1 (MQ=255)

                                                               tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                          <  2:531963/110‑1 (MQ=255)

                                                                 ttATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTc                                                                                                                                                    >  2:1319275/1‑82 (MQ=255)

                                                                   attcGGCCTACTGCGTCTTATCCGGCCTACGTGTTCCGCCA‑TTTTGTTTAACGGGTTATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGAAGGTGAGCTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                               <  4:1065307/146‑1 (MQ=255)

                                                                   aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                               >  1:829502/1‑149 (MQ=255)

                                                                    tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                         <  3:700843/138‑1 (MQ=255)

                                                                     ccGGCCTACAGCGTCTTATCCCGCCTACATGTCCCGCCA‑TTGTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgc                                                                             <  4:1624011/149‑1 (MQ=255)

                                                                      cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                            >  1:333375/1‑149 (MQ=255)

                                                                        gCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGcc                                                                          >  2:695610/1‑149 (MQ=255)

                                                                                cGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                                                                  <  2:341876/149‑1 (MQ=255)

                                                                                  tCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                                >  3:553883/1‑149 (MQ=255)

                                                                                   cTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACTGTCAAAGCGCCAGCGCGCCTGa                                                               >  3:1058689/1‑149 (MQ=255)

                                                                                   cTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGa                                                               >  3:1418530/1‑149 (MQ=255)

                                                                                    ttATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                              >  1:377426/1‑149 (MQ=255)

                                                                                       tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGGCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc                                                           <  2:1403127/149‑1 (MQ=255)

                                                                                       tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                     >  3:1579762/1‑107 (MQ=255)

                                                                                       tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                     <  4:1579762/107‑1 (MQ=255)

                                                                                       tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgt                                                                                                                                          >  3:1545444/1‑67 (MQ=37)

                                                                                       tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgt                                                                                                                                          <  4:1545444/70‑4 (MQ=37)

                                                                                         cGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCca                                                         >  4:1585560/1‑149 (MQ=255)

                                                                                             cTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAATGCGCCAGCGCGCCTGAATGCCACAAt                                                     <  4:697712/149‑1 (MQ=255)

                                                                                              tACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                    >  1:448612/1‑149 (MQ=255)

                                                                                               aCATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                   >  3:1249375/1‑149 (MQ=255)

                                                                                                cATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                  <  3:318541/149‑1 (MQ=255)

                                                                                                cATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                  <  3:246767/149‑1 (MQ=255)

                                                                                                      ccGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                            >  1:1187726/1‑149 (MQ=255)

                                                                                                       cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                           >  1:1296544/1‑149 (MQ=255)

                                                                                                       cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                           >  3:1020758/1‑149 (MQ=255)

                                                                                                       cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                           >  2:368897/1‑149 (MQ=255)

                                                                                                             ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                         <  3:1431033/98‑1 (MQ=255)

                                                                                                             ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                         >  4:1431033/1‑98 (MQ=255)

                                                                                                             ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                      >  3:1445157/1‑149 (MQ=255)

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

                                                                                                               ttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCt                                    <  4:287658/149‑1 (MQ=255)

                                                                                                                  tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTc                                 <  3:1613388/149‑1 (MQ=255)

                                                                                                                   ttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCg                                >  1:1565848/1‑149 (MQ=255)

                                                                                                                    tAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGc                               >  2:975202/1‑149 (MQ=255)

                                                                                                                     aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCa                                                                                                                                        >  4:1606252/1‑43 (MQ=255)

                                                                                                                     aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCt                              >  1:1408464/1‑149 (MQ=255)

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

                                                                                                                          gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGt                         <  4:931991/149‑1 (MQ=255)

                                                                                                                             aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                      >  3:457809/1‑149 (MQ=255)

                                                                                                                             aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                      <  1:1670416/149‑1 (MQ=255)

                                                                                                                             aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                      <  1:1186869/149‑1 (MQ=255)

                                                                                                                              tCCACAACGTGGGCCAGGCGTCTGGTCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAcc                     >  3:511474/1‑149 (MQ=255)

                                                                                                                              tCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAcc                     >  1:968815/1‑149 (MQ=255)

                                                                                                                              tCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAcc                     >  1:316364/1‑149 (MQ=255)

                                                                                                                                 acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                  >  2:860665/1‑149 (MQ=255)

                                                                                                                                    gcatGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCa               >  4:74425/4‑149 (MQ=255)

                                                                                                                                     catGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACACTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc              >  3:1655648/3‑149 (MQ=255)

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

                                                                                                                                     cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc              <  4:656132/149‑1 (MQ=255)

                                                                                                                                     cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc              >  4:1304686/1‑149 (MQ=255)

                                                                                                                                      gTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg             <  4:374469/149‑1 (MQ=255)

                                                                                                                                      atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg             >  4:504500/2‑149 (MQ=255)

                                                                                                                                        gggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                   >  4:826762/1‑77 (MQ=255)

                                                                                                                                         ggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                  <  1:53408/109‑1 (MQ=255)

                                                                                                                                         ggCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGa          <  1:1451720/149‑1 (MQ=255)

                                                                                                                                            cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGAttt       >  1:1601393/1‑149 (MQ=255)

                                                                                                                                             aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCgcgc                                       <  2:1684135/116‑1 (MQ=255)

                                                                                                                                             aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtc      >  2:288319/1‑149 (MQ=255)

                                                                                                                                               gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc              >  4:1366504/1‑139 (MQ=255)

                                                                                                                                               gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc    <  1:325554/149‑1 (MQ=255)

                                                                                                                                               gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc    >  1:948559/1‑149 (MQ=255)

                                                                                                                                               gCGTCTGGCCCAGGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc              <  3:1366504/139‑1 (MQ=255)

                                                                                                                                                 gTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCtcat  <  1:860665/149‑1 (MQ=255)

TGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTCAT  >  NC_000913/2581620‑2581912

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