BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	380,496	(C)_5→6	coding (16/444 nt)	b0359 ←	putative transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,491	1	.	C	90.1%	1009.4 / 104.3	356	coding (21/444 nt)	b0359	putative transferase
Reads supporting (aligned to +/- strand): ref base . (12/23); new base C (152/168); total (165/191)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-02

Only 100 of 422 total aligned reads displayed.

ATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAC  >  NC_000913/380399‑380572

aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                      >  1:1881259/1‑93 (MQ=255)

 tAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑c                                                                                    <  1:1173620/93‑1 (MQ=255)

  aaCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc                                                                                   <  1:5153361/93‑1 (MQ=255)

   aCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑c                                                                                 <  1:5511243/93‑1 (MQ=255)

   aCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑c                                                                                 <  2:7546572/93‑1 (MQ=255)

   aCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑c                                                                                 <  2:6041848/93‑1 (MQ=255)

      gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccc                                                                              >  1:9947798/1‑92 (MQ=255)

      gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccc                                                                              >  2:3339133/1‑92 (MQ=255)

      gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccc                                                                              >  2:4058190/1‑92 (MQ=255)

      gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccc                                                                              >  2:6525678/1‑92 (MQ=255)

       cATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccct                                                                             <  2:9069212/93‑3 (MQ=255)

         tttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccctga                                                                           <  2:678368/93‑5 (MQ=255)

         tttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑cc‑cccctga                                                                           >  1:3537358/1‑89 (MQ=255)

            ttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAAtc                                                                        >  2:5698923/1‑93 (MQ=255)

            ttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAAtc                                                                        >  1:8386139/1‑93 (MQ=255)

             tAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAAtct                                                                       >  2:8353689/1‑93 (MQ=255)

             tAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAAtct                                                                       >  2:1916218/1‑93 (MQ=255)

               aTCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTcc                                                                     >  2:3488356/1‑93 (MQ=255)

               aTCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTcc                                                                     >  1:9191651/1‑93 (MQ=255)

                  tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGa                                                                  >  1:1779935/1‑93 (MQ=255)

                  tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGa                                                                  >  1:9789383/1‑93 (MQ=255)

                      tGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAAc                                                              >  1:3918322/1‑93 (MQ=255)

                       ggCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAAca                                                             >  2:9297562/1‑92 (MQ=255)

                          aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAat                                                          >  1:9275325/1‑93 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAata                                                         <  1:7427305/93‑1 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAata                                                         <  2:942326/93‑1 (MQ=255)

                           aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAata                                                         >  2:8537586/1‑93 (MQ=255)

                            aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAatat                                                        >  2:9854400/1‑93 (MQ=255)

                               cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCAc                                                     <  1:2148614/93‑1 (MQ=255)

                                attttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACtt                                                   <  2:4674094/91‑1 (MQ=255)

                                aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACt                                                    <  2:3690863/93‑1 (MQ=255)

                                aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACt                                                    <  2:7487072/93‑1 (MQ=255)

                                aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACtt                                                   <  2:3598707/93‑1 (MQ=255)

                                 aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACtt                                                   <  2:176161/93‑1 (MQ=255)

                                  tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTa                                                  >  1:2832221/1‑93 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTaaa                                                <  2:5273864/93‑1 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTaaa                                                <  1:8854216/93‑1 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTaaa                                                <  2:6174888/93‑1 (MQ=255)

                                    tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTaaa                                                <  1:5644412/93‑1 (MQ=255)

                                     ttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAAt                                               <  2:8712251/93‑1 (MQ=255)

                                     ttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAAt                                               <  2:6984663/93‑1 (MQ=255)

                                     ttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATa                                              <  1:1022596/93‑1 (MQ=255)

                                        cGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAg                                            >  1:4365740/1‑93 (MQ=255)

                                         gCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGt                                           <  1:7210823/93‑1 (MQ=255)

                                          ccATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATACGTg                                          >  1:83150/1‑93 (MQ=255)

                                           cATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                         <  2:5879582/93‑1 (MQ=255)

                                            aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCCTTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                        <  2:5029814/93‑1 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       <  2:6377723/93‑1 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       <  2:3905429/93‑1 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       >  1:6791367/1‑93 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       >  2:9539631/1‑93 (MQ=255)

                                             ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                       >  1:9429271/1‑93 (MQ=255)

                                                cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                    <  2:7257953/93‑1 (MQ=255)

                                                  gATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                  >  1:8623075/1‑93 (MQ=255)

                                                  gATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                  <  1:8117611/93‑1 (MQ=255)

                                                   attattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                 >  1:3409818/1‑93 (MQ=255)

                                                    ttattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                <  2:9633755/93‑1 (MQ=255)

                                                     tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                               <  1:5410982/93‑1 (MQ=255)

                                                      attGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                              >  1:650766/1‑93 (MQ=255)

                                                      attGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                              >  1:664952/1‑93 (MQ=255)

                                                      attGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                             >  2:4466744/1‑93 (MQ=255)

                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                           >  1:5922732/1‑93 (MQ=255)

                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                           <  2:5574063/93‑1 (MQ=255)

                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                           <  1:4841694/93‑1 (MQ=255)

                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                           <  2:8772479/93‑1 (MQ=255)

                                                          tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                          <  2:134220/93‑1 (MQ=255)

                                                           gCCAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                         <  1:2324538/93‑1 (MQ=255)

                                                            ccAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                        >  2:8025253/1‑93 (MQ=255)

                                                            ccAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                        <  1:1138691/93‑1 (MQ=255)

                                                            ccAGGACGAATAAACGCATAATTACTAGACTGG‑CCACCC‑TGAATCTCCAGACAACCAATATCACTTAAATAATTGATAGTCTTAATACTAGtt                        <  2:8897886/93‑1 (MQ=255)

                                                             cAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                       >  1:466441/1‑93 (MQ=255)

                                                             cAGGACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                       <  2:7808826/93‑1 (MQ=255)

                                                               ggACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                     >  1:2562842/1‑93 (MQ=255)

                                                                gACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                    <  2:3093403/93‑1 (MQ=255)

                                                                gACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                    <  2:1581080/93‑1 (MQ=255)

                                                                gACGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                    <  1:4177659/93‑1 (MQ=255)

                                                                 aCGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                   >  2:2195204/1‑93 (MQ=255)

                                                                 aCGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                   >  2:1992323/1‑93 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  <  1:1275710/93‑1 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  <  2:1125327/93‑1 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  >  2:9099791/1‑93 (MQ=255)

                                                                  cGAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                  >  1:9665430/1‑93 (MQ=255)

                                                                   gAATAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                 >  2:6713378/1‑93 (MQ=255)

                                                                    aaTAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                <  2:2406161/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               >  1:640015/1‑93 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               <  1:5758390/93‑1 (MQ=255)

                                                                     aTAAACGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa               <  2:6466957/93‑1 (MQ=255)

                                                                        aaCGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc            >  1:8002397/1‑93 (MQ=255)

                                                                         aCGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCa           >  1:3974519/1‑93 (MQ=255)

                                                                          cGCATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt          <  2:10043652/93‑1 (MQ=255)

                                                                            cATAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg        >  1:1452403/1‑93 (MQ=255)

                                                                             aTAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg       >  1:6257242/1‑93 (MQ=255)

                                                                             aTAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg       <  1:6141021/93‑1 (MQ=255)

                                                                             aTAATTACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg       <  1:2961320/93‑1 (MQ=255)

                                                                                 ttACTAGACTGGCCC‑CCC‑TGAATCTCCAGACCACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa   >  1:9050814/1‑93 (MQ=255)

                                                                                 ttACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa   <  2:8264999/93‑1 (MQ=12)

                                                                                 ttACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa   >  2:2173475/1‑93 (MQ=12)

                                                                                  tACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc  >  1:6998980/1‑93 (MQ=9)

                                                                                  tACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc  <  1:5076007/93‑1 (MQ=9)

                                                                                  tACTAGACTGGCCC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc  >  2:7106464/1‑93 (MQ=9)

ATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGG‑CC‑CCC‑TGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAC  >  NC_000913/380399‑380572

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