BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,445	T→C	N265S (AAC→AGC)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,445	0	T	C	99.8%	1413.3 / ‑5.6	446	N265S (AAC→AGC)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (232/213); total (233/213)
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.64e-01

Only 100 of 458 total aligned reads displayed.

TGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCC  >  NC_000913/3897355‑3897534

tGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGATAATCAGCGGGATAATATCGCtc                                                                                         >  2:4294508/1‑93 (MQ=255)

   tACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtcatc                                                                                      >  2:4457531/1‑93 (MQ=255)

    aCCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtcatca                                                                                     >  1:4964540/1‑93 (MQ=255)

    aCCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtcatca                                                                                     <  1:2369132/93‑1 (MQ=255)

      ccGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGt                                                                                   <  2:876661/93‑1 (MQ=255)

      ccGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGt                                                                                   >  2:4363676/1‑93 (MQ=255)

        gcattGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 <  1:6133616/90‑1 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 >  1:1189240/1‑93 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 >  2:4028183/1‑93 (MQ=255)

         ccTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcg                                                                                >  1:7900087/1‑93 (MQ=255)

           ttGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgcg                                                                              >  2:1290218/1‑93 (MQ=255)

           ttGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgcg                                                                              >  2:1859428/1‑93 (MQ=255)

           ttGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgcg                                                                              >  1:4769416/1‑93 (MQ=255)

           ttGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgcg                                                                              >  2:5182471/1‑93 (MQ=255)

            tgtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgcgt                                                                             >  1:7552533/1‑92 (MQ=255)

             gtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGAt                                                                            >  1:7714794/1‑93 (MQ=255)

              tgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATc                                                                           <  2:2116941/93‑1 (MQ=255)

               gAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCt                                                                          >  1:5111490/1‑93 (MQ=255)

                aaaGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTg                                                                         >  2:5113506/1‑93 (MQ=255)

                 aaGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAATGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGc                                                                        <  2:4675580/93‑1 (MQ=255)

                 aaGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGc                                                                        >  1:6486729/1‑93 (MQ=255)

                 aaGCCCTGCCAGCGCGCGGCCATCGACATCTTTACGGACATAAAAACGCTTGCTAGTCAGCGGGATATTATCGCTCATCAGTTCGCGATCTGc                                                                        <  2:7025503/93‑1 (MQ=255)

                      cTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCtt                                                                   <  1:2793696/93‑1 (MQ=255)

                       tGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCttt                                                                  >  1:5773944/1‑93 (MQ=255)

                        gCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTa                                                                 >  1:4422413/1‑93 (MQ=255)

                        gCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTa                                                                 <  2:1419533/93‑1 (MQ=255)

                         ccAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAg                                                                >  2:5905589/1‑93 (MQ=255)

                           aGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              <  2:5749555/93‑1 (MQ=255)

                           aGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              <  1:7851894/93‑1 (MQ=255)

                           aGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              >  1:5680581/1‑93 (MQ=255)

                            gcgcgcGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCAt                                                             >  1:2161115/1‑93 (MQ=255)

                            gcgcgcGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCAt                                                             <  2:3181184/93‑1 (MQ=255)

                                   gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGCTAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCgg                                                      >  1:3590544/1‑93 (MQ=255)

                                    ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGa                                                     <  2:269783/93‑1 (MQ=255)

                                    ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGa                                                     >  1:3308407/1‑93 (MQ=255)

                                    ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGa                                                     >  1:401020/1‑93 (MQ=255)

                                     cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGAt                                                    <  2:6554270/93‑1 (MQ=255)

                                     cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGAt                                                    <  1:372462/93‑1 (MQ=255)

                                     cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGAt                                                    <  1:4321279/93‑1 (MQ=255)

                                      aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATa                                                   >  1:2881731/1‑93 (MQ=255)

                                      aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATa                                                   >  1:7157178/1‑93 (MQ=255)

                                      aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATa                                                   >  2:5529855/1‑93 (MQ=255)

                                       tCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaa                                                  >  1:1383923/1‑93 (MQ=255)

                                        cGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaaa                                                 <  1:6264051/93‑1 (MQ=255)

                                         ggCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaaaa                                                >  2:2828002/1‑93 (MQ=255)

                                         ggCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaaaa                                                <  1:1062378/93‑1 (MQ=255)

                                         ggCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaaaa                                                <  1:3456362/93‑1 (MQ=255)

                                           cATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGa                                              <  2:5999420/93‑1 (MQ=255)

                                           cATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGa                                              <  2:3689765/93‑1 (MQ=255)

                                            aTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGAc                                             >  1:247401/1‑93 (MQ=255)

                                            aTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGAc                                             >  2:7608013/1‑93 (MQ=255)

                                            aTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGAc                                             >  2:1738319/1‑93 (MQ=255)

                                              cTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           <  1:7399913/93‑1 (MQ=255)

                                              cTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           >  2:4996587/1‑93 (MQ=255)

                                              cTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           <  1:1870966/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          <  1:1290218/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          >  1:5552181/1‑93 (MQ=255)

                                                ttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTa                                         <  1:668702/93‑1 (MQ=255)

                                                ttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTa                                         >  1:4049968/1‑93 (MQ=255)

                                                 tACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTAc                                        <  2:6699395/93‑1 (MQ=255)

                                                 tACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTAc                                        <  2:5268538/93‑1 (MQ=255)

                                                   cGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGt                                      <  1:422094/93‑1 (MQ=255)

                                                        aTTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTc                                 >  2:4864863/1‑93 (MQ=255)

                                                        aTTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTc                                 <  2:4196158/93‑1 (MQ=255)

                                                        aTTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTc                                 >  2:3273214/1‑93 (MQ=255)

                                                         ttAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTcc                                >  2:7318955/1‑93 (MQ=255)

                                                           aaaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTg                              >  2:638777/1‑93 (MQ=255)

                                                           aaaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTg                              <  1:1337521/93‑1 (MQ=255)

                                                            aaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGa                             <  2:4621267/93‑1 (MQ=255)

                                                             aaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAg                            >  2:6199861/1‑93 (MQ=255)

                                                             aaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAg                            <  1:7094681/93‑1 (MQ=255)

                                                               cGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGc                          >  1:1497976/1‑93 (MQ=255)

                                                               cGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGc                          >  2:2006098/1‑93 (MQ=255)

                                                                 cTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCAc                        >  2:3107752/1‑93 (MQ=255)

                                                                 cTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCAc                        <  1:7925800/93‑1 (MQ=255)

                                                                  ttGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACg                       >  2:5831567/1‑93 (MQ=255)

                                                                   tGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGa                      >  1:3841774/1‑93 (MQ=255)

                                                                   ggCTAATCAGCGGGATAATATCGCTCATCAGTTCACGATCTGCCTCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGa                      <  2:109673/92‑1 (MQ=255)

                                                                      tAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAAt                   <  2:7387463/93‑1 (MQ=255)

                                                                      tAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAAt                   >  2:3509092/1‑93 (MQ=255)

                                                                        aTCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCt                 <  2:7680243/93‑1 (MQ=255)

                                                                         tCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCtt                >  1:7278091/1‑93 (MQ=255)

                                                                          cAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTc               <  2:1493665/93‑1 (MQ=255)

                                                                          cAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTc               <  2:2991614/93‑1 (MQ=255)

                                                                          cAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTc               <  2:2550762/93‑1 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             <  2:6798848/93‑1 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             >  2:2356102/1‑93 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             <  2:7814705/93‑1 (MQ=255)

                                                                             cGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg            >  2:3643243/1‑93 (MQ=255)

                                                                              gggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGa           >  1:7484230/1‑93 (MQ=255)

                                                                               ggATAATATCTCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          <  1:4526682/93‑1 (MQ=255)

                                                                               ggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          <  1:2008377/93‑1 (MQ=255)

                                                                               ggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          >  1:4652846/1‑93 (MQ=255)

                                                                                gATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaat         >  1:583858/1‑93 (MQ=255)

                                                                                 ataataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaata        >  1:5903520/1‑93 (MQ=255)

                                                                                  taataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaataa       >  1:5394781/1‑93 (MQ=255)

                                                                                   aataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaataat      >  1:3060964/1‑93 (MQ=255)

                                                                                   aataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTAACGTCTTTCCTGAGGCACGAAATCTTCGGGaataat      >  1:661969/1‑93 (MQ=255)

                                                                                     tatCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGc    >  1:5237668/1‑93 (MQ=255)

                                                                                      atCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGcc   <  1:6799553/93‑1 (MQ=255)

TGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCC  >  NC_000913/3897355‑3897534

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