BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,665,031	A→G	V279A (GTG→GCG)	hcaT ←	HcaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,665,031	0	A	G	99.6%	702.1 / ‑4.9	224	V279A (GTG→GCG)	hcaT	HcaT
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (91/132); total (92/132)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 240 total aligned reads displayed.

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCC  >  NC_000913/2664941‑2665120

 ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTa                                                                                        >  2:6667776/1‑93 (MQ=255)

 ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTa                                                                                        >  1:5491232/1‑93 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                       >  1:8527491/1‑93 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                       <  2:2121586/93‑1 (MQ=255)

      cAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAcgccg                                                                                   >  2:8671362/1‑93 (MQ=255)

       aGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAcgccgc                                                                                  <  2:6220859/93‑1 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                 >  1:8622749/1‑93 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                 >  2:334289/1‑93 (MQ=255)

         aCCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCaa                                                                                >  2:6255824/1‑93 (MQ=255)

          ccGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCaaa                                                                               <  1:9421737/93‑1 (MQ=255)

          ccGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCaaa                                                                               <  2:3590090/93‑1 (MQ=255)

           cGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAAt                                                                              <  2:2875081/93‑1 (MQ=255)

            gTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATc                                                                             >  2:8410534/1‑93 (MQ=255)

             tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCg                                                                            >  2:2613080/1‑93 (MQ=255)

                aGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCg                                                                         <  1:6408251/93‑1 (MQ=255)

                 ggTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCTAATCGCCga                                                                        <  2:6738158/93‑1 (MQ=255)

                 ggTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCga                                                                        <  2:1263378/93‑1 (MQ=255)

                  gTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCgag                                                                       >  2:7589869/1‑93 (MQ=255)

                   tGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCgaga                                                                      <  1:4896040/93‑1 (MQ=255)

                       gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCaa                                                                  <  1:1008735/93‑1 (MQ=255)

                       gCAAGGCAGAATTTCCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCaa                                                                  <  1:8217271/93‑1 (MQ=255)

                        cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAAc                                                                 <  2:10043360/93‑1 (MQ=255)

                         aaTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGCTCAACa                                                                >  1:8109176/1‑93 (MQ=255)

                           tGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGc                                                              >  1:6858669/1‑93 (MQ=255)

                            gCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCa                                                             >  1:9909877/1‑93 (MQ=255)

                            gCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCa                                                             <  2:7360685/93‑1 (MQ=255)

                              aGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCata                                                           <  1:3026446/93‑1 (MQ=255)

                              aGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCata                                                           <  2:9219267/93‑1 (MQ=255)

                              aGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATG‑CCCAGCGCGCTACGCCGCCAATCGCCGAGATCAACAGCatat                                                          >  1:3074590/1‑93 (MQ=255)

                               gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                          <  1:923634/93‑1 (MQ=255)

                               gAATTTGCACCACTATCAACCACGACAACGCCATAGTTACACCCATAATGCCCCAGCGCGCAACGCAGCAAATCACCGAGATCAACCGCaaat                                                          <  2:7718539/93‑1 (MQ=255)

                               gAATTGGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                          <  1:4857127/93‑1 (MQ=255)

                                aaTTTTCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATc                                                         <  2:4990413/93‑1 (MQ=255)

                                aaTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATc                                                         <  2:187689/93‑1 (MQ=255)

                                 aTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcg                                                        <  1:2628879/93‑1 (MQ=255)

                                  tttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                       >  2:9834506/1‑93 (MQ=255)

                                  tttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                       >  1:4814772/1‑93 (MQ=255)

                                       cccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                  <  2:2230483/92‑1 (MQ=255)

                                       cccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                  <  1:4195026/92‑1 (MQ=255)

                                       accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                  >  1:7184408/1‑93 (MQ=255)

                                        ccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCAc                                                 >  2:1799772/1‑93 (MQ=255)

                                        ccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCAc                                                 <  1:2878247/93‑1 (MQ=255)

                                           cTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                              >  2:9673967/1‑93 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                             >  1:5038681/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                            >  1:27359/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                            <  2:6458034/93‑1 (MQ=255)

                                               cAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcg                                          <  2:8145192/93‑1 (MQ=255)

                                                 aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                        >  1:2085766/1‑93 (MQ=255)

                                                  ccACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcg                                       >  2:8012453/1‑93 (MQ=255)

                                                   cACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcgg                                      >  2:6787008/1‑93 (MQ=255)

                                                    aCGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGa                                     <  1:4316833/93‑1 (MQ=255)

                                                      ggCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaa                                   >  1:9648813/1‑93 (MQ=255)

                                                       gCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaaa                                  >  1:1197179/1‑93 (MQ=255)

                                                         aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGt                                <  1:446837/93‑1 (MQ=255)

                                                         aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGt                                >  2:9597/1‑93 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                               <  2:5720565/93‑1 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                               >  2:7270641/1‑93 (MQ=255)

                                                           cGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGttt                              <  1:7858304/93‑1 (MQ=255)

                                                           cGCCGTAGTGCCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGttt                              <  2:3041217/93‑1 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                            <  1:1734022/93‑1 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttattc                          >  2:5274523/1‑92 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                          >  1:6550829/1‑93 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                          <  1:2662558/93‑1 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                          >  1:9148140/1‑93 (MQ=255)

                                                                tAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTAc                         <  2:4325327/93‑1 (MQ=255)

                                                                tAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTAc                         >  2:911546/1‑93 (MQ=255)

                                                                 aGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                        <  2:4429418/93‑1 (MQ=255)

                                                                 aGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                        >  1:3378493/1‑93 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                       >  1:4576507/1‑93 (MQ=255)

                                                                   ttGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTATTcc                      >  1:3142838/1‑92 (MQ=255)

                                                                    tGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAg                     <  1:5284442/93‑1 (MQ=255)

                                                                     gCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgc                    <  1:6075327/93‑1 (MQ=255)

                                                                       taccATAATGCCTCAGCGCGCTGCGCCGCAGATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgcgc                  <  1:9962101/91‑1 (MQ=255)

                                                                        cccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCa                 >  1:8728325/1‑93 (MQ=255)

                                                                         ccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaa                >  2:7887815/1‑93 (MQ=255)

                                                                         ccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaa                <  2:5010531/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  2:1660977/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  2:6622807/93‑1 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGa             >  2:5879902/1‑93 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAAGGCGGCAAAGTTTATTACTCAGCGCAAAGa             >  2:5413545/1‑93 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa           <  1:4417456/93‑1 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa           <  2:852234/93‑1 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa           <  2:3814370/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  1:2571841/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  1:2418394/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  1:5591586/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  2:8805112/93‑1 (MQ=255)

                                                                                 ccccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATg        <  1:2917113/93‑1 (MQ=255)

                                                                                   ccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGAc      >  2:7603938/1‑93 (MQ=255)

                                                                                   ccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGAc      <  2:1773908/93‑1 (MQ=255)

                                                                                   ccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGAc      <  2:5455239/93‑1 (MQ=255)

                                                                                    cAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACt     <  1:7614661/93‑1 (MQ=255)

                                                                                    cAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACt     >  1:1480072/1‑93 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt    <  2:4541482/93‑1 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt    >  2:8237751/1‑93 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt    <  1:1296552/93‑1 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt    <  2:4862039/93‑1 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc   >  1:3218001/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc   <  2:2322052/93‑1 (MQ=255)

                                                                                       cgcgCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTcc  <  1:742484/93‑1 (MQ=255)

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCC  >  NC_000913/2664941‑2665120

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