BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,793,668	C→A	6.2%	E148* (GAA→TAA)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,668	0	C	A	6.2%	436.2 / 12.5	125	E148* (GAA→TAA)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand): ref base C (61/56); new base A (2/6); total (63/62)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

AGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAAC  >  NC_000913/1793520‑1793815

aGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTc                                                                                                                                                     <  2:426420/149‑1 (MQ=255)

 gATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCa                                                                                                                                                    >  2:834666/1‑149 (MQ=255)

  aTTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaa                                                                                                                                                   <  4:391532/149‑1 (MQ=255)

         aaCATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcct                                                                                                                                            <  2:595394/149‑2 (MQ=255)

          aCATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctg                                                                                                                                           <  1:730096/149‑3 (MQ=255)

          aCATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctg                                                                                                                                           <  2:683195/149‑3 (MQ=255)

            aTCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctggt                                                                                                                                         <  1:1150735/149‑5 (MQ=255)

              cTTGTAACTGCATCACAATAGACTTTATTCTGCCAACGCCAGTTTAATGCTTTCAATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGCCCCTGCCAACCAGGAATTTTTCAAAATTCCATaaaa                                                                                                                                       <  1:1226415/149‑1 (MQ=255)

                tGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAata                                                                                                                                     <  3:792678/149‑1 (MQ=255)

                     cTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcatc                                                                                                                                >  4:18227/1‑149 (MQ=255)

                     cTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcatc                                                                                                                                <  3:1158400/149‑1 (MQ=255)

                     cTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcatc                                                                                                                                <  1:917278/149‑1 (MQ=255)

                       ccacaacaAAAAACATTATTTTGCAACCCCCGGTTTAAGCCTTTCCACCAC‑ATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCgg                                                                                                                             <  1:572678/144‑1 (MQ=255)

                         aTCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCggg                                                                                                                            >  1:701887/1‑149 (MQ=255)

                              aaTAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAg                                                                                                                       >  1:1516779/1‑149 (MQ=255)

                               aTAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                                                                                                                      >  4:43473/1‑149 (MQ=255)

                                  gACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGt                                                                                                                   <  1:241339/149‑1 (MQ=255)

                                  gACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGt                                                                                                                   <  4:1157582/149‑1 (MQ=255)

                                  gACATCATTTTGCCAACGCCAGTTTAACGCTTTCCCTCACAATGGGATCTTCCGGCGTAATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATTCGGGTACAGCGGt                                                                                                                   <  4:1357714/149‑1 (MQ=255)

                                     attattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCa                                                                                                                >  3:1338582/1‑149 (MQ=255)

                                     attattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCa                                                                                                                >  2:869241/1‑149 (MQ=255)

                                       tattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACg                                                                                                              >  1:359324/1‑149 (MQ=255)

                                         ttttGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGc                                                                                                            >  2:1192894/1‑149 (MQ=255)

                                           ttGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCt                                                                                                          <  2:815954/149‑1 (MQ=255)

                                           ttGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCt                                                                                                          <  1:1467197/149‑1 (MQ=255)

                                             gCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCttt                                                                                                        <  2:1516779/149‑1 (MQ=255)

                                              ccAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       <  3:749395/149‑1 (MQ=255)

                                              ccAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       <  1:1192894/149‑1 (MQ=255)

                                               cAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGc                                                                                                      >  1:1344637/1‑149 (MQ=255)

                                                aaCGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                     >  3:1262042/1‑149 (MQ=255)

                                                aaCGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                     <  1:430127/149‑1 (MQ=255)

                                                aaCGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                     >  1:930424/1‑149 (MQ=255)

                                                aaCGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                     <  2:701887/149‑1 (MQ=255)

                                                   gCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGt                                                                                                                   >  2:1205451/1‑132 (MQ=255)

                                                   gCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGAc                                                                                                  >  4:189811/1‑149 (MQ=255)

                                                     cAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCa                                                                                                >  2:420679/1‑149 (MQ=255)

                                                        tttAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGa                                                                                                   >  3:618859/1‑143 (MQ=255)

                                                        tttAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGa                                                                                                   <  4:618859/143‑1 (MQ=255)

                                                        tttAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATAc                                                                                             <  3:313117/149‑1 (MQ=255)

                                                              gCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                       <  3:628939/149‑1 (MQ=255)

                                                                   ccATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAAt                                                                                  <  2:225265/149‑1 (MQ=255)

                                                                    cATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATc                                                                                 <  2:1498615/149‑1 (MQ=255)

                                                                         caATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctc                                                                            <  1:1289845/149‑1 (MQ=255)

                                                                         caATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctc                                                                            <  3:43473/149‑1 (MQ=255)

                                                                          aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctct                                                                           >  1:1501893/1‑149 (MQ=255)

                                                                          aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctct                                                                           >  2:455304/1‑149 (MQ=255)

                                                                          aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctct                                                                           <  4:1211893/149‑1 (MQ=255)

                                                                          aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctgg                                                                                                                                          >  3:1100321/1‑83 (MQ=255)

                                                                          aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctgg                                                                                                                                          <  4:1100321/86‑4 (MQ=255)

                                                                            tGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGAcc                                                                                                 <  2:1089142/125‑1 (MQ=255)

                                                                            tGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGAcc                                                                                                 >  1:1089142/1‑125 (MQ=255)

                                                                            tGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctctct                                                                         >  4:158290/1‑149 (MQ=255)

                                                                             gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       >  1:939297/1‑118 (MQ=255)

                                                                             gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       <  2:939297/118‑1 (MQ=255)

                                                                                aTCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCg                                                                     >  3:173633/1‑149 (MQ=255)

                                                                                aTCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCg                                                                     <  1:206061/149‑1 (MQ=255)

                                                                                aTCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCg                                                                     >  2:983182/1‑149 (MQ=255)

                                                                                aTCTTCCGGCGTCATATCCGGGGAAAAACGCCGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCg                                                                     >  4:113470/1‑149 (MQ=255)

                                                                                   ttCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcg                                                                  >  1:51547/1‑149 (MQ=255)

                                                                                        gCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACt                                                             <  1:566825/149‑1 (MQ=255)

                                                                                                  cGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGgcg                                                   >  4:1346789/1‑149 (MQ=255)

                                                                                                         aaaCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGc                                            <  4:579770/149‑1 (MQ=255)

                                                                                                           aCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaa                                          <  2:572887/149‑1 (MQ=255)

                                                                                                           aCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaa                                          >  1:1574133/1‑149 (MQ=255)

                                                                                                            cGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaat                                         <  2:1344637/149‑1 (MQ=255)

                                                                                                               tGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaa                                      >  4:581240/1‑149 (MQ=255)

                                                                                                                 gATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       >  4:527717/1‑82 (MQ=255)

                                                                                                                 gATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTg                                                                                                       <  3:527717/82‑1 (MQ=255)

                                                                                                                   tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAtttt                                  >  3:253331/1‑149 (MQ=255)

                                                                                                                   tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGAATGCGGTCGGCGCTGCGGCAATCAAtttt                                  <  3:1346789/149‑1 (MQ=255)

                                                                                                                    gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                 <  2:930424/149‑1 (MQ=255)

                                                                                                                    gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                 <  2:661184/149‑1 (MQ=255)

                                                                                                                    gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                 <  1:412015/149‑1 (MQ=255)

                                                                                                                      cTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGa                               <  2:603931/149‑1 (MQ=255)

                                                                                                                      cTTTTCCGTCCCTGCCAACCAGGAATTTTGCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGa                               <  3:231461/149‑1 (MQ=255)

                                                                                                                       ttttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGAt                              >  4:1383339/1‑149 (MQ=255)

                                                                                                                       ttttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGAt                              >  3:657766/1‑149 (MQ=255)

                                                                                                                        tttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATa                             <  2:1175940/149‑1 (MQ=255)

                                                                                                                         ttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATAc                            >  3:585710/1‑149 (MQ=255)

                                                                                                                             gTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCg                        >  3:1475222/1‑149 (MQ=255)

                                                                                                                              tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCgg                       >  4:364093/1‑149 (MQ=255)

                                                                                                                              tCCCTGCCAACCAGGAATTTTTAAAAATTCAACAAAATATCATCCGTGCACAGCAGAGCACGGCCTAAGCTGACCATACGGTTATACAATCAGCAATATGCCGGCCCGATTCCGTTCGCCGCTGCGGCAATCAATTTTTGATAcctcgg                       >  4:272935/1‑144 (MQ=255)

                                                                                                                               cccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTCCTGACCATACGGGCATAGGCTCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACCGTGGa                      <  4:1338582/149‑1 (MQ=255)

                                                                                                                                  tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATgg                   >  1:203830/1‑149 (MQ=255)

                                                                                                                                  tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATgg                   >  2:1424094/1‑149 (MQ=255)

                                                                                                                                   gCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGc                  <  4:1475222/149‑1 (MQ=255)

                                                                                                                                   gCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGc                  >  2:538857/1‑149 (MQ=255)

                                                                                                                                          agaaaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           <  1:1019820/146‑1 (MQ=255)

                                                                                                                                          aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGAGCATACGGGCATAGAATCCGCTCTCTTCCGGCGAGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           >  4:1496471/1‑149 (MQ=255)

                                                                                                                                          aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           >  3:701940/1‑149 (MQ=255)

                                                                                                                                          aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           >  1:406890/1‑149 (MQ=255)

                                                                                                                                          aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           >  1:715098/1‑149 (MQ=255)

                                                                                                                                          aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc           >  4:239845/1‑149 (MQ=255)

                                                                                                                                           ggAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATAAATTTTTGATACAGCGGATGGCGTCCTTCg          >  3:826350/1‑149 (MQ=255)

                                                                                                                                           ggAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTACCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGCTACAGCGGATGGCGTCCTTCg          >  1:196793/1‑149 (MQ=255)

                                                                                                                                             aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGcc        >  2:783028/1‑149 (MQ=255)

                                                                                                                                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTaa   >  1:122498/1‑149 (MQ=255)

                                                                                                                                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTaa   <  1:513431/149‑1 (MQ=255)

                                                                                                                                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTaa   >  4:1524980/1‑149 (MQ=255)

                                                                                                                                                   tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAAc  <  2:825119/149‑1 (MQ=255)

AGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAAC  >  NC_000913/1793520‑1793815

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