BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	97.0%	517.2 / 7.9	166	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (100/60); total (105/61)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 178 total aligned reads displayed.

GCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006389‑4006567

gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c                                                                                        <  2:2069229/93‑1 (MQ=255)

    gcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cccca                                                                                    >  1:2069778/1‑91 (MQ=255)

     ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ccccaa                                                                                   <  2:4125101/93‑4 (MQ=255)

      cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCACCACTTCAAACGCAATAACTGAGCGGGAAAATAGCG‑T‑ccccaac                                                                                  >  1:3517937/1‑89 (MQ=255)

         tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAAcccc                                                                                >  2:1579363/1‑93 (MQ=255)

         tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAAcccc                                                                                >  1:992512/1‑93 (MQ=255)

         tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAAcccc                                                                                <  1:2560446/93‑1 (MQ=255)

          gctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCa                                                                               >  2:2877703/1‑93 (MQ=255)

           ctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCAc                                                                              >  1:4111504/1‑93 (MQ=255)

              gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGcc                                                                           >  1:4593937/1‑93 (MQ=255)

               aTTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCa                                                                          >  2:4504940/1‑93 (MQ=255)

                ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCac                                                                         >  1:491128/1‑92 (MQ=255)

                ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCaa                                                                         <  1:1222161/93‑1 (MQ=255)

                 tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAAc                                                                        >  2:876206/1‑93 (MQ=255)

                 tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAAc                                                                        <  2:171338/93‑1 (MQ=255)

                 tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAAc                                                                        <  2:3862284/93‑1 (MQ=255)

                  cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAAcc                                                                       <  2:3446297/93‑1 (MQ=255)

                   ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCa                                                                      >  1:1121031/1‑93 (MQ=255)

                    cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAg                                                                     >  1:3611798/1‑93 (MQ=255)

                     aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGc                                                                    >  1:1637878/1‑93 (MQ=255)

                         aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACaa                                                                >  2:4323062/1‑93 (MQ=255)

                          aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAat                                                               >  1:4287548/1‑93 (MQ=255)

                            tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAatat                                                             <  1:146405/93‑1 (MQ=255)

                             cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATg                                                            <  1:775487/93‑1 (MQ=255)

                               aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGaa                                                          >  1:2168516/1‑93 (MQ=255)

                                gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  1:681323/1‑93 (MQ=255)

                                gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  2:2452733/93‑1 (MQ=255)

                                gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  2:465984/93‑1 (MQ=255)

                                gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  1:1408088/93‑1 (MQ=255)

                                  aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       >  1:3051936/1‑93 (MQ=255)

                                  aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       >  2:4571570/1‑93 (MQ=255)

                                  aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       >  1:504843/1‑93 (MQ=255)

                                    gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  2:2989136/1‑93 (MQ=255)

                                      cccgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAg                                                   <  1:3701145/91‑1 (MQ=255)

                                       gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  2:1125263/1‑93 (MQ=255)

                                        cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:1019604/93‑1 (MQ=255)

                                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  2:595850/1‑93 (MQ=255)

                                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  2:652610/1‑93 (MQ=255)

                                          ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                                               <  2:1366189/93‑1 (MQ=255)

                                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:3801753/1‑93 (MQ=255)

                                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:1792799/1‑93 (MQ=255)

                                            tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:4629378/1‑93 (MQ=255)

                                            cgCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             <  1:863349/92‑1 (MQ=255)

                                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  1:4474514/93‑1 (MQ=255)

                                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  1:1937126/93‑1 (MQ=255)

                                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGGATCGCCAGTCCGGTCTg                                          <  1:1002977/93‑1 (MQ=255)

                                                cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  2:3469534/1‑93 (MQ=255)

                                                cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  1:4471828/93‑1 (MQ=255)

                                                cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  2:1837278/1‑93 (MQ=255)

                                                  cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  1:4487468/1‑93 (MQ=255)

                                                   ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:3438479/1‑93 (MQ=255)

                                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  1:932116/93‑1 (MQ=255)

                                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:2768378/1‑93 (MQ=255)

                                                       acacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  2:408727/1‑93 (MQ=255)

                                                         acttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGCAGCCCAGCa                                >  1:4547213/1‑93 (MQ=255)

                                                         acttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:3885014/1‑93 (MQ=255)

                                                         acttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:4163999/1‑93 (MQ=255)

                                                          cttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:2372056/1‑93 (MQ=255)

                                                           ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  2:3104688/1‑93 (MQ=255)

                                                           ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  1:1816676/93‑1 (MQ=255)

                                                           ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  1:1804662/1‑93 (MQ=255)

                                                            tcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  2:1153965/93‑1 (MQ=255)

                                                            tcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  1:1598251/93‑1 (MQ=255)

                                                             caaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCGGACCGGTATGAAGCCCAGAAAttg                            >  1:4045595/1‑91 (MQ=255)

                                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  2:4695893/1‑93 (MQ=255)

                                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:1551765/93‑1 (MQ=255)

                                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  2:1270604/1‑93 (MQ=255)

                                                               aaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  2:1325350/93‑1 (MQ=255)

                                                               aaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:3488599/1‑93 (MQ=255)

                                                               aaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  1:4072304/93‑1 (MQ=255)

                                                                 cGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        <  1:2476018/93‑1 (MQ=255)

                                                                  gCAATCACTGAGCGGGAAAATAGCGTT‑CCCAAACCCACGCCAACCAGCACAATATGAATCGGCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  2:2660056/1‑93 (MQ=255)

                                                                  gCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  1:3100527/1‑93 (MQ=255)

                                                                  gCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  1:2017191/1‑93 (MQ=255)

                                                                    aaTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:1057981/1‑93 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  2:2702052/1‑93 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  1:1280982/93‑1 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:2975237/1‑93 (MQ=255)

                                                                      tCACTGAGCGGGAAAATATCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  1:1542892/1‑93 (MQ=255)

                                                                        aCTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:4037339/1‑93 (MQ=255)

                                                                         cTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                >  2:2180224/1‑93 (MQ=255)

                                                                         cTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                <  2:1911133/93‑1 (MQ=255)

                                                                           gAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              <  1:904488/93‑1 (MQ=255)

                                                                             gCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  2:338909/1‑93 (MQ=255)

                                                                             gCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  1:3781383/1‑93 (MQ=255)

                                                                               gggAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  2:920390/1‑93 (MQ=255)

                                                                                 gAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAACAGACGCCACCGCGCCGCGATaa        <  2:4172342/93‑1 (MQ=255)

                                                                                 gAAAATAGCG‑TCCCCAACCACACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  2:2369524/93‑1 (MQ=255)

                                                                                  aaaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  2:3694728/1‑93 (MQ=255)

                                                                                  aaaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  2:680814/1‑93 (MQ=255)

                                                                                  aaaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       <  1:3374629/93‑1 (MQ=255)

                                                                                  aaaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       <  1:3712990/93‑1 (MQ=255)

                                                                                   aaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      <  2:2980262/93‑1 (MQ=255)

                                                                                   aaaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      >  2:3715171/1‑93 (MQ=255)

                                                                                    aaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  1:4448673/1‑93 (MQ=255)

                                                                                    aaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  1:2866552/1‑93 (MQ=255)

                                                                                      tAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  1:2302345/93‑1 (MQ=255)

                                                                                      tAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCAATAACCGTg   <  2:2148809/93‑1 (MQ=255)

                                                                                       aGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  >  1:2957505/1‑93 (MQ=255)

                                                                                       aGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  >  2:4151686/1‑93 (MQ=255)

GCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006389‑4006567

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