BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,406,298	A→G	Q689Q (CAA→CAG)	b4179 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,406,298	0	A	G	96.0%	526.0 / 7.4	175	Q689Q (CAA→CAG)	b4179	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/3); new base G (104/64); total (108/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406208‑4406387

cttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCg                                                                                         >  2:2310000/1‑93 (MQ=255)

cttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCg                                                                                         >  1:3948342/1‑93 (MQ=255)

     tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  2:833737/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                   <  2:4108780/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                   >  2:2942854/1‑93 (MQ=255)

        gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 >  1:3693004/1‑93 (MQ=255)

        gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 >  1:4401595/1‑93 (MQ=255)

         tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggg                                                                                >  2:2348266/1‑93 (MQ=255)

         tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggg                                                                                >  1:1799605/1‑93 (MQ=255)

         tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggg                                                                                <  2:547213/93‑1 (MQ=255)

          cTGGACTACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgggg                                                                               <  2:3224398/93‑1 (MQ=255)

            ggACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgggggg                                                                             >  1:4152713/1‑92 (MQ=255)

            ggACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGa                                                                             >  2:3243788/1‑93 (MQ=255)

               cgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATc                                                                          >  1:3006288/1‑93 (MQ=255)

                gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATcc                                                                         >  2:306906/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        >  1:4467516/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        >  2:667743/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        >  1:864870/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        >  1:3766832/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        >  1:1255955/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                        <  2:890779/93‑1 (MQ=255)

                   ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      >  2:543836/1‑93 (MQ=255)

                    tGTTCATTGATGGTCTGGTCCATGTCTCTTCTCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCg                                                                     >  2:4086655/1‑93 (MQ=255)

                    tGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCg                                                                     >  2:3897498/1‑93 (MQ=255)

                     gTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCgg                                                                    >  2:2815715/1‑93 (MQ=255)

                       tCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGcc                                                                  <  1:3485616/93‑1 (MQ=255)

                       tCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGcc                                                                  <  2:4211944/93‑1 (MQ=255)

                        cATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGAAGGGCAGCGACTGATGGGGGTATCCAGCGGCCa                                                                 >  2:3674243/1‑93 (MQ=255)

                         aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAg                                                                >  1:275287/1‑93 (MQ=255)

                         aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAg                                                                >  1:1076471/1‑93 (MQ=255)

                          ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               >  2:2254178/1‑93 (MQ=255)

                          ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               >  1:1653351/1‑93 (MQ=255)

                          ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               <  1:2133533/93‑1 (MQ=255)

                          ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               >  1:2158116/1‑93 (MQ=255)

                           tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGAc                                                              >  1:1407289/1‑93 (MQ=255)

                            gATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGCT‑GGGCAATCCAGCGGCCAGACtt                                                            >  1:1853364/1‑93 (MQ=255)

                            gATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACt                                                             >  1:3779994/1‑93 (MQ=255)

                             aTGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATTCAGCGGCCAGACtt                                                            >  2:1372749/1‑93 (MQ=255)

                               ggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTAt                                                          <  2:3369266/93‑1 (MQ=255)

                                gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                         <  1:1448451/93‑1 (MQ=255)

                                     gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    <  2:4266685/93‑1 (MQ=255)

                                     gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    <  1:2444807/93‑1 (MQ=255)

                                     gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    >  1:1595470/1‑93 (MQ=255)

                                       ccATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTggg                                                  >  2:1414501/1‑93 (MQ=255)

                                        cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTgggg                                                 >  1:387606/1‑92 (MQ=255)

                                        cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGc                                                 >  2:1219335/1‑93 (MQ=255)

                                        cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGc                                                 >  1:2006014/1‑93 (MQ=255)

                                        cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGc                                                 >  1:1763253/1‑93 (MQ=255)

                                         aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCg                                                >  2:3976212/1‑93 (MQ=255)

                                          tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                               >  1:255955/1‑93 (MQ=255)

                                          tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                               >  2:3297031/1‑93 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                              <  2:2034055/93‑1 (MQ=255)

                                            tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATc                                             <  1:3532367/93‑1 (MQ=255)

                                             ctctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcg                                            >  1:1580487/1‑93 (MQ=255)

                                              tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgc                                           >  2:1288053/1‑93 (MQ=255)

                                              tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgc                                           >  1:1009539/1‑93 (MQ=255)

                                              tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgc                                           >  2:3461460/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  2:2476703/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  1:988434/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         <  2:4460977/93‑1 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  1:345190/1‑93 (MQ=255)

                                                     tGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAg                                    >  1:1033895/1‑93 (MQ=255)

                                                     tGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGGGGAAg                                    >  1:603315/1‑93 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATTGCCTGGGCGATCGCGTGGCAGt                                   >  2:1294106/1‑93 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                   >  2:1738538/1‑93 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                   <  1:2331003/93‑1 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                   >  1:4040511/1‑93 (MQ=255)

                                                        aCAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTc                                 <  1:2867741/93‑1 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               >  2:111504/1‑93 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               >  2:1881374/1‑93 (MQ=255)

                                                            tGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGGAGTTCGCGt                             >  2:3341921/1‑93 (MQ=255)

                                                            tGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGt                             <  2:2555111/93‑1 (MQ=255)

                                                            tGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGACCGCGTGGAAGTTCGCGt                             >  1:4210935/1‑93 (MQ=255)

                                                             gACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTc                            >  2:2699769/1‑93 (MQ=255)

                                                               ctactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGa                          >  2:2170002/1‑93 (MQ=255)

                                                                tactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGaa                         <  1:239032/93‑1 (MQ=255)

                                                                 actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                        >  1:559763/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      >  1:722159/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      <  1:781473/93‑1 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      <  2:1625420/93‑1 (MQ=255)

                                                                    aTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                     >  1:1432626/1‑93 (MQ=255)

                                                                    aTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                     >  2:645199/1‑93 (MQ=255)

                                                                    aTCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGAATTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                     >  1:3088217/1‑93 (MQ=255)

                                                                       gCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTa                  >  2:2281006/1‑93 (MQ=255)

                                                                        cTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                 >  2:1466345/1‑93 (MQ=255)

                                                                          ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               <  2:865185/93‑1 (MQ=255)

                                                                            gACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATg             <  1:15325/93‑1 (MQ=255)

                                                                            gACCAGGTAGGGCAACACCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATg             <  1:4326226/93‑1 (MQ=255)

                                                                             aCCAGGTGGGGCAGCGCCTGATGGTGGAAGCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATgg            <  2:4376204/93‑1 (MQ=255)

                                                                             aCCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATgg            <  1:2263181/93‑1 (MQ=255)

                                                                               cAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGAc          >  1:4380139/1‑93 (MQ=255)

                                                                                aGGTAGGGCAGCGCCTGATGGGGTAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACg         <  1:1959169/93‑1 (MQ=255)

                                                                                aGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGGCTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACg         <  1:2876254/93‑1 (MQ=255)

                                                                                aGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACg         >  2:3176855/1‑93 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        <  1:2689425/93‑1 (MQ=255)

                                                                                  gTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg       >  2:4134371/1‑93 (MQ=255)

                                                                                  gTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg       >  2:4603659/1‑93 (MQ=255)

                                                                                  gTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg       <  2:3230084/93‑1 (MQ=255)

                                                                                   tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATTGACGAgc      <  2:3462553/93‑1 (MQ=255)

                                                                                     gggCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgcgc    >  1:3189855/1‑93 (MQ=255)

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406208‑4406387

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