BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,675,076	+G	coding (920/2076 nt)	yhjG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,675,072	1	.	G	98.6%	1224.6 / 6.2	360	coding (924/2076 nt)	yhjG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/2); new base G (182/173); total (185/175)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 415 total aligned reads displayed.

TCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGCCC  >  NC_000913/3674980‑3675159

tCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc                                                                                            <  2:7321892/93‑1 (MQ=255)

tCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc                                                                                            <  1:5408834/93‑1 (MQ=255)

 cACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑g                                                                                          >  1:4391124/1‑93 (MQ=255)

 cACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑g                                                                                          >  1:3420577/1‑93 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                         <  2:8399919/93‑1 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                         <  2:1057054/93‑1 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                         >  1:4230672/1‑93 (MQ=255)

    cAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCCACCAGCCGACCATCCGTTTCAAACGGC‑gg‑gg                                                                                      >  1:1104355/1‑93 (MQ=255)

     aaTTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gg‑ggg                                                                                     <  2:8785284/93‑2 (MQ=255)

     aaTTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gg‑ggg                                                                                     >  2:3319552/1‑92 (MQ=255)

      aTTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gg‑gggt                                                                                    <  1:7590719/93‑3 (MQ=255)

      aTTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gg‑gggg                                                                                    >  1:2998687/1‑91 (MQ=255)

       ttCGCCCATTAAAACCGCGATAATCAAAGGCCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gg‑gggta                                                                                   <  2:2065440/93‑4 (MQ=255)

          gCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCg                                                                                <  1:5911372/93‑1 (MQ=255)

          gCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCg                                                                                <  2:7901513/93‑1 (MQ=255)

             cATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCgggg                                                                             >  2:576803/1‑92 (MQ=255)

             cATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGgc                                                                             >  1:5600041/1‑93 (MQ=255)

               ttAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGgcag                                                                           <  2:5096006/93‑1 (MQ=255)

               ttAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGgcag                                                                           >  1:5774205/1‑93 (MQ=255)

                 aaaaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGgcagca                                                                         >  2:7159984/1‑93 (MQ=255)

                 aaaaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGgcagca                                                                         >  2:2925384/1‑93 (MQ=255)

                   aaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGa                                                                       >  1:871570/1‑93 (MQ=255)

                    aCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGTAGaa                                                                      >  1:3540816/1‑93 (MQ=255)

                       gcgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGc                                                                   <  1:6333044/93‑1 (MQ=255)

                       gcgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGc                                                                   >  2:1379613/1‑93 (MQ=255)

                         gATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGccc                                                                 <  1:1948184/93‑1 (MQ=255)

                         gATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGccc                                                                 <  1:6565285/93‑1 (MQ=255)

                         gATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGccc                                                                 <  1:3361729/93‑1 (MQ=255)

                          aTAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCg                                                                <  2:6983215/93‑1 (MQ=255)

                           tAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGt                                                               <  2:1649208/93‑1 (MQ=255)

                           tAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGt                                                               >  2:5989026/1‑93 (MQ=255)

                           tAATCAAAGACCGACGATTTTTCAGTGTCGATTGTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGt                                                               <  1:8507368/93‑1 (MQ=255)

                            aaTCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTc                                                              >  1:4531243/1‑93 (MQ=255)

                            aaTCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTc                                                              >  1:14437/1‑93 (MQ=255)

                               cAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGt                                                           <  1:1266574/93‑1 (MQ=255)

                               cAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGt                                                           <  2:4151348/93‑1 (MQ=255)

                               aaaaGGCCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGt                                                           <  2:7595113/92‑1 (MQ=255)

                                aaaGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACCACCGTCAGtt                                                          >  1:2834727/1‑93 (MQ=255)

                                  aGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCa                                                        <  1:2075886/93‑1 (MQ=255)

                                  aGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCa                                                        <  1:8570870/93‑1 (MQ=255)

                                  aGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTCAGTTCa                                                        >  2:3746673/1‑93 (MQ=255)

                                   gACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCAt                                                       >  2:2608062/1‑93 (MQ=255)

                                   gACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCAt                                                       <  2:388640/93‑1 (MQ=255)

                                     ccGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATag                                                     >  2:8118793/1‑93 (MQ=255)

                                     ccGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTCAGTTCATaa                                                     >  2:6659030/1‑92 (MQ=255)

                                       gacgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATagag                                                   >  2:6144317/1‑93 (MQ=255)

                                       gacgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTCAGTTCATagag                                                   >  2:8420980/1‑93 (MQ=255)

                                         cgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGAt                                                 <  1:3928339/93‑1 (MQ=255)

                                          gaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATATAGATc                                                >  2:1227475/1‑93 (MQ=255)

                                          gaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATc                                                <  2:4051269/93‑1 (MQ=255)

                                           aTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCa                                               <  1:7090791/93‑1 (MQ=255)

                                            tttttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATATAGATCAc                                              <  2:2837278/93‑1 (MQ=255)

                                            tttttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCAc                                              <  1:8154532/93‑1 (MQ=255)

                                              tttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCAccc                                            <  2:5024723/93‑1 (MQ=255)

                                               ttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCa                                           <  1:2225163/93‑1 (MQ=255)

                                                tCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAg                                          <  1:7152780/93‑1 (MQ=255)

                                                tCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAg                                          <  1:549584/93‑1 (MQ=255)

                                                tCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAg                                          >  2:4775879/1‑93 (MQ=255)

                                                        gATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGc                                  >  2:5458878/1‑93 (MQ=255)

                                                         aTTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGcc                                 >  2:5383246/1‑93 (MQ=255)

                                                          ttttCGCTCCCCGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCa                                <  1:7043400/93‑1 (MQ=255)

                                                             tCGCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGaa                             >  2:4529149/1‑93 (MQ=255)

                                                               gCTCCCAGCCGACCATCCGTGTCAAACGGC‑GGAGG‑TATCGGGTAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGaaaa                           <  1:1193355/93‑1 (MQ=255)

                                                               gCTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGaaaa                           >  1:1626720/1‑93 (MQ=255)

                                                                cTACCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAt                          >  1:1392910/1‑93 (MQ=255)

                                                                  aCCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAttt                        <  2:2727495/93‑1 (MQ=255)

                                                                  aCCAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAttt                        >  2:19970/1‑93 (MQ=255)

                                                                   ccAGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTa                       <  2:7041661/93‑1 (MQ=255)

                                                                     aGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGGAGAACGCCCGTCAGTTCATAGAGACCACCCAGGGAATGGCCAGAAAATTTAAg                     >  2:7424326/1‑93 (MQ=255)

                                                                     aGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAg                     <  1:7490530/93‑1 (MQ=255)

                                                                     aGCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAg                     <  2:3561746/93‑1 (MQ=255)

                                                                      gCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGc                    <  1:1005388/93‑1 (MQ=255)

                                                                      gCCGACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGc                    <  1:6692601/93‑1 (MQ=255)

                                                                         gACCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGt                 >  1:7273028/1‑93 (MQ=255)

                                                                          aCCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTa                >  2:3446958/1‑93 (MQ=255)

                                                                          aCCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTa                >  1:2510012/1‑93 (MQ=255)

                                                                          aCCATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTa                >  1:7387384/1‑93 (MQ=255)

                                                                           ccATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaa               <  1:8298504/93‑1 (MQ=255)

                                                                           ccATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaa               <  1:7946070/93‑1 (MQ=255)

                                                                            cATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa              >  2:1955091/1‑93 (MQ=255)

                                                                            cATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa              >  2:5600164/1‑93 (MQ=255)

                                                                            cATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa              <  1:7766277/93‑1 (MQ=255)

                                                                            cATCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa              <  1:5977472/93‑1 (MQ=255)

                                                                             aTCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTCAAt             >  2:4830459/1‑93 (MQ=255)

                                                                             aTCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAAt             <  2:8248668/93‑1 (MQ=255)

                                                                             aTCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAAt             <  1:4656065/93‑1 (MQ=255)

                                                                             aTCCGTTTCAAACGGCGGG‑GG‑AATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAAt             >  1:324516/1‑93 (MQ=255)

                                                                              tCCGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc            >  1:552136/1‑93 (MQ=255)

                                                                              tCCGTTTCAAACGGCGGG‑GG‑GATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc            >  2:422880/1‑93 (MQ=255)

                                                                              tCCGTTTCAAACGGC‑GG‑GGTTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc            <  1:4882884/93‑1 (MQ=255)

                                                                               ccGTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCg           >  1:2560621/1‑93 (MQ=255)

                                                                                 gTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGac         <  2:3158093/93‑1 (MQ=255)

                                                                                 gTTTCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGCCAGCTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGac         >  2:4653902/1‑93 (MQ=255)

                                                                                  tttCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGaca        <  1:6041667/93‑1 (MQ=255)

                                                                                  tttCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGaca        <  1:5193332/93‑1 (MQ=255)

                                                                                  tttCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGaca        >  1:4029654/1‑93 (MQ=255)

                                                                                  gttCAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGaca        <  2:6232393/92‑1 (MQ=255)

                                                                                     cAAACGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCg     <  2:2340945/93‑1 (MQ=255)

                                                                                       aaCGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGcc   <  2:7761659/93‑1 (MQ=255)

                                                                                       aaCGGCGGG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATAGACACCGcc   >  1:6848066/1‑93 (MQ=255)

TCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GG‑GG‑TATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGCCC  >  NC_000913/3674980‑3675159

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