BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,964,245	T→C	Y72H (TAC→CAC)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,964,245	0	T	C	100.0%	613.0 / NA	197	Y72H (TAC→CAC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (98/99); total (98/99)

Only 100 of 213 total aligned reads displayed.

CACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATT  >  NC_000913/3964155‑3964334

caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCAc                                                                                         <  1:2477913/93‑1 (MQ=255)

 aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCAcc                                                                                        >  1:1817680/1‑93 (MQ=255)

  cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCt                                                                                       <  2:372239/93‑1 (MQ=255)

  cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCt                                                                                       <  1:3746649/93‑1 (MQ=255)

  cGCAAAGAGTGGCGAAGATATCTTGGGTGATGGCGTACTGGAGATATTGCAGGATGGATTGGGTTTCCACCTTTCCGCAGACAGCTCCCACCt                                                                                       <  2:425126/93‑1 (MQ=255)

   gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTc                                                                                      <  1:1306825/93‑1 (MQ=255)

    cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCg                                                                                     >  1:1211013/1‑93 (MQ=255)

     aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGc                                                                                    <  2:345336/93‑1 (MQ=255)

      aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGcc                                                                                   <  2:3123101/93‑1 (MQ=255)

        gagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCgg                                                                                 >  2:3918949/1‑93 (MQ=255)

          gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTc                                                                               >  2:2059865/1‑93 (MQ=255)

           tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTcc                                                                              >  1:542669/1‑93 (MQ=255)

              cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtga                                                                           <  1:3444088/93‑1 (MQ=255)

               gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgat                                                                          >  2:3244874/1‑93 (MQ=255)

               gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgat                                                                          >  1:4676980/1‑93 (MQ=255)

                aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgatg                                                                         >  2:1228496/1‑93 (MQ=255)

                aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgatg                                                                         >  1:1028317/1‑93 (MQ=255)

                   atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACa                                                                      <  1:972769/93‑1 (MQ=255)

                    tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACAt                                                                     >  1:2667535/1‑93 (MQ=255)

                    tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACAt                                                                     <  2:2018661/93‑1 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  1:2662286/1‑93 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    <  2:1490090/93‑1 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  2:3844332/1‑93 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    <  1:3380194/93‑1 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  1:3411082/1‑93 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    <  1:4401359/93‑1 (MQ=255)

                     atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  1:887008/1‑93 (MQ=255)

                       cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATAACATCTa                                                                  <  1:2794359/93‑1 (MQ=255)

                        tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTAc                                                                 <  2:2184262/93‑1 (MQ=255)

                        tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTAc                                                                 >  2:817912/1‑93 (MQ=255)

                          tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGt                                                               >  2:76147/1‑93 (MQ=255)

                           ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGtt                                                              <  1:2162954/93‑1 (MQ=255)

                             tGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTc                                                            <  2:3582315/93‑1 (MQ=255)

                                  gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCGCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTa                                                       <  2:1213413/93‑1 (MQ=255)

                                     tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                    >  2:2823062/1‑93 (MQ=255)

                                      aCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCa                                                   >  2:4484103/1‑93 (MQ=255)

                                       cTGGAGATATTGCAGGATGGATTTGGTTTACTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaa                                                  >  2:4636333/1‑93 (MQ=255)

                                        tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaaa                                                 >  1:3092585/1‑93 (MQ=255)

                                        tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaaa                                                 >  1:1533195/1‑93 (MQ=255)

                                        tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaaa                                                 <  1:3543555/93‑1 (MQ=255)

                                         ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                <  1:588042/93‑1 (MQ=255)

                                         ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                >  1:612306/1‑93 (MQ=255)

                                          gagaTATTGCAGAATGGATTTGGGTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                               <  2:461720/93‑1 (MQ=255)

                                           agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATcc                                              >  2:3507281/1‑93 (MQ=255)

                                             atatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                                            <  2:100924/93‑1 (MQ=255)

                                             atatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                                            >  2:876243/1‑93 (MQ=255)

                                              tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcg                                           >  1:1468894/1‑92 (MQ=255)

                                              tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                                           <  2:514847/93‑1 (MQ=255)

                                               atTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgccg                                          >  1:565753/1‑93 (MQ=255)

                                                ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGt                                         <  2:1794193/93‑1 (MQ=255)

                                                 tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGtt                                        >  1:4302628/1‑93 (MQ=255)

                                                 tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGtt                                        <  2:153036/93‑1 (MQ=255)

                                                   cAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc                                      <  2:245350/93‑1 (MQ=255)

                                                    aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCa                                     >  1:1195456/1‑93 (MQ=255)

                                                    aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCa                                     >  2:2145212/1‑93 (MQ=255)

                                                    aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCa                                     >  2:4341951/1‑93 (MQ=255)

                                                    aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCCGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCa                                     >  2:3402498/1‑93 (MQ=255)

                                                     ggatggatTTGGTTTTCTTCGTTCCGCAGACAGCTGCCACCTCGGCGCTCCTGATGACATCTACGTTTGCCCTCGACAAATACCGCGGTTCaa                                    >  2:4413889/1‑93 (MQ=255)

                                                     ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCaa                                    >  2:2626816/1‑93 (MQ=255)

                                                     ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCaa                                    >  1:85916/1‑93 (MQ=255)

                                                      gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc                                   <  2:2157624/93‑1 (MQ=255)

                                                      gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc                                   >  1:2245761/1‑93 (MQ=255)

                                                       atggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                  >  1:4021050/1‑93 (MQ=255)

                                                        tggatTTGGTTTCCTCCGTTCCGCAGCCAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt                                 <  2:2329943/93‑1 (MQ=255)

                                                          gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                               >  2:3905353/1‑93 (MQ=255)

                                                          gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                               <  1:2233113/93‑1 (MQ=255)

                                                           atTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                              >  2:3705254/1‑93 (MQ=255)

                                                            tttGTTTTCCTCAGTTCCGTAGCCACCTCCCGCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGc                             <  2:4060003/93‑1 (MQ=255)

                                                             ttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                            <  2:1539090/93‑1 (MQ=255)

                                                             ttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                            <  1:1959626/93‑1 (MQ=255)

                                                             ttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                            <  2:2445370/93‑1 (MQ=255)

                                                             ttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCACCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                            >  2:580551/1‑93 (MQ=255)

                                                                gTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTg                         >  2:3097870/1‑93 (MQ=255)

                                                                  ttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGt                       <  2:1924508/93‑1 (MQ=255)

                                                                  ttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGt                       <  2:3259407/93‑1 (MQ=255)

                                                                    cctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     >  1:2010931/1‑93 (MQ=255)

                                                                    cctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     >  2:1038183/1‑93 (MQ=255)

                                                                    cctccGTCCCGCAGACAGCTCCCTCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     <  1:2401507/93‑1 (MQ=255)

                                                                     ctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGAt                    >  2:3774189/1‑93 (MQ=255)

                                                                     ctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGAt                    <  1:2315789/93‑1 (MQ=255)

                                                                      tccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATa                   <  2:2748971/93‑1 (MQ=255)

                                                                      tccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATa                   >  2:1374664/1‑93 (MQ=255)

                                                                       ccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAc                  <  2:490907/93‑1 (MQ=255)

                                                                         gTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                <  2:33426/93‑1 (MQ=255)

                                                                         gTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                <  1:4506719/93‑1 (MQ=255)

                                                                         gTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                >  1:3032302/1‑93 (MQ=255)

                                                                          ttCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACCTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt               >  2:600083/1‑93 (MQ=255)

                                                                          ttCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt               >  2:4395994/1‑93 (MQ=255)

                                                                           tCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCAtc              <  2:4187660/93‑1 (MQ=255)

                                                                             cGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtctc            <  1:1434047/93‑1 (MQ=255)

                                                                               cAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTg          <  2:1883017/93‑1 (MQ=255)

                                                                                aGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg         <  1:1622347/93‑1 (MQ=255)

                                                                                 gACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt        <  2:2837194/93‑1 (MQ=255)

                                                                                  acacCTCCCACCCCGCCGTTCCTGATGCCATCTACGTTCCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa       <  1:3218012/93‑1 (MQ=255)

                                                                                  aCAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa       <  1:3844284/93‑1 (MQ=255)

                                                                                  aCAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa       >  2:1272663/1‑93 (MQ=255)

                                                                                   cAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa      <  1:72344/93‑1 (MQ=255)

                                                                                   cAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa      <  1:936147/93‑1 (MQ=255)

                                                                                       tCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt  <  1:2951830/93‑1 (MQ=255)

                                                                                       tCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGGTTCAACCTCCGCACTGGTGATCCCATCTCTGGTAAGAtt  >  2:1053360/1‑93 (MQ=255)

CACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATT  >  NC_000913/3964155‑3964334

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`