BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,125,339	C→T	D46N (GAT→AAT)	yiiX ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,125,339	0	C	T	99.2%	1323.5 / ‑3.6	368	D46N (GAT→AAT)	yiiX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/2); new base T (187/177); total (188/180)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 374 total aligned reads displayed.

GCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGG  >  NC_000913/4125249‑4125428

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       >  1:5938929/1‑93 (MQ=255)

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       <  1:2236493/93‑1 (MQ=255)

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       <  2:1934897/93‑1 (MQ=255)

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       <  2:2511001/93‑1 (MQ=255)

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       <  2:4087744/93‑1 (MQ=255)

  ttGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTa                                                                                       >  1:335682/1‑93 (MQ=255)

   tGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTAt                                                                                      <  1:4809422/93‑1 (MQ=255)

   tGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTAt                                                                                      <  2:816334/93‑1 (MQ=255)

   tGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTAt                                                                                      <  1:242741/93‑1 (MQ=255)

   tGAGCGGGGGGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTAt                                                                                      <  2:4376830/93‑1 (MQ=255)

    gAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATg                                                                                     <  2:6940923/93‑1 (MQ=255)

     aGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATgg                                                                                    <  1:7358921/93‑1 (MQ=255)

     aGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATgg                                                                                    <  2:6383952/93‑1 (MQ=255)

       cGGGGTGTATTCCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGt                                                                                  <  2:585343/93‑1 (MQ=255)

        ggggTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTc                                                                                 >  1:1329847/1‑93 (MQ=255)

          ggTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGc                                                                               >  1:1023761/1‑93 (MQ=255)

          ggTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGc                                                                               <  2:5712907/93‑1 (MQ=255)

           gtgtATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGcc                                                                              >  2:1701439/1‑93 (MQ=255)

           gtgtATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGcc                                                                              >  2:1472231/1‑93 (MQ=255)

            tgtATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCa                                                                             >  1:400499/1‑93 (MQ=255)

              tATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGt                                                                           >  1:4410105/1‑93 (MQ=255)

              tATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGt                                                                           >  1:78528/1‑93 (MQ=255)

              tATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGt                                                                           >  1:4661074/1‑93 (MQ=255)

              tATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGt                                                                           >  1:5745420/1‑93 (MQ=255)

               aTTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGtt                                                                          <  2:942498/93‑1 (MQ=255)

                tttCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTg                                                                         >  2:4645436/1‑93 (MQ=255)

                   cACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCTAGTTGGAt                                                                      <  2:7635848/93‑1 (MQ=255)

                   cACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGAt                                                                      >  2:3035210/1‑93 (MQ=255)

                   cACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGAt                                                                      <  2:968576/93‑1 (MQ=255)

                      cGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGc                                                                   <  2:5349697/93‑1 (MQ=255)

                        ggCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCtt                                                                 >  1:5307916/1‑93 (MQ=255)

                        ggCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCtt                                                                 >  2:724356/1‑93 (MQ=255)

                         gCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCttt                                                                >  1:7908465/1‑93 (MQ=255)

                          ccGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTa                                                               <  2:2827533/93‑1 (MQ=255)

                            gACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACt                                                             >  2:5208480/1‑93 (MQ=255)

                             aCTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACtt                                                            <  2:1998543/93‑1 (MQ=255)

                               tGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTg                                                          <  1:7605794/93‑1 (MQ=255)

                                gCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCGCCGGCATTCCGGTGTGGCTATAATTGGTATGGATCGCCAGTTGGATCGCTTTACTTTGc                                                         <  2:7023282/93‑1 (MQ=255)

                                gCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGTGTCGCCAGTTGGATCGCTTTACTTTGc                                                         >  1:6654986/1‑93 (MQ=255)

                                gCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGc                                                         >  2:4599070/1‑93 (MQ=255)

                                gCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGc                                                         >  2:7301002/1‑93 (MQ=255)

                                  ttCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGa                                                       >  2:4452525/1‑93 (MQ=255)

                                  ttCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGa                                                       <  1:6395125/93‑1 (MQ=255)

                                  ttCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGa                                                       >  2:6700552/1‑93 (MQ=255)

                                    cAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGCGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATg                                                     >  1:559450/1‑93 (MQ=255)

                                     aaaaaCGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGa                                                    <  2:3459474/93‑1 (MQ=255)

                                      aaaaCGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAg                                                   <  2:1623119/93‑1 (MQ=255)

                                      aaaaCGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAg                                                   <  2:7445891/93‑1 (MQ=255)

                                       aaaCGTAGTGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGc                                                  <  2:4791798/93‑1 (MQ=255)

                                       aaaCGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGc                                                  <  2:2254883/93‑1 (MQ=255)

                                           gTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCgaga                                              <  1:4587510/93‑1 (MQ=255)

                                           gTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCgaga                                              >  1:472865/1‑93 (MQ=255)

                                           gTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCgaga                                              >  1:6306966/1‑93 (MQ=255)

                                            tAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCgagag                                             >  1:5217629/1‑93 (MQ=255)

                                             aGGGTTTTTTGTTGCGTATCCCCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCgagaga                                            <  2:493120/93‑1 (MQ=255)

                                              gggCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGAt                                           >  1:251298/1‑93 (MQ=255)

                                              gggCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGAt                                           <  2:2666356/93‑1 (MQ=255)

                                               ggCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATc                                          <  1:2766961/93‑1 (MQ=255)

                                               ggCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATc                                          <  1:3722465/93‑1 (MQ=255)

                                                gCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCt                                         >  1:2323026/1‑93 (MQ=255)

                                                 cTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTg                                        <  2:7902845/93‑1 (MQ=255)

                                                  tttttGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGa                                       >  2:5584117/1‑93 (MQ=255)

                                                  tttttGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGa                                       >  1:1668755/1‑93 (MQ=255)

                                                   ttttGTTGCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGaa                                      >  2:3845241/1‑93 (MQ=255)

                                                     ttgttgCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAg                                    >  1:5417483/1‑93 (MQ=255)

                                                       gttgCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGCTCGCTTTACTCTGCGATGAGCGAGAGATCTGAAAgat                                  >  2:6975434/1‑93 (MQ=255)

                                                       gttgCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAgat                                  <  1:3507747/93‑1 (MQ=255)

                                                        ttgCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAgatg                                 >  1:2941973/1‑93 (MQ=255)

                                                        ttgCGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAgatg                                 <  1:983125/93‑1 (MQ=255)

                                                           cGTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGatgatg                              <  2:6249318/93‑1 (MQ=255)

                                                            gTATCACCAGCATACCGGTGTTTCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGt                             >  1:1412939/1‑93 (MQ=255)

                                                            gTATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGt                             >  1:6969644/1‑93 (MQ=255)

                                                             tATCACCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTc                            <  2:2362434/93‑1 (MQ=255)

                                                                 aCCAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCg                        <  1:1925065/93‑1 (MQ=255)

                                                                  ccAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCgg                       <  2:1217432/93‑1 (MQ=255)

                                                                  ccAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCgg                       >  1:6783533/1‑93 (MQ=255)

                                                                   cAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGt                      >  2:6557399/1‑93 (MQ=255)

                                                                   cAGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGt                      >  1:3043058/1‑93 (MQ=255)

                                                                    aGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGtt                     >  1:3857488/1‑93 (MQ=255)

                                                                    aGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGtt                     >  1:6046327/1‑93 (MQ=255)

                                                                    aGCATACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGtt                     <  2:6582007/93‑1 (MQ=255)

                                                                       aTACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTtgt                  >  1:7160508/1‑93 (MQ=255)

                                                                        tACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTtgtg                 >  1:4922608/1‑93 (MQ=255)

                                                                        tACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTtgtg                 >  2:1904871/1‑93 (MQ=255)

                                                                        tACCGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTtgtg                 >  2:3648155/1‑93 (MQ=255)

                                                                          ccGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGc               <  1:6212274/93‑1 (MQ=255)

                                                                          ccGGTGTGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGc               <  1:1236268/93‑1 (MQ=255)

                                                                             gtgtGGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGc            >  2:5478197/1‑93 (MQ=255)

                                                                              tgtgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGcc           >  1:2530649/1‑93 (MQ=255)

                                                                              tgtgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGcc           >  1:660759/1‑93 (MQ=255)

                                                                                tgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAg         >  1:1141756/1‑93 (MQ=255)

                                                                                tgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAg         >  2:6108584/1‑93 (MQ=255)

                                                                                tgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAg         >  2:3153945/1‑93 (MQ=255)

                                                                                tgGCTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAg         >  2:3331519/1‑93 (MQ=255)

                                                                                   cTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCa      >  1:7416887/1‑93 (MQ=255)

                                                                                   cTATAATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCa      >  2:5010870/1‑93 (MQ=255)

                                                                                    tataATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCaa     >  1:1807359/1‑93 (MQ=255)

                                                                                    tataATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCaa     >  2:6906697/1‑93 (MQ=255)

                                                                                     ataATTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCaaa    >  1:4186227/1‑93 (MQ=255)

                                                                                       aaTTGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAgg  <  2:667754/93‑1 (MQ=255)

GCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGG  >  NC_000913/4125249‑4125428

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