BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,717	1	.	C	99.8%	1593.5 / NA	476	T29R (ACG→AGG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (217/258); minor base A (0/1); total (217/259)
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.99e-01

Only 100 of 571 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                                         >  1:606216/1‑93 (MQ=255)

tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                                         >  2:1594497/1‑93 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACgc                                                                                        >  2:3928397/1‑92 (MQ=255)

  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  1:428961/93‑3 (MQ=255)

  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  1:7876240/93‑3 (MQ=255)

  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  2:8685270/93‑3 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTGAATCAACTTGAAGTACgctg                                                                                      >  2:7539991/1‑90 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctg                                                                                      <  2:2364543/93‑4 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctg                                                                                      >  1:407419/1‑90 (MQ=255)

      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                   >  2:7421010/1‑93 (MQ=255)

       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                                  <  1:6639602/93‑1 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                <  2:2758111/93‑1 (MQ=255)

         acaccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                <  1:4450915/92‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCtata                                                                              >  2:2045760/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             <  1:3781040/93‑1 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATAc                                                                             >  2:6797637/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTCCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             >  1:221215/1‑93 (MQ=255)

             cgccaTAAAGAAAAACGAGCGGATCACCAGATGCGTCAGGATTTCATCGGCCGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcc                                                                            <  2:1250771/91‑1 (MQ=255)

              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                                           <  1:1771126/93‑1 (MQ=255)

               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcccc                                                                          <  1:7121817/93‑1 (MQ=255)

               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcccc                                                                          >  1:1421074/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCa                                                                         >  2:4318249/1‑93 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaa                                                                        >  2:3977942/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                                       >  2:8514636/1‑93 (MQ=255)

                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGAACGTAGTAAAACAAATTGAAGTACGCTTGCGCTATACCCCaaata                                                                     >  2:4431569/1‑93 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaa                                                                   >  2:4242386/1‑93 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaat                                                                 >  1:5860784/1‑93 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                >  1:8009056/1‑93 (MQ=255)

                         aaaCGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                <  2:8753790/93‑1 (MQ=255)

                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                              <  2:179982/93‑1 (MQ=255)

                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                              <  2:3632217/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                             >  2:901046/1‑93 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:3338333/93‑1 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAATAAGCg                                                            >  2:3208707/1‑93 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                            <  2:5958166/93‑1 (MQ=255)

                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                           <  1:1999578/93‑1 (MQ=255)

                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGc                                                          <  2:3109539/93‑1 (MQ=255)

                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGACGCTGGCGCTATACCCCAAATAAAAAAAGCGGc                                                          >  1:8796430/1‑93 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                         <  1:4509/93‑1 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        >  1:5221800/1‑93 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        <  1:4710765/93‑1 (MQ=255)

                                 aGATCACGCGATGCGTCAGGAGTTCATCGGCGGGCAAGGGGTTAATCAACTTGAAGGAAGCTCGCGCTATAGCCCAAATAAAATAAGCGGCaa                                                        >  2:5833031/1‑93 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      >  1:5101189/1‑93 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                     >  2:5095906/1‑93 (MQ=255)

                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                    <  1:5480691/93‑1 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                   >  2:2680806/1‑93 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGACGCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAACg                                                   >  1:7741377/1‑93 (MQ=255)

                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  2:1409984/93‑1 (MQ=255)

                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                 >  2:8608108/1‑93 (MQ=255)

                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCAAATAAAATAAGCGGGAAGAGCGAg                                                 >  2:5301888/1‑93 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                               <  2:8326572/93‑1 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGGACGTAGTAAATCAACTTGAAGAAAGCTGGCGCCATACCCCAAATAAACTAAACGGCAAGAGCGAGGaa                                              >  2:1142604/1‑93 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                              <  1:8557309/93‑1 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                              >  1:2591056/1‑93 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:6559230/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  2:4903714/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:3892588/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            >  1:6961279/1‑93 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGGAAAACAACTTTAAGTACGCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGAGGAATa                                            >  2:7690825/1‑93 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  1:16210/93‑1 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:7210646/93‑1 (MQ=255)

                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                      >  1:7412154/1‑93 (MQ=255)

                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                    <  1:7650762/93‑1 (MQ=255)

                                                      tttCATCGGCGGGCACGTAGGAAATCAATTTTAAGGACGTTGTCGCTTTAGCCCAAATTGAATAAGCGGTCAGAGCGAGTGATACGACCTggc                                   >  2:6414182/1‑91 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  <  1:3629290/93‑1 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  <  1:4197907/93‑1 (MQ=255)

                                                        tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc                                 >  2:7765027/1‑93 (MQ=255)

                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               <  1:11848/93‑1 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  1:8307197/1‑93 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  2:384240/93‑1 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  1:6796228/1‑93 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  1:2999762/93‑1 (MQ=255)

                                                            cggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                             >  2:4555271/1‑93 (MQ=255)

                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            <  2:1892610/93‑1 (MQ=255)

                                                              gcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt                           <  2:343085/93‑1 (MQ=255)

                                                                gggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                         >  2:2598404/1‑93 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  1:977583/1‑93 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  1:8114732/1‑93 (MQ=255)

                                                                     cGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                    <  1:5425607/93‑1 (MQ=255)

                                                                     cGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                    >  2:6700510/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  2:7944002/93‑1 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  1:2803561/93‑1 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  2:2967884/93‑1 (MQ=255)

                                                                        agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                 >  2:6326451/1‑93 (MQ=255)

                                                                        agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                 <  1:3205665/93‑1 (MQ=255)

                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  2:5388937/93‑1 (MQ=255)

                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  1:2926325/93‑1 (MQ=255)

                                                                          taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGCAATACGCCCTGCCGCGTTTGTTTTGcatccatc               >  1:8493124/1‑93 (MQ=255)

                                                                            aaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  2:8136717/1‑93 (MQ=255)

                                                                            aaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             <  1:4983836/93‑1 (MQ=255)

                                                                              tCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa           <  2:7039735/93‑1 (MQ=255)

                                                                              tCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa           >  2:3317240/1‑93 (MQ=255)

                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         >  2:4695689/1‑93 (MQ=255)

                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  2:610621/93‑1 (MQ=255)

                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAAGAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  1:7245600/93‑1 (MQ=255)

                                                                                 aCTTGAAGTACGCTGGCGCTATTCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        <  1:1043402/93‑1 (MQ=255)

                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       <  1:6196670/93‑1 (MQ=255)

                                                                                   ttGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      >  2:6535669/1‑93 (MQ=255)

                                                                                   ttGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      <  2:8278048/93‑1 (MQ=255)

                                                                                   ttGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTATTGCATCAATCGGGAATACt      >  1:2432787/1‑93 (MQ=255)

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

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