BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,453,670	T→C	T152T (ACT→ACC)	yheE →	YheE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,453,670	0	T	C	99.0%	740.2 / NA	208	T152T (ACT→ACC)	yheE	YheE
Reads supporting (aligned to +/- strand): ref base T (0/1); major base C (117/87); minor base A (0/2); total (118/90)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

AGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATAC  >  NC_000913/3453580‑3453758

aGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAAccct                                                                                        >  1:1097464/1‑93 (MQ=255)

    cTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCtccat                                                                                    <  1:5955892/93‑1 (MQ=255)

     ttGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCtccatc                                                                                   >  1:3350500/1‑93 (MQ=255)

     ttGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCtccatc                                                                                   >  2:3818795/1‑93 (MQ=255)

         gTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAgg                                                                               <  1:5705369/93‑1 (MQ=255)

          tCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAggg                                                                              >  2:1772494/1‑93 (MQ=255)

           cAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCTCGGGCATTACGAAAAAGTAACCCTCCATCCAGGGc                                                                             >  1:852546/1‑93 (MQ=255)

            aGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCt                                                                            <  2:5447753/93‑1 (MQ=255)

            aGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCt                                                                            <  1:3050554/93‑1 (MQ=255)

             gcgcCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTa                                                                           >  2:5627219/1‑93 (MQ=255)

               gcagggAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTAcc                                                                         <  2:4624926/90‑1 (MQ=255)

               gcCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTAcc                                                                         >  2:5986672/1‑93 (MQ=255)

               gcCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTAcc                                                                         >  1:6141486/1‑93 (MQ=255)

               gcCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTAcc                                                                         >  2:157707/1‑93 (MQ=255)

                 cGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTg                                                                       >  2:2529475/1‑93 (MQ=255)

                  gggAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGa                                                                      >  2:1532022/1‑93 (MQ=255)

                  gggAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGCTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGa                                                                      >  2:26209/1‑93 (MQ=255)

                   ggAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAc                                                                     >  1:3588490/1‑93 (MQ=255)

                    gAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAca                                                                    >  2:1005965/1‑93 (MQ=255)

                    gAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAca                                                                    >  1:2920624/1‑93 (MQ=255)

                    gAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAca                                                                    >  1:2239005/1‑93 (MQ=255)

                     aaCGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAcat                                                                   <  1:700400/93‑1 (MQ=255)

                      aCGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAcatc                                                                  >  2:5052543/1‑93 (MQ=255)

                       cGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAcatca                                                                 <  2:2685307/93‑1 (MQ=255)

                       cGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAcatca                                                                 >  2:1009990/1‑93 (MQ=255)

                        gTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTGCCTGAcatcat                                                                <  1:5333120/93‑1 (MQ=255)

                        gTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGAcatcat                                                                <  1:4624172/93‑1 (MQ=255)

                         tttATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACGCTCCATCCAGGGCTACCTGAcatcata                                                               >  2:1248563/1‑92 (MQ=255)

                         tttATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACatcatc                                                               <  2:3441018/93‑1 (MQ=255)

                         tttATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACatcatc                                                               >  2:2064621/1‑93 (MQ=255)

                         tttATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACatcatc                                                               >  1:5375858/1‑93 (MQ=255)

                          ttattaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACAtcatca                                                              >  1:4724445/1‑93 (MQ=255)

                          ttattaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACAtcatca                                                              >  1:5268129/1‑93 (MQ=255)

                           tattaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCaa                                                             <  1:6045735/93‑1 (MQ=255)

                            attaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAag                                                            >  2:2099907/1‑93 (MQ=255)

                             ttaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagc                                                           >  1:3936728/1‑93 (MQ=255)

                             ttaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagc                                                           >  2:3769201/1‑93 (MQ=255)

                             ttaGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagc                                                           >  1:4083579/1‑93 (MQ=255)

                              taGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagca                                                          >  1:2279456/1‑93 (MQ=255)

                              taGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagca                                                          >  2:1444464/1‑93 (MQ=255)

                                gAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagcagc                                                        >  1:4334279/1‑93 (MQ=255)

                                gAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagcagc                                                        >  1:4118863/1‑93 (MQ=255)

                                gAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAagcagc                                                        >  1:864952/1‑93 (MQ=255)

                                 aaaaaTAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGcc                                                       <  2:5303526/93‑1 (MQ=255)

                                  aaaaTAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCt                                                      >  1:26481/1‑93 (MQ=255)

                                  aaaaTAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCt                                                      >  1:4823303/1‑93 (MQ=255)

                                  aaaaTAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCt                                                      <  2:5514893/93‑1 (MQ=255)

                                    aataaataAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGa                                                    <  2:2319775/93‑1 (MQ=255)

                                     ataaataaAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGAt                                                   <  2:609107/93‑1 (MQ=255)

                                     ataaataaAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGAt                                                   >  1:127964/1‑93 (MQ=255)

                                      taaataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGAtt                                                  >  2:4424715/1‑93 (MQ=255)

                                       aaataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTc                                                 <  2:4269677/93‑1 (MQ=255)

                                       aaataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTc                                                 <  2:3164580/93‑1 (MQ=255)

                                        aataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCa                                                >  1:1944427/1‑93 (MQ=255)

                                        aataaaACCAGCGTTGGGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTAGCTGACATCATCAAGCAGTCTGATTCa                                                >  1:781802/1‑93 (MQ=255)

                                         ataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAg                                               >  2:5750077/1‑93 (MQ=255)

                                          taaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGa                                              >  1:5625649/1‑93 (MQ=255)

                                          taaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGa                                              >  2:2018541/1‑93 (MQ=255)

                                          taaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATAATCAAGCAGCCTGATTCAGa                                              >  1:4680837/1‑93 (MQ=255)

                                            aaaCCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGaaa                                            <  2:273369/93‑1 (MQ=255)

                                               ccAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTc                                         >  2:6014660/1‑93 (MQ=255)

                                                 aGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCaa                                       <  2:1871432/93‑1 (MQ=255)

                                                  gCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCaaa                                      >  1:1604475/1‑93 (MQ=255)

                                                     ttGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATg                                   >  1:662432/1‑93 (MQ=255)

                                                     ttGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATg                                   >  1:4713147/1‑93 (MQ=255)

                                                     ttGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATg                                   >  1:5539355/1‑93 (MQ=255)

                                                     ttGTGTTTGAAACTCACGGGAATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTAAAAATg                                   >  1:26153/1‑93 (MQ=255)

                                                      tgtgTTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGt                                  >  2:5086437/1‑93 (MQ=255)

                                                       gtgtTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTa                                 <  2:2177960/93‑1 (MQ=255)

                                                       gtgtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTa                                 <  2:3658713/93‑1 (MQ=255)

                                                       gtgtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTa                                 >  2:1636752/1‑93 (MQ=255)

                                                       gtgtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTa                                 <  2:1728904/93‑1 (MQ=255)

                                                        tgtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTAc                                >  2:2586047/1‑93 (MQ=255)

                                                        tgtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTAc                                <  1:385105/93‑1 (MQ=255)

                                                         gtTTGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACt                               >  2:4103693/1‑93 (MQ=255)

                                                          tttGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTc                              >  1:4778435/1‑93 (MQ=255)

                                                          tttGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTc                              >  2:4817571/1‑93 (MQ=255)

                                                           ttGAAACTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCg                             <  2:2851455/93‑1 (MQ=255)

                                                              aaaCTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTg                          >  2:3032603/1‑93 (MQ=255)

                                                               aaCTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGa                         <  2:1604050/93‑1 (MQ=255)

                                                                aCTCACGGGCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGat                        <  1:3993225/93‑1 (MQ=255)

                                                                    gcGGGCATTACGAAAAAGTATCCATCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattata                    <  1:3528522/92‑1 (MQ=255)

                                                                       ggCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattatatta                 >  1:2990496/1‑93 (MQ=255)

                                                                       ggCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattatatta                 >  2:3152135/1‑93 (MQ=255)

                                                                        gCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattatattat                >  1:1192851/1‑93 (MQ=255)

                                                                        gCATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattatattat                <  2:51260/93‑1 (MQ=255)

                                                                         cATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATc               <  2:4512147/93‑1 (MQ=255)

                                                                         cATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATc               <  2:4914062/93‑1 (MQ=255)

                                                                         cATTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATc               >  2:5930460/1‑93 (MQ=255)

                                                                           ttACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGACAAAATGTACTCGCTGATTATAATATAGc             >  2:4138646/1‑93 (MQ=255)

                                                                             aCGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCa           <  1:5074910/93‑1 (MQ=255)

                                                                              cGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCAc          <  2:2971816/93‑1 (MQ=255)

                                                                               gAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACg         >  1:1142032/1‑93 (MQ=255)

                                                                                aaaaaGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGc        >  2:3529298/1‑93 (MQ=255)

                                                                                aaaaaGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGc        >  1:2040761/1‑93 (MQ=255)

                                                                                   aaGTACCCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCta     <  2:584647/93‑1 (MQ=255)

                                                                                   aaGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCta     >  1:2427856/1‑93 (MQ=255)

                                                                                    aGTAACCCTCCAGCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCtat    <  1:4915806/93‑1 (MQ=255)

                                                                                     gTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCtata   >  1:5714955/1‑93 (MQ=255)

                                                                                      tAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATAc  <  1:1495589/93‑1 (MQ=255)

AGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATAC  >  NC_000913/3453580‑3453758

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