BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,044,033	(C)_5→6	coding (1149/7152 nt)	b1978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,044,028	1	.	C	93.3%	632.4 / 36.3	208	coding (1144/7152 nt)	b1978	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/13); new base C (119/75); total (120/88)
Fisher's exact test for biased strand distribution p-value = 7.23e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 234 total aligned reads displayed.

TCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043936‑2044115

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       >  2:3381348/1‑93 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       >  2:73525/1‑93 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       >  1:321443/1‑93 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       >  2:1618521/1‑93 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       <  2:5666233/93‑1 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       <  1:1874187/93‑1 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cc                                                                                       >  2:4399159/1‑93 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cccccct                                                                                  >  1:5739769/1‑91 (MQ=255)

        cAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑cccccctt                                                                                 >  2:1837779/1‑90 (MQ=255)

         aGTATTATGGCGATGAAGTGGCTCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑ccccccttc                                                                                <  1:4352865/93‑5 (MQ=255)

           tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTccc                                                                               >  1:1537072/1‑93 (MQ=255)

           tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTccc                                                                               >  1:5628214/1‑93 (MQ=255)

           tataatGGCGATGAAGTGGCCCTGATAGATACAGACGTACGGCACAGTTATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTccc                                                                               <  1:6307652/93‑1 (MQ=255)

            attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCg                                                                              >  1:5199015/1‑93 (MQ=255)

             ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGc                                                                             <  2:4009419/93‑1 (MQ=255)

              tatGGGGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCt                                                                            <  2:5239179/93‑1 (MQ=255)

                 ggCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCtgat                                                                         >  2:3474098/1‑93 (MQ=255)

                  gCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCtgatg                                                                        >  2:3030546/1‑93 (MQ=255)

                    gATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGAc                                                                      >  2:326516/1‑93 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  1:4682969/1‑93 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  1:2937784/1‑93 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  1:4565450/1‑93 (MQ=255)

                       gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCtt                                                                   >  1:291274/1‑93 (MQ=255)

                       gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCtt                                                                   >  1:3703404/1‑93 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  >  2:1405685/1‑93 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  <  2:3918924/93‑1 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  <  1:6132321/93‑1 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCAGTAGAGTAATCCTCATGCCATAACCGCTGGACTTAACTATT‑C‑CCCCCTTCCCGCTGATGACCTTc                                                                  <  2:5795987/93‑1 (MQ=255)

                         aGTGGCACTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAATCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCa                                                                 <  1:3379995/93‑1 (MQ=255)

                          gTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAg                                                                >  2:3456308/1‑93 (MQ=255)

                            ggCCCTGTTCGCTAAAGCAGATTGGAAAAGTAATCTTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAgcg                                                              <  1:2036867/93‑1 (MQ=255)

                            ggCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAgcg                                                              >  1:6194294/1‑93 (MQ=255)

                             gCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAgcgc                                                             <  1:1041526/93‑1 (MQ=255)

                               ccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCgg                                                           >  1:2338017/1‑93 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          <  2:1261600/93‑1 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          >  2:316154/1‑93 (MQ=255)

                                 tGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAg                                                         >  2:2739706/1‑93 (MQ=255)

                                 tGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAg                                                         <  2:5035277/93‑1 (MQ=255)

                                 tGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAg                                                         >  1:2522265/1‑93 (MQ=255)

                                  gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGc                                                        >  1:3281438/1‑93 (MQ=255)

                                   ttCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCa                                                       >  2:4876/1‑93 (MQ=255)

                                   ttCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCa                                                       >  1:5312380/1‑93 (MQ=255)

                                    tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCaa                                                      <  1:4790092/93‑1 (MQ=255)

                                    tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCaa                                                      >  2:4788224/1‑93 (MQ=255)

                                    tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCaa                                                      <  1:1741101/93‑1 (MQ=255)

                                    tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGAGCGGAGCaa                                                      >  1:5916138/1‑93 (MQ=255)

                                     cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAAc                                                     <  1:3710623/93‑1 (MQ=255)

                                       aTAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGc                                                   <  1:4780524/93‑1 (MQ=255)

                                        tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGcc                                                  >  1:264967/1‑93 (MQ=255)

                                          aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                >  1:2094745/1‑93 (MQ=255)

                                             ccactCGCCAAAGTACTCCTCATCCCATACCCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGt                                             <  2:556018/89‑1 (MQ=255)

                                              cGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTa                                            <  1:1457348/93‑1 (MQ=255)

                                              cGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTa                                            <  1:1285373/93‑1 (MQ=255)

                                               gATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTaa                                           <  1:1373921/93‑1 (MQ=255)

                                                 tCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑TCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAAc                                         <  1:1146155/93‑1 (MQ=255)

                                                  cGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACa                                        <  2:2047263/93‑1 (MQ=255)

                                                   ggCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       >  1:3979160/1‑93 (MQ=255)

                                                   ggCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       >  1:1583497/1‑93 (MQ=255)

                                                   ggCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       >  2:2301038/1‑93 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:2968798/93‑1 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:5297514/93‑1 (MQ=255)

                                                      aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑TCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGc                                    <  2:3430912/93‑1 (MQ=255)

                                                      aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGc                                    >  1:6238733/1‑93 (MQ=255)

                                                       aaGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCg                                   <  1:5438434/93‑1 (MQ=255)

                                                           aaTCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaaa                               >  2:2139404/1‑93 (MQ=255)

                                                            aTCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAAt                              >  1:4967066/1‑93 (MQ=255)

                                                             tCCTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGATCATCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATg                             >  2:6289994/1‑93 (MQ=255)

                                                               cTCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGCAAATGac                           >  1:4487139/1‑93 (MQ=255)

                                                                tCATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGaca                          >  1:4068620/1‑93 (MQ=255)

                                                                 cATGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGacac                         <  1:2255764/93‑1 (MQ=255)

                                                                  aTGCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACAcc                        >  1:4895190/1‑93 (MQ=255)

                                                                    gCCATAACCGCTGGACTTAACTA‑TCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAACGACACCCg                      <  2:6152851/93‑1 (MQ=255)

                                                                    gCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCg                      >  2:403998/1‑93 (MQ=255)

                                                                    gCCATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCg                      >  1:3802409/1‑93 (MQ=255)

                                                                     ccATAACCGCTGGACTTAACTATT‑C‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGt                     <  1:2052572/93‑1 (MQ=255)

                                                                      cATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCgat                    >  2:1160753/1‑91 (MQ=255)

                                                                      cATAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtt                    <  1:850496/93‑1 (MQ=255)

                                                                        tAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  <  1:985865/93‑1 (MQ=255)

                                                                        tAACCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  >  2:1509008/1‑93 (MQ=255)

                                                                         aaCCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTg                 <  2:2641462/93‑1 (MQ=255)

                                                                          aCCGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                >  1:2873047/1‑93 (MQ=255)

                                                                           ccGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               <  2:813526/93‑1 (MQ=255)

                                                                           ccGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               >  2:681571/1‑93 (MQ=255)

                                                                           ccGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               >  2:4378624/1‑93 (MQ=255)

                                                                           ccGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               <  1:126012/93‑1 (MQ=255)

                                                                           ccGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGcc               >  1:1496983/1‑93 (MQ=255)

                                                                            cGCTGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCg              >  1:2957284/1‑93 (MQ=255)

                                                                               tGGACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTAccgccg           >  1:4178072/1‑93 (MQ=255)

                                                                                ggACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGa          >  2:4651771/1‑93 (MQ=255)

                                                                                ggACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGa          >  2:4247935/1‑93 (MQ=255)

                                                                                ggACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGa          >  1:1220519/1‑93 (MQ=255)

                                                                                 gACTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAt         >  2:1332834/1‑93 (MQ=255)

                                                                                  aCTTAACTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtt        >  2:5151229/1‑93 (MQ=255)

                                                                                      accTA‑T‑ACCCCCCTTCCCGCTGATGGCCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACCCCCGTTTTGCCGTCGATTTTac    <  2:311799/91‑1 (MQ=255)

                                                                                      aaCTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTac    <  2:2896059/93‑1 (MQ=255)

                                                                                      aaCTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTac    <  2:5886048/93‑1 (MQ=255)

                                                                                      aaCTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTac    >  2:601249/1‑93 (MQ=255)

                                                                                       aCTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacc   <  1:5680446/93‑1 (MQ=255)

                                                                                       aCTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacc   >  2:3768635/1‑93 (MQ=255)

                                                                                        cTA‑T‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacct  <  2:5821732/93‑1 (MQ=255)

TCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTA‑T‑A‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043936‑2044115

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