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.3%	471.4 / ‑4.3	150	coding (86/903 nt)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (62/87); total (63/87)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 177 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑g                                                                                         <  2:1301765/93‑1 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gc                                                                                        <  1:126629/93‑2 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctg                                                                                      >  2:3686481/1‑90 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctg                                                                                      <  1:2070522/93‑4 (MQ=255)

   tcatcaCCACCACCAT‑AAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctgg                                                                                     <  1:3812565/93‑5 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctgg                                                                                     <  2:2840594/93‑5 (MQ=255)

      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGgcg                                                                                   <  2:4622384/93‑1 (MQ=255)

       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGgcgc                                                                                  <  2:4098809/93‑1 (MQ=255)

          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtat                                                                               <  2:2299046/93‑1 (MQ=255)

          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtat                                                                               <  2:3957079/93‑1 (MQ=255)

          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtat                                                                               <  2:1187715/93‑1 (MQ=255)

          caccacCAGAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtat                                                                               <  2:3957080/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCCACTTGAAGTAC‑GCTGGCGCtata                                                                              >  2:1941652/1‑93 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtata                                                                              >  2:2987805/1‑93 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtata                                                                              >  1:2882471/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAc                                                                             <  1:42768/93‑1 (MQ=255)

             caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAcc                                                                            >  2:3779127/1‑93 (MQ=255)

              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAccc                                                                           >  1:4579231/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCa                                                                         >  1:2569841/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCa                                                                         >  1:2709173/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCa                                                                         >  1:1581637/1‑93 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAccccca                                                                        >  1:1793869/1‑91 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaa                                                                       <  2:3836976/93‑1 (MQ=255)

                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataa                                                                    <  1:4586807/93‑1 (MQ=255)

                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataa                                                                    >  2:4275400/1‑93 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaa                                                                   >  1:3860144/1‑93 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGAAGTAAATCAACTTGAAGAACAC‑TGGCGCTATTCCCCaaataaa                                                                   >  2:3521479/1‑93 (MQ=255)

                       aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaa                                                                  <  2:4278788/93‑1 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaaa                                                                 >  1:1411802/1‑92 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaaa                                                                 >  2:1489427/1‑92 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTAA‑GCTGGCGCTATACCCCAaataaaataa                                                               >  1:1313183/1‑93 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAaataaaataa                                                               <  1:3040152/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAaataaaataa                                                               <  2:26064/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  2:3478645/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:3289447/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  2:3250573/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:1796482/93‑1 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAAAAAgag                                                            >  1:3794131/1‑91 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAA‑GCTGGCGCTATACCCCAAATAAAAAAAGCggaa                                                         >  1:4078523/1‑91 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        <  1:2895138/93‑1 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        <  1:1140273/93‑1 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                       >  1:1647920/1‑93 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGAAGAAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAAAAAGCGGCAag                                                       >  1:3527994/1‑93 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  2:900289/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:2837478/93‑1 (MQ=255)

                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                    <  1:3001346/93‑1 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATATCCCAAATATAATAAGCGGGAAGAGCg                                                   >  2:4213567/1‑93 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCg                                                   >  1:1180386/1‑93 (MQ=255)

                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  1:3906037/93‑1 (MQ=255)

                                        gcgATGCGTCAGGATTTCATCTGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                 <  2:2511539/93‑1 (MQ=255)

                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCCAATAAAATAAGCGGCAAGAGCGAga                                                >  1:1000075/1‑92 (MQ=255)

                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                <  2:2554040/93‑1 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                               <  2:2125147/93‑1 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGAGTa                                               >  2:4437140/1‑93 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             >  2:4460521/1‑93 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             >  1:1848728/1‑93 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             <  2:1483143/93‑1 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             <  1:3315642/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:2390680/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           >  1:1935766/1‑93 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  2:464692/1‑93 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGG‑TGGCGCTATACCCCAAATAAAATAAGAAGGACAAGCGAGTAATTCGc                                         >  2:2015683/1‑93 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         >  1:2561678/1‑93 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         >  2:3647489/1‑93 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         >  2:1140031/1‑93 (MQ=255)

                                                 cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcc                                        <  1:2797986/93‑1 (MQ=255)

                                                 cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcc                                        <  2:3270590/93‑1 (MQ=255)

                                                 cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGGAATACGcc                                        >  1:3751287/1‑93 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:2431559/93‑1 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:241344/93‑1 (MQ=255)

                                                    gATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                                     >  2:693635/1‑93 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  <  1:4007263/93‑1 (MQ=255)

                                                        tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc                                 >  1:2280684/1‑93 (MQ=255)

                                                        tCATCGGCGGGCACGTAGGAAATCAACTTGAAGTAC‑GCTGGGGCTATACCCCAAATAAAATAACCGGGAAGAGCGACAAATACGCCCTGCcgg                                 >  2:1007421/1‑92 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  2:3902444/1‑93 (MQ=255)

                                                                gggCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                         <  2:346173/93‑1 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:4682397/1‑93 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     <  1:4250160/93‑1 (MQ=255)

                                                                      gtagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                   >  1:2969370/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  1:3949986/93‑1 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  >  2:2772439/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  >  1:594606/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  >  2:2215025/1‑93 (MQ=255)

                                                                        agtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                 <  1:4731507/93‑1 (MQ=255)

                                                                         gtaAATCAACTTGAAGTAC‑GCTGGCTCTATTCCCCAAATAAAATAAGCGGCAAGAGCGATTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  1:3221111/93‑1 (MQ=255)

                                                                         gtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  2:3975373/93‑1 (MQ=255)

                                                                          taAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               <  2:1420432/93‑1 (MQ=255)

                                                                           aaaTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg              <  2:2920289/93‑1 (MQ=255)

                                                                            aaTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  1:3729852/1‑93 (MQ=255)

                                                                               cAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa          <  1:3472294/93‑1 (MQ=255)

                                                                               cAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa          <  1:2075287/93‑1 (MQ=255)

                                                                               cAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa          <  2:2348308/93‑1 (MQ=255)

                                                                                aaCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         >  2:339575/1‑93 (MQ=255)

                                                                                aaCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         >  2:2750413/1‑93 (MQ=255)

                                                                                  cTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       <  2:1634682/93‑1 (MQ=255)

                                                                                     gAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    <  1:1968784/93‑1 (MQ=255)

                                                                                     gAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    >  1:4510625/1‑93 (MQ=255)

                                                                                      aaGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt   <  1:4444913/93‑1 (MQ=255)

                                                                                      aaGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt   >  1:2861259/1‑93 (MQ=255)

                                                                                       aGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt  <  1:3187192/93‑1 (MQ=255)

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