BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,953,910	A→G	N320S (AAC→AGC)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,953,910	0	A	G	99.4%	1120.1 / ‑4.9	355	N320S (AAC→AGC)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (167/185); total (167/188)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 383 total aligned reads displayed.

CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999

cgctggcgctggCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCt                                                                                         <  2:8668621/93‑1 (MQ=255)

 gctggcgctggcGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCtt                                                                                        >  1:3990154/1‑93 (MQ=255)

  ctggcgctgtggGGAATGAAAAAATATATCGCACTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTc                                                                                       <  2:5547188/93‑1 (MQ=255)

  ctggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTc                                                                                       <  1:4841205/93‑1 (MQ=255)

   tggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTcc                                                                                      >  1:5805059/1‑93 (MQ=255)

    ggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCa                                                                                     <  2:2383292/93‑1 (MQ=255)

     gcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCAc                                                                                    <  1:3288482/93‑1 (MQ=255)

     gcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCAc                                                                                    >  2:2605466/1‑93 (MQ=255)

     gcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATAATCTTTCCGGTGCCAACGTGAGCTTCCAc                                                                                    >  2:6454095/1‑93 (MQ=255)

       gctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACgg                                                                                  <  2:9831922/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  1:8061531/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                >  2:2030372/1‑93 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  1:4952933/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  2:1682937/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                >  2:7919721/1‑93 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  1:9838108/93‑1 (MQ=255)

           gcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTg                                                                              <  1:6477380/93‑1 (MQ=255)

            cgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgc                                                                             <  2:6941019/93‑1 (MQ=255)

              ggAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgcgc                                                                           >  2:604518/1‑93 (MQ=255)

              ggAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgcgc                                                                           <  1:3618188/93‑1 (MQ=255)

               gAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCt                                                                          >  2:7413961/1‑93 (MQ=255)

                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         <  1:9926187/93‑1 (MQ=255)

                 aTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGCACGGCTGGCGCATATTTTTTTCTGCGCCAACCTGGGCTTCCACGTCCTTCGCTAc                                                                        >  2:7179748/1‑93 (MQ=255)

                 aTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTAc                                                                        <  1:8569157/93‑1 (MQ=255)

                  tGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACg                                                                       >  1:1065565/1‑93 (MQ=255)

                  tGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACg                                                                       >  2:675207/1‑93 (MQ=255)

                  tGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCGACg                                                                       >  2:3510313/1‑93 (MQ=255)

                     aaaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtct                                                                    >  1:4418746/1‑93 (MQ=255)

                      aaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtctc                                                                   >  1:1420372/1‑93 (MQ=255)

                        aaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATCTTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTCCGtctctg                                                                 >  2:9616724/1‑91 (MQ=255)

                        aaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAg                                                                 >  2:2701632/1‑93 (MQ=255)

                         atatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGa                                                                <  1:5570263/93‑1 (MQ=255)

                         atatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCACATTCTTTCCGGTGCCAACGTGAGCTTCCACAGCCTGCGCTACGTCTCCGa                                                                >  1:4095884/1‑93 (MQ=255)

                          tatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGaa                                                               <  1:5307097/93‑1 (MQ=255)

                          tatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGaa                                                               >  2:1760436/1‑93 (MQ=255)

                          tatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGaa                                                               <  2:5131439/93‑1 (MQ=255)

                           atatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAAc                                                              >  1:8997976/1‑93 (MQ=255)

                            tatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACg                                                             >  2:2893135/1‑93 (MQ=255)

                             atCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGc                                                            <  2:6297547/93‑1 (MQ=255)

                              tCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCt                                                           <  2:5226824/93‑1 (MQ=255)

                                 cccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCg                                                        >  2:3183966/1‑93 (MQ=255)

                                 cccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGACCTTCCACGGCCTCCGCTACGTCTCAGAACGCTGCg                                                        <  1:1833500/93‑1 (MQ=255)

                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  2:5935021/1‑93 (MQ=255)

                                   cTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGaa                                                      <  2:9441979/93‑1 (MQ=255)

                                   cTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGaa                                                      <  1:4260763/93‑1 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGCAc                                                     >  2:3707788/1‑93 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     <  1:10041960/93‑1 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     >  2:3307645/1‑93 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     >  1:9569092/1‑93 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     >  1:8330989/1‑93 (MQ=255)

                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     <  1:8597058/93‑1 (MQ=255)

                                     gCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACt                                                    <  1:6908418/93‑1 (MQ=255)

                                     gCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACt                                                    <  2:4494562/93‑1 (MQ=255)

                                      cacaACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTg                                                   <  1:877949/93‑1 (MQ=255)

                                       acaacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTgg                                                  <  2:7984556/93‑1 (MQ=255)

                                       acaacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTgg                                                  >  2:4568958/1‑93 (MQ=255)

                                             ttCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAAc                                            <  2:9065540/93‑1 (MQ=255)

                                              tCGCGGCGAACGGCTGGCGCATATTCTTCCCGGTGCCAACGTGAGTTTCCACGGCCTGCGCTACGTCCCAGAACGCTGCGAACTGGGCGAACa                                           <  2:1123094/93‑1 (MQ=255)

                                               cgcgGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAg                                          <  2:4177739/93‑1 (MQ=255)

                                                  ggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGACCAGCgg                                       >  1:593556/1‑92 (MQ=255)

                                                  ggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGt                                       >  2:1233261/1‑93 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGTTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    <  1:2071962/93‑1 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    >  2:43209/1‑93 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    <  2:3869006/93‑1 (MQ=255)

                                                      accGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   <  1:5766190/91‑1 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   >  2:8936384/1‑93 (MQ=255)

                                                        cGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCg                                 <  1:280593/93‑1 (MQ=255)

                                                        cGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGATTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCg                                 >  1:169512/1‑93 (MQ=255)

                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  1:7915873/93‑1 (MQ=255)

                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  1:1518357/93‑1 (MQ=255)

                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgtt                            >  1:9793701/1‑93 (MQ=255)

                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgtt                            >  1:126073/1‑93 (MQ=255)

                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGGGAAGCgttgtt                            <  2:841117/93‑1 (MQ=255)

                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           >  1:9097806/1‑93 (MQ=255)

                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           <  2:9270770/93‑1 (MQ=255)

                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           >  2:2197230/1‑93 (MQ=255)

                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           <  1:7262915/93‑1 (MQ=255)

                                                                gcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGc                         <  2:8153063/93‑1 (MQ=255)

                                                                gcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACCTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTggg                         >  2:4882581/1‑92 (MQ=255)

                                                                 cATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCg                        >  1:2949775/1‑93 (MQ=255)

                                                                 cATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCg                        >  1:343091/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      >  2:879982/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      >  1:9364702/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      >  1:6690872/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGAGAACTGGGCGAAAAGCGATAAGCGTTGTTGGCggg                      >  2:5705500/1‑92 (MQ=255)

                                                                      tCTTTCCGGTGCCAACGTGAGCTGCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGAc                   <  2:9149019/93‑1 (MQ=255)

                                                                        tttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCa                 <  2:1473668/93‑1 (MQ=255)

                                                                         ttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAt                <  1:7637985/93‑1 (MQ=255)

                                                                          tCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAtt               <  1:2590347/93‑1 (MQ=255)

                                                                           ccGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTc              <  2:7763557/93‑1 (MQ=255)

                                                                            cGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTcc             <  1:9264077/93‑1 (MQ=255)

                                                                                 ccAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaag        <  1:6985727/93‑1 (MQ=255)

                                                                                 ccAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGGGAAGCGTTGTTGGCGGTGACCATTCCGgaag        >  1:2310230/1‑93 (MQ=255)

                                                                                  cAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaaga       <  2:5385322/93‑1 (MQ=255)

                                                                                  cAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaaga       >  1:6596072/1‑93 (MQ=255)

                                                                                    aCGTGATCTCACATGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     <  1:3700757/93‑1 (MQ=255)

                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   <  1:7422224/93‑1 (MQ=255)

                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   <  1:135019/93‑1 (MQ=255)

                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   <  2:5831364/93‑1 (MQ=255)

                                                                                       tGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAg  <  2:4732559/93‑1 (MQ=255)

CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999

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