BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	376,777	A→G	N272N (AAT→AAC)	yaiM ←	putative esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	376,777	0	A	G	99.5%	675.6 / ‑5.3	212	N272N (AAT→AAC)	yaiM	putative esterase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (92/119); total (92/120)
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.95e-01

Only 100 of 219 total aligned reads displayed.

CAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/376687‑376866

cAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGtt                                                                                         <  1:4003489/93‑1 (MQ=255)

 aGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTttgttg                                                                                        >  2:2866856/1‑93 (MQ=255)

 aGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTAACTttgttg                                                                                        >  1:256433/1‑93 (MQ=18)

  ggCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTgg                                                                                       >  1:3123183/1‑93 (MQ=255)

  ggCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTgg                                                                                       <  2:1651934/93‑1 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGt                                                                                    >  1:3708145/1‑93 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGt                                                                                    >  1:3470000/1‑93 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGt                                                                                    >  1:2125738/1‑93 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTg                                                                                   <  2:1871449/93‑1 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGGTGGCGTg                                                                                   >  2:2591492/1‑93 (MQ=255)

        cGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGt                                                                                 <  1:2774758/93‑1 (MQ=255)

         ggCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTa                                                                                <  1:1884921/93‑1 (MQ=255)

         ggCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTa                                                                                >  2:465923/1‑93 (MQ=255)

          gCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAg                                                                               >  1:1749347/1‑93 (MQ=255)

            ttgcGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGc                                                                             <  2:2737219/90‑1 (MQ=255)

            ttACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGc                                                                             >  2:66367/1‑93 (MQ=255)

             tACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCa                                                                            >  1:1715346/1‑93 (MQ=255)

              aCGGGCGAAGTGTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCaa                                                                           <  2:2767001/93‑1 (MQ=255)

                 ggCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAatat                                                                        <  2:877356/93‑1 (MQ=255)

                  gCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATg                                                                       >  1:2537333/1‑93 (MQ=255)

                   cGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGc                                                                      <  2:1582884/93‑1 (MQ=255)

                   cGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGc                                                                      <  2:1945641/93‑1 (MQ=255)

                    gAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCt                                                                     >  1:1728111/1‑93 (MQ=255)

                         gtAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCa                                                                >  1:1942005/1‑93 (MQ=255)

                          tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCaa                                                               >  1:1039452/1‑93 (MQ=255)

                              ccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATaaa                                                           <  2:1904090/93‑1 (MQ=255)

                              ccGGAGAAGGCGTTCACGCGGCAGCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTGGCCGTGGTAGGCAATATGCTCGCCaaaaaa                                                           <  2:2940503/93‑1 (MQ=255)

                                ggATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAAc                                                         <  1:4144236/93‑1 (MQ=255)

                                ggATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAAc                                                         >  1:3928123/1‑93 (MQ=255)

                                ggATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAAc                                                         >  1:4409488/1‑93 (MQ=255)

                                 gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACt                                                        <  1:2583504/93‑1 (MQ=255)

                                 gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACt                                                        <  1:3384617/93‑1 (MQ=255)

                                   tAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTgg                                                      >  2:1244472/1‑93 (MQ=255)

                                    aaGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                     <  2:3610144/93‑1 (MQ=255)

                                    aaGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                     <  1:2235164/93‑1 (MQ=255)

                                    aaGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                     <  1:3039680/93‑1 (MQ=255)

                                    aaGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                     <  2:177034/93‑1 (MQ=255)

                                       gCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGAc                                                  <  1:3569963/93‑1 (MQ=255)

                                       gCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGAc                                                  <  1:704251/93‑1 (MQ=255)

                                       gCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTACGCAATATGCTCGCCAATAAAACTGGAGAc                                                  >  1:147240/1‑93 (MQ=255)

                                        ggTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACa                                                 <  2:1610534/92‑1 (MQ=255)

                                        cGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACa                                                 >  1:1324164/1‑93 (MQ=255)

                                         gtccACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCCCCAATAAAACTGGAGACaa                                                <  2:1473924/90‑1 (MQ=255)

                                         gTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaa                                                <  2:959704/93‑1 (MQ=255)

                                           tCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaaa                                              <  2:4051710/93‑1 (MQ=255)

                                           tCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaaa                                              >  1:578994/1‑93 (MQ=255)

                                           tCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaaa                                              <  1:2340860/93‑1 (MQ=255)

                                            cACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaat                                             <  1:1652212/93‑1 (MQ=255)

                                             aCGCCCCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaata                                            <  2:1835959/93‑1 (MQ=255)

                                              cgccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaataa                                           <  1:3872785/93‑1 (MQ=255)

                                               gccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaataat                                          <  1:3806491/93‑1 (MQ=255)

                                                ccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                         <  1:2661106/93‑1 (MQ=255)

                                                ccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                         <  1:226108/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  2:70568/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  2:1711786/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  1:75923/93‑1 (MQ=255)

                                                   cATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCt                                      <  1:3200787/93‑1 (MQ=255)

                                                     tCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgt                                    >  1:832215/1‑93 (MQ=255)

                                                     tCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgt                                    <  2:1045940/93‑1 (MQ=255)

                                                      ccGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgtg                                   >  2:1224432/1‑93 (MQ=255)

                                                      ccGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgtg                                   >  1:2804109/1‑93 (MQ=255)

                                                       cGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGa                                  >  2:3665093/1‑93 (MQ=255)

                                                         gCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATc                                >  1:4657647/1‑93 (MQ=255)

                                                          cAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCa                               >  2:2500852/1‑93 (MQ=255)

                                                          cAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGCGATCa                               <  1:1031502/93‑1 (MQ=255)

                                                           aGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCAt                              <  1:2510462/93‑1 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:594585/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  2:1270385/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:2724544/1‑93 (MQ=255)

                                                              cGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAc                           <  1:1072979/93‑1 (MQ=255)

                                                              cGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAc                           <  1:781170/93‑1 (MQ=255)

                                                               gTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTTGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAcc                          >  2:2419213/1‑93 (MQ=255)

                                                               gTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGGCAAAATAATAGCTGTGATCATAAcc                          >  2:2749188/1‑93 (MQ=255)

                                                               gTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAcc                          >  1:3955532/1‑93 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  2:1074530/1‑93 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         <  2:4646954/93‑1 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  2:3309089/1‑93 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  1:4632166/1‑93 (MQ=255)

                                                                 gCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                        >  2:1159163/1‑93 (MQ=255)

                                                                 gCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                        >  1:877763/1‑93 (MQ=255)

                                                                  cTATTATCAACACATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTc                       <  1:3177471/93‑1 (MQ=255)

                                                                   tattatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCg                      <  1:4356649/93‑1 (MQ=255)

                                                                    attatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGt                     <  1:3096456/93‑1 (MQ=255)

                                                                      tatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGa                   <  1:2068068/93‑1 (MQ=255)

                                                                       atCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGAt                  >  1:117800/1‑93 (MQ=255)

                                                                          aaCGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAAc               <  1:3504615/93‑1 (MQ=255)

                                                                           aCGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACg              <  2:3759431/93‑1 (MQ=255)

                                                                            cGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACgg             >  2:4393038/1‑93 (MQ=255)

                                                                             gCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGa            <  2:2249879/93‑1 (MQ=255)

                                                                             gCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGa            <  2:4549452/93‑1 (MQ=255)

                                                                              cATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAt           <  1:2052717/93‑1 (MQ=255)

                                                                                tatTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTa         <  1:3003610/93‑1 (MQ=255)

                                                                                tatTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTa         <  1:3227228/93‑1 (MQ=255)

                                                                                 atTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa        >  1:2909179/1‑93 (MQ=255)

                                                                                  ttCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       <  2:4592005/93‑1 (MQ=255)

                                                                                  ttCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       >  2:581395/1‑93 (MQ=255)

                                                                                  ttCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       <  1:4306206/93‑1 (MQ=255)

                                                                                     aGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc    >  2:3003295/1‑93 (MQ=255)

                                                                                       tttGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCtt  >  2:1260888/1‑93 (MQ=255)

                                                                                       tttGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCtt  <  2:2127659/93‑1 (MQ=255)

CAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/376687‑376866

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