BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,061,601	T→C	C194R (TGC→CGC)	torD →	TorD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,061,601	0	T	C	99.3%	447.3 / ‑5.0	143	C194R (TGC→CGC)	torD	TorD
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (53/88); total (55/88)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 147 total aligned reads displayed.

TGGTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCG  >  NC_000913/1061511‑1061690

tGGTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAgcgc                                                                                         >  2:679504/1‑93 (MQ=255)

  gTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGa                                                                                       >  1:4312294/1‑93 (MQ=255)

  gTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGa                                                                                       >  2:477836/1‑93 (MQ=255)

      ccAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGaccac                                                                                   >  1:2037015/1‑93 (MQ=255)

       cAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGaccacc                                                                                  >  2:4044955/1‑93 (MQ=255)

       cAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGaccacc                                                                                  >  1:2407375/1‑93 (MQ=255)

        agagTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGAccacca                                                                                 >  2:2466294/1‑93 (MQ=255)

        agagTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGAccacca                                                                                 >  1:1305938/1‑93 (MQ=255)

         gagTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCaa                                                                                <  1:3403886/93‑1 (MQ=255)

          agTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCaaa                                                                               >  2:3480284/1‑93 (MQ=255)

           gTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCaaaa                                                                              <  2:510628/93‑1 (MQ=255)

            tttGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAAc                                                                             >  1:515327/1‑93 (MQ=255)

            tttGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAAc                                                                             <  2:2924856/93‑1 (MQ=255)

            tttGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAAc                                                                             <  2:1526904/93‑1 (MQ=255)

             ttgttgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACa                                                                            >  1:192593/1‑93 (MQ=255)

                 tgCGCGTTGTCGTCAGTATGTCGGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATa                                                                        <  2:2650277/93‑1 (MQ=255)

                 tgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATa                                                                        <  2:3759113/93‑1 (MQ=255)

                  gcgcGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATaa                                                                       <  2:878514/93‑1 (MQ=255)

                  gcgcGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATaa                                                                       <  1:1749822/93‑1 (MQ=255)

                   cgcgTTGTCGTCCGTATGACAGCTTTGGTTTTTACGCGGCACTACGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAc                                                                      <  2:4442190/93‑1 (MQ=255)

                    gcgTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcg                                                                     >  1:626848/1‑93 (MQ=255)

                     cgTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgt                                                                    >  1:508838/1‑93 (MQ=255)

                       ttGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgtcg                                                                  <  1:4145809/93‑1 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                >  2:3412612/1‑93 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  1:3214171/93‑1 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                >  2:1425455/1‑93 (MQ=255)

                         ggcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  2:1429777/91‑1 (MQ=255)

                         ggcgtcAGTATGACAGCTTTGATTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  2:1531245/91‑1 (MQ=255)

                          tcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGttt                                                               <  2:597798/93‑1 (MQ=255)

                          tcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGttt                                                               <  2:1492442/93‑1 (MQ=255)

                          tcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGttt                                                               <  2:1479365/93‑1 (MQ=255)

                              cAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgc                                                           <  1:4648534/93‑1 (MQ=255)

                              cAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgc                                                           >  2:1299986/1‑93 (MQ=255)

                               aGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcg                                                          >  1:60610/1‑93 (MQ=255)

                               aGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcg                                                          >  1:4238056/1‑93 (MQ=255)

                                gTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcgc                                                         >  1:2688911/1‑93 (MQ=255)

                                   tGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTg                                                      >  1:2555616/1‑93 (MQ=255)

                                    gACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGa                                                     <  1:384034/93‑1 (MQ=255)

                                    gACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGa                                                     <  1:2944994/93‑1 (MQ=255)

                                       aGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaaaa                                                  >  2:1818133/1‑93 (MQ=255)

                                       aGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaaaa                                                  <  1:3549937/93‑1 (MQ=255)

                                         cTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGa                                                >  2:2173942/1‑93 (MQ=255)

                                         cTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGa                                                <  1:2727041/93‑1 (MQ=255)

                                          tttGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAc                                               <  2:4662697/93‑1 (MQ=255)

                                          tttGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAc                                               <  1:143643/93‑1 (MQ=255)

                                            tGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgc                                             <  2:243481/93‑1 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgcg                                            <  2:4342842/93‑1 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgcg                                            <  1:1031162/93‑1 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGT‑GAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGt                                           >  1:1122128/1‑93 (MQ=255)

                                               tttttACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGtt                                          <  2:3249786/93‑1 (MQ=255)

                                               tttttACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGtt                                          >  2:3105739/1‑93 (MQ=255)

                                                  ttACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAg                                       <  2:909155/93‑1 (MQ=255)

                                                   tACGCGGCACTAAGCCAGGTATTGCTGGTGTTAGTGGAGGGCGGCCCCCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGc                                      <  2:3941707/93‑1 (MQ=255)

                                                    aCGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCg                                     <  1:1087945/93‑1 (MQ=255)

                                                     cgcgGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGt                                    <  1:4435905/93‑1 (MQ=255)

                                                      gcgGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTc                                   <  1:168177/93‑1 (MQ=255)

                                                      gcgGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTc                                   <  1:4066671/93‑1 (MQ=255)

                                                        ggCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGc                                 <  2:4318131/93‑1 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAACGACGCGTTTAGCGTCGCATc                              >  2:1944642/1‑93 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATc                              >  2:4563324/1‑93 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATc                              >  1:408548/1‑93 (MQ=255)

                                                            cTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATcg                             <  1:3695962/93‑2 (MQ=255)

                                                            cTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCa                             >  1:949182/1‑93 (MQ=255)

                                                            cTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCa                             >  1:4397150/1‑93 (MQ=255)

                                                             tAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCATCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAg                            >  1:2064422/1‑93 (MQ=255)

                                                             tAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAg                            <  1:354349/93‑1 (MQ=255)

                                                             tAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAg                            >  2:3047823/1‑93 (MQ=255)

                                                               aGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGc                          >  2:1929322/1‑93 (MQ=255)

                                                                gCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCa                         <  1:1613235/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  1:1803522/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  1:1841831/93‑1 (MQ=255)

                                                                  cAGTTATTGCTGGTGTTAGTGTAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAtt                       <  1:3106808/93‑1 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      <  1:1747334/93‑1 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      <  2:1659127/93‑1 (MQ=255)

                                                                    gTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTAt                     >  2:2764699/1‑93 (MQ=255)

                                                                    gTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTAt                     >  2:196218/1‑93 (MQ=255)

                                                                    gTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTAt                     >  1:259454/1‑93 (MQ=255)

                                                                     ttATTGCTGGTGTTAGTGGAGCGCGACCACCAAAATAGATAACGTCGTTTGTGCGCCTGAGAAGACGCGTTTAGCGTCGCATCAGGCATTATg                    <  1:2738701/93‑1 (MQ=255)

                                                                     ttATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATg                    <  1:863592/93‑1 (MQ=255)

                                                                     ttATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATg                    <  1:2125322/93‑1 (MQ=255)

                                                                     ttATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATg                    <  1:2941765/93‑1 (MQ=255)

                                                                         tGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGgcgc                <  2:2348766/93‑1 (MQ=255)

                                                                          gCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCa               >  2:121317/1‑93 (MQ=255)

                                                                          gCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCa               >  1:4683679/1‑93 (MQ=255)

                                                                           cTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAg              <  2:3799566/93‑1 (MQ=255)

                                                                             ggTGTTAGTGGAGCGCGCCCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGtt            <  1:1969077/93‑1 (MQ=255)

                                                                             ggTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGtt            <  1:706275/93‑1 (MQ=255)

                                                                             ggTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGtt            <  2:3055172/93‑1 (MQ=255)

                                                                             ggTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGtg            >  1:3706269/1‑92 (MQ=255)

                                                                              gtgtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTg           <  1:2077313/93‑1 (MQ=255)

                                                                                tttAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:336075/92‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:3863817/93‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:462334/93‑1 (MQ=255)

                                                                                  tAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCgg       <  1:2193567/93‑1 (MQ=255)

                                                                                   aGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGa      <  1:4294013/93‑1 (MQ=255)

                                                                                   aGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGa      >  2:289300/1‑93 (MQ=255)

                                                                                    gTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGAt     <  2:4155334/93‑1 (MQ=255)

                                                                                    gTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGAt     <  1:3384262/93‑1 (MQ=255)

                                                                                      ggAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATgc   <  2:3885696/93‑1 (MQ=255)

                                                                                       gAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATgcg  <  2:3121886/93‑1 (MQ=255)

TGGTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCG  >  NC_000913/1061511‑1061690

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