BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,500,011	G→A	R12Q (CGG→CAG)	b3372 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,500,011	0	G	A	82.4%	731.2 / 108.6	284	R12Q (CGG→CAG)	b3372	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (34/16); new base A (124/110); total (158/126)
Fisher's exact test for biased strand distribution p-value = 6.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 294 total aligned reads displayed.

ATATTAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTCGCGCCGG  >  NC_000913/3499921‑3500103

atatTAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCagcagc                                                                                            <  2:5073981/93‑1 (MQ=255)

 tatTAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCt                                                                                           >  2:1645008/1‑93 (MQ=255)

  atTAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTg                                                                                          >  1:646072/1‑93 (MQ=255)

   ttAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGc                                                                                         <  1:1052523/93‑1 (MQ=255)

   ttAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGc                                                                                         <  1:5336635/93‑1 (MQ=255)

    tAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGcc                                                                                        >  1:7279137/1‑93 (MQ=255)

    tAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGcc                                                                                        <  1:4363124/93‑1 (MQ=255)

       tCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTAt                                                                                     >  1:7358621/1‑93 (MQ=255)

       tCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTAt                                                                                     <  1:2022972/93‑1 (MQ=255)

       tCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTAt                                                                                     >  1:6330867/1‑93 (MQ=255)

          ctctCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGa                                                                                  >  2:1002131/1‑93 (MQ=255)

          ctctCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGa                                                                                  >  1:6019962/1‑93 (MQ=255)

          ctctCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGa                                                                                  <  2:3412599/93‑1 (MQ=255)

           tctcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGaa                                                                                 >  2:937608/1‑93 (MQ=255)

           tctcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGaa                                                                                 <  2:2448663/93‑1 (MQ=255)

           tctcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGaa                                                                                 >  2:5997738/1‑93 (MQ=255)

           tctcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGaa                                                                                 >  1:3902299/1‑93 (MQ=255)

             tcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAaca                                                                               >  1:4578817/1‑93 (MQ=255)

             tcCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAaca                                                                               >  2:1754282/1‑93 (MQ=255)

              cccggcccggCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAacac                                                                              >  1:1622588/1‑93 (MQ=255)

               ccggcccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACg                                                                             >  2:5821966/1‑93 (MQ=255)

               ccggcccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACg                                                                             >  1:4972150/1‑93 (MQ=255)

                cggcccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGc                                                                            >  1:2962771/1‑93 (MQ=255)

                cggcccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGc                                                                            >  2:2707709/1‑93 (MQ=255)

                 ggcccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCa                                                                           <  1:6135640/93‑1 (MQ=255)

                   cccggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCAtt                                                                         <  2:1217368/93‑1 (MQ=255)

                     cggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTc                                                                       <  1:1388808/93‑1 (MQ=255)

                     cggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTc                                                                       >  2:2956838/1‑93 (MQ=255)

                     cggcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTc                                                                       >  1:6738642/1‑93 (MQ=255)

                       gcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGt                                                                     >  2:3747737/1‑93 (MQ=255)

                       gcAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGt                                                                     <  1:6690669/93‑1 (MQ=255)

                         aaCGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGa                                                                   >  2:6122791/1‑93 (MQ=255)

                         aaCGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGa                                                                   <  1:5480049/93‑1 (MQ=255)

                            gggCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGc                                                                <  1:3708294/93‑1 (MQ=255)

                             ggCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCa                                                               <  1:6408530/93‑1 (MQ=255)

                              gCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCaa                                                              >  2:1832488/1‑93 (MQ=255)

                              gCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCaa                                                              <  2:7192256/93‑1 (MQ=255)

                               ccGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAg                                                             >  1:3788902/1‑93 (MQ=255)

                               ccGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAg                                                             <  1:6188583/93‑1 (MQ=255)

                               ccGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAg                                                             <  1:5009130/93‑1 (MQ=255)

                               ccGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAg                                                             <  1:4755926/93‑1 (MQ=255)

                               ccGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAg                                                             <  1:2666623/93‑1 (MQ=255)

                                  gATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGa                                                          <  1:3403797/93‑1 (MQ=255)

                                    tttACGCAAGGAGTTACCCGATGCAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGc                                                        >  2:171194/1‑93 (MQ=255)

                                    tttACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGc                                                        <  2:6216221/93‑1 (MQ=255)

                                    tttACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGc                                                        >  1:6267833/1‑93 (MQ=255)

                                     ttACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCAAGCt                                                       <  2:7374193/93‑1 (MQ=255)

                                      tACGCAAGGAGTTCCCCGAGGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTg                                                      <  1:6929298/93‑1 (MQ=255)

                                      tACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTg                                                      <  2:5344948/93‑1 (MQ=255)

                                       aCGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTgg                                                     >  2:2055606/1‑93 (MQ=255)

                                       aCGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTgg                                                     <  1:3241902/93‑1 (MQ=255)

                                           aaGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTa                                                 >  2:1304031/1‑93 (MQ=255)

                                           aaGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTa                                                 <  2:1391476/93‑1 (MQ=255)

                                            aGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTac                                                >  2:6301001/1‑93 (MQ=255)

                                             ggAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTacg                                               <  1:3957485/93‑1 (MQ=255)

                                              gAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTacga                                              <  2:6429364/93‑1 (MQ=255)

                                              gAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTacga                                              >  2:3011016/1‑93 (MQ=255)

                                                gTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTacgacg                                            >  2:2565943/1‑93 (MQ=255)

                                                gTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTacgacg                                            <  1:6645390/93‑1 (MQ=255)

                                                  tACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGc                                          >  2:4897161/1‑93 (MQ=255)

                                                  tACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTATGACGGc                                          <  2:4721141/93‑1 (MQ=255)

                                                   aCCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCa                                         <  1:2506795/93‑1 (MQ=255)

                                                     ccGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATc                                       >  2:5914734/1‑93 (MQ=255)

                                                       gATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGa                                     >  1:210605/1‑93 (MQ=255)

                                                       gATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGa                                     >  1:1229273/1‑93 (MQ=255)

                                                       gATGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGa                                     >  1:7443491/1‑93 (MQ=255)

                                                         tGAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGaaa                                   >  2:3082856/1‑93 (MQ=255)

                                                          gAAAACAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAAt                                  <  1:1030357/93‑1 (MQ=255)

                                                           aaaaCAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAAtt                                 <  2:6612675/93‑1 (MQ=255)

                                                            aaaCAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAAttt                                >  2:4852046/1‑93 (MQ=255)

                                                             aaCAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTg                               >  2:5497695/1‑93 (MQ=255)

                                                              aCAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTgg                              <  2:3642127/93‑1 (MQ=255)

                                                               cAGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTggg                             >  2:5162476/1‑93 (MQ=255)

                                                                aGGTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTgggg                            >  2:4607990/1‑93 (MQ=255)

                                                                 ggTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGc                           >  2:4684817/1‑93 (MQ=255)

                                                                 ggTATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGc                           >  2:321514/1‑93 (MQ=255)

                                                                   tATGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCgg                         <  1:5628586/93‑1 (MQ=255)

                                                                    aTGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGt                        >  1:5587925/1‑93 (MQ=255)

                                                                    aTGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGt                        >  1:5849232/1‑93 (MQ=255)

                                                                     tGTTTACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTc                       >  2:2208537/1‑93 (MQ=255)

                                                                        ttACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGcc                    >  2:4868114/1‑93 (MQ=255)

                                                                        ttACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGcc                    >  1:1201650/1‑93 (MQ=255)

                                                                         tACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGccc                   <  2:3952501/93‑1 (MQ=255)

                                                                         tACCTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGccc                   <  2:6051340/93‑1 (MQ=255)

                                                                          aCCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCg                  >  2:3135229/1‑93 (MQ=255)

                                                                           ccTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGc                 >  1:1876657/1‑93 (MQ=255)

                                                                            cTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCa                >  2:2360656/1‑93 (MQ=255)

                                                                            cTGCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCa                <  1:3584180/93‑1 (MQ=255)

                                                                             tGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCAc               >  2:5655498/1‑93 (MQ=255)

                                                                              gCGGCCACCAGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCAcg              <  1:858700/93‑1 (MQ=255)

                                                                                  ccaccaGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGtt          <  1:4175187/93‑1 (MQ=255)

                                                                                   caccaGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTc         >  2:4304742/1‑93 (MQ=255)

                                                                                   caccaGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTc         <  2:7432266/93‑1 (MQ=255)

                                                                                    accaGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTcg        >  1:7113433/1‑93 (MQ=255)

                                                                                    accaGCAGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTcg        >  2:3366602/1‑93 (MQ=255)

                                                                                      cagcagCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTcgtg      >  1:7188143/1‑91 (MQ=255)

                                                                                      cagcagCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTcgtg      >  1:2292821/1‑91 (MQ=255)

                                                                                       agcagcTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTCgtgc     <  2:3614224/93‑1 (MQ=255)

                                                                                       agcagcTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTCgtgc     <  2:6742979/93‑1 (MQ=255)

                                                                                          ggCTGCCCATTGAACACGCATTTCGTGCGGCAACCGAGCTGGGTTACGACGGCCTCGAAATTTGGGGCGGTCGCCCGCACGCGGTCGCGCCgg  >  1:4451652/1‑93 (MQ=255)

ATATTAATCCCTCTCCCGGCCCGGCAACGGGCCGGATTTACGCAAGGAGTTACCCGATGAAAACAGGTATGTTTACCTGCGGCCACCAGCGGCTGCCTATTGAACACGCATTTCGTGATGCAAGCGAGCTGGGTTACGACGGCATCGAAATTTGGGGCGGTCGCCCGCACGCGTTCGCGCCGG  >  NC_000913/3499921‑3500103

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