BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,224	A→G	T157T (ACT→ACC)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,224	0	A	G	99.8%	1705.1 / NA	530	T157T (ACT→ACC)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (328/200); total (329/201)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 605 total aligned reads displayed.

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

gCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGt                                                                                         <  1:3118531/93‑1 (MQ=255)

   aTAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTg                                                                                      >  1:7107914/1‑93 (MQ=255)

   aTAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTg                                                                                      >  1:3243673/1‑93 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   >  2:9545631/1‑93 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   >  1:8693513/1‑93 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   >  1:91524/1‑93 (MQ=255)

       agttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  <  2:6831595/91‑1 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTTATaa                                                                                  >  1:3966886/1‑93 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  >  1:7857635/1‑93 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  <  2:7985102/93‑1 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 >  1:8018585/1‑93 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 >  1:8289692/1‑93 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 <  2:1151591/93‑1 (MQ=255)

         ttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATg                                                                                <  1:2703181/93‑1 (MQ=255)

           cACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCt                                                                              >  2:4604418/1‑93 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCtt                                                                             <  2:4860625/93‑1 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTTGGTTTGATAATGCTTCt                                                                           >  1:8989296/1‑93 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCt                                                                           <  2:9278461/93‑1 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCt                                                                           <  1:8379033/93‑1 (MQ=255)

               gCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGACTTTGAAAACGTTCAGGTTGGGTTTGATAAGGCTtgtg                                                                          >  1:3615770/1‑93 (MQ=255)

               gCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTg                                                                          >  2:2535279/1‑93 (MQ=255)

               gCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTg                                                                          >  1:4698513/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAACCGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      <  2:7142463/93‑1 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  1:1549634/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  1:8161913/1‑93 (MQ=255)

                    ttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAAt                                                                     <  1:9834379/93‑1 (MQ=255)

                     tGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATc                                                                    >  1:1156875/1‑93 (MQ=255)

                      gggATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCg                                                                   <  2:9037635/93‑1 (MQ=255)

                      gggATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCg                                                                   >  2:5721734/1‑93 (MQ=255)

                         aTTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCg                                                                >  2:6139449/1‑93 (MQ=255)

                           ttGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAt                                                              <  1:4837257/93‑1 (MQ=255)

                            tGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAtt                                                             <  1:5588985/93‑1 (MQ=255)

                            tGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAtt                                                             >  2:3866519/1‑93 (MQ=255)

                             gCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAttt                                                            >  2:6459450/1‑93 (MQ=255)

                              cTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTg                                                           >  1:9599981/1‑93 (MQ=255)

                               tCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGc                                                          >  1:4104992/1‑93 (MQ=255)

                               tCACTGTATTTTCTGCCCGGATAAATACCGGGAGGAAGAATTTAAAAACGTTCAGGTTTGGTTTGATAATGCTCTTGGAAATCGACGATTTGc                                                          >  2:768516/1‑93 (MQ=255)

                                cACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCt                                                         >  2:9903047/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        <  1:2562329/93‑1 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        <  1:5145203/93‑1 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       >  2:8854964/1‑93 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       >  2:8233139/1‑93 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   <  2:2318541/64‑1 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      <  1:5975908/93‑1 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      >  2:2170289/1‑93 (MQ=255)

                                       ttttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGt                                                  >  2:3440242/1‑93 (MQ=255)

                                       ttttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGt                                                  >  1:2226595/1‑93 (MQ=255)

                                        tttCTGTCCGGATAAATACCGGAAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTc                                                 <  1:5765279/93‑1 (MQ=255)

                                          tCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAg                                               >  1:6666091/1‑93 (MQ=255)

                                          tCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAg                                               >  2:7014400/1‑93 (MQ=255)

                                            tGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGca                                             <  1:3789990/93‑1 (MQ=255)

                                               ccGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatca                                          >  1:2881489/1‑93 (MQ=255)

                                                cGGATAAATAACGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatcat                                         <  1:5526164/93‑1 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  2:10028651/1‑93 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  2:7400777/1‑93 (MQ=255)

                                                   aTAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAg                                      <  1:430779/93‑1 (MQ=255)

                                                   aTAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAg                                      >  2:7877234/1‑93 (MQ=255)

                                                   aTAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAg                                      <  1:7408324/93‑1 (MQ=255)

                                                     aaaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGtt                                    <  1:1773375/93‑1 (MQ=255)

                                                     aaaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGtt                                    >  2:3263596/1‑93 (MQ=255)

                                                        tACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcg                                 >  2:1847953/1‑93 (MQ=255)

                                                        tACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcg                                 >  2:4916606/1‑93 (MQ=255)

                                                         aCCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgt                                >  1:10027715/1‑93 (MQ=255)

                                                          cggtgaggaAGAATTTGAAAACGGTCAGGTTCGGTGTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtc                               <  2:845209/89‑1 (MQ=255)

                                                          ccGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtc                               <  1:4373169/93‑1 (MQ=255)

                                                          ccGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtc                               <  2:9397975/93‑1 (MQ=255)

                                                           cGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtcg                              >  2:286018/1‑93 (MQ=255)

                                                            gggAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGa                             >  1:5364607/1‑93 (MQ=255)

                                                             ggaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            >  1:3340325/1‑93 (MQ=255)

                                                             ggaggaAGAATTTGAAAACGTTCAGGTTCGCTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            <  2:5164090/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           <  1:6510951/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           <  2:4090448/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           >  2:1268882/1‑93 (MQ=255)

                                                               aggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAttt                          >  2:9716017/1‑93 (MQ=255)

                                                               aggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAttt                          >  2:8643672/1‑93 (MQ=255)

                                                               aggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGAGCAGCATCATAGGTTTAGTCTAttt                          >  1:3445495/1‑93 (MQ=255)

                                                                ggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTa                         <  2:7297393/93‑1 (MQ=255)

                                                                 gaagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAc                        >  1:7749634/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  1:1372896/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  1:9376911/93‑1 (MQ=255)

                                                                     aaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTg                    <  2:1661690/93‑1 (MQ=255)

                                                                      aTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGc                   >  2:4738525/1‑93 (MQ=255)

                                                                       tttGAAAACGTTCAGGTTCTGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCt                  >  2:5377435/1‑93 (MQ=255)

                                                                        ttGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTc                 >  2:7545226/1‑93 (MQ=255)

                                                                         tGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCg                >  1:7674161/1‑93 (MQ=255)

                                                                          gAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGa               >  2:4516769/1‑93 (MQ=255)

                                                                          gAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGa               <  2:3712473/93‑1 (MQ=255)

                                                                           aaaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaa              <  1:2015056/93‑1 (MQ=255)

                                                                           aaaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaa              <  1:1651083/93‑1 (MQ=255)

                                                                           aaaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaa              >  1:6416865/1‑93 (MQ=255)

                                                                             aaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAg            >  1:4545408/1‑93 (MQ=255)

                                                                                gttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcg         >  2:1481746/1‑93 (MQ=255)

                                                                                  tcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgct       <  2:6007609/93‑1 (MQ=255)

                                                                                  tcaggttcggTTGGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgct       <  2:3098873/93‑1 (MQ=255)

                                                                                    aggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctcg     >  2:5605871/1‑93 (MQ=255)

                                                                                     ggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctcgc    >  1:2158011/1‑93 (MQ=255)

                                                                                     ggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctcgc    >  1:8517864/1‑93 (MQ=255)

                                                                                      gttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTcgctcgct   >  1:4986906/1‑93 (MQ=255)

                                                                                      gttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTcgctcgct   >  1:8817062/1‑93 (MQ=255)

                                                                                       ttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCgctcgctc  >  2:1044315/1‑93 (MQ=255)

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

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: `‑`