BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,047,558	T→C	V4V (GTA→GTG)	gcvH ←	GcvH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,047,558	0	T	C	97.9%	757.4 / 1.8	242	V4V (GTA→GTG)	gcvH	GcvH
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (125/112); total (127/115)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 249 total aligned reads displayed.

CTGAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAACCGGC  >  NC_000913/3047468‑3047647

cTGAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCAc                                                                                         >  2:1620626/1‑93 (MQ=255)

cTGAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCAc                                                                                         >  1:1790435/1‑93 (MQ=255)

cTGAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCAc                                                                                         <  1:1755401/93‑1 (MQ=255)

  gAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGt                                                                                       <  1:3528977/93‑1 (MQ=255)

  gAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGt                                                                                       >  1:45802/1‑93 (MQ=255)

  gAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGt                                                                                       <  2:2094126/93‑1 (MQ=255)

  gAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGt                                                                                       <  2:660100/93‑1 (MQ=255)

   aGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGtt                                                                                      >  2:3738937/1‑93 (MQ=255)

       tGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtc                                                                                  >  2:2471705/1‑93 (MQ=255)

       tGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtc                                                                                  <  2:1041808/93‑1 (MQ=255)

       tGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtc                                                                                  >  2:945894/1‑93 (MQ=255)

         ttCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtcat                                                                                >  2:4379978/1‑93 (MQ=255)

         ttCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtcat                                                                                <  2:793816/93‑1 (MQ=255)

         ttCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCtcat                                                                                <  1:971671/93‑1 (MQ=255)

            ggTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCatcaa                                                                             <  2:4140903/93‑1 (MQ=255)

            ggTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCatcaa                                                                             >  1:2172375/1‑93 (MQ=255)

             gTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCatcaat                                                                            >  1:1337648/1‑93 (MQ=255)

              tAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCatcaatc                                                                           <  2:1088300/93‑1 (MQ=255)

               aaTACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCatcaatca                                                                          >  1:783608/1‑93 (MQ=255)

                aTACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCAtcaatcaa                                                                         <  1:3912275/93‑1 (MQ=255)

                  aCCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCaatcaatc                                                                       <  2:2922720/93‑1 (MQ=255)

                  aCCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCaatcaatc                                                                       <  1:2341172/93‑1 (MQ=255)

                   ccAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAAtct                                                                      <  1:3126158/93‑1 (MQ=255)

                    cAACGGTGTAAGTGCCGTCAGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAAtctc                                                                     >  1:2468220/1‑93 (MQ=255)

                     aaCGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTcc                                                                    >  2:3699110/1‑93 (MQ=255)

                       cGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaa                                                                  >  2:4014061/1‑93 (MQ=255)

                        ggTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaaa                                                                 >  2:3261113/1‑93 (MQ=255)

                        ggTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaaa                                                                 >  1:3023829/1‑93 (MQ=255)

                         gtgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaaaa                                                                >  2:2237632/1‑93 (MQ=255)

                         gtgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaaaa                                                                >  2:2691664/1‑93 (MQ=255)

                         gtgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCaaaa                                                                >  2:354404/1‑93 (MQ=255)

                            tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAg                                                             <  1:405072/93‑1 (MQ=255)

                             aaGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGt                                                            <  2:1648218/93‑1 (MQ=255)

                              aGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                                           >  1:3041345/1‑93 (MQ=255)

                               gTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTaa                                                          <  2:1239429/93‑1 (MQ=255)

                                ttccGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTaaa                                                         <  2:1359253/91‑1 (MQ=255)

                                 gCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAAt                                                        <  1:3966752/93‑1 (MQ=255)

                                  ccGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCATTCAATCTCCAAAAAAGTAAATc                                                       >  1:3310206/1‑93 (MQ=255)

                                  ccGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATc                                                       >  1:1658987/1‑93 (MQ=255)

                                   cGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                                      <  1:4231896/93‑1 (MQ=255)

                                    gTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAc                                                     >  2:254586/1‑93 (MQ=255)

                                    gTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAc                                                     <  2:4236689/93‑1 (MQ=255)

                                    gTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAc                                                     <  2:4722938/93‑1 (MQ=255)

                                    gTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAc                                                     >  1:1935873/1‑93 (MQ=255)

                                       ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgcg                                                  >  2:1715759/1‑93 (MQ=255)

                                        gCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGa                                                 >  2:4057298/1‑93 (MQ=255)

                                        gCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGa                                                 >  1:1454891/1‑93 (MQ=255)

                                        gCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGa                                                 >  2:694513/1‑93 (MQ=255)

                                          ttcttTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACg                                               >  1:4176574/1‑93 (MQ=255)

                                           tcttTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                                              >  1:2792302/1‑93 (MQ=255)

                                            cttTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCTATCAATCTCCAAAAAAGTAAATCACGCGACGGc                                             >  1:74780/1‑93 (MQ=255)

                                            cttTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGc                                             >  1:2200409/1‑93 (MQ=255)

                                            cttTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGc                                             <  2:2625411/93‑1 (MQ=255)

                                               tACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCttt                                          >  2:3088926/1‑93 (MQ=255)

                                                 cGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGc                                        <  2:930169/93‑1 (MQ=255)

                                                  gCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGcc                                       >  2:541753/1‑93 (MQ=255)

                                                   cAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCg                                      >  2:2443983/1‑93 (MQ=255)

                                                    aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGt                                     <  2:1544036/93‑1 (MQ=255)

                                                    aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGt                                     <  1:2005223/93‑1 (MQ=255)

                                                     gCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGtt                                    >  2:1150093/1‑93 (MQ=255)

                                                       cATTCGTGTTCTTTGCTGTATTTCAGTGCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTAc                                  <  1:1550031/93‑1 (MQ=255)

                                                       cATTCGTGTTCTTTGCTGTATTTCAATTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTAc                                  >  1:327512/1‑93 (MQ=255)

                                                        aTTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACg                                 >  1:2690724/1‑93 (MQ=255)

                                                        aTTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACg                                 >  1:1358879/1‑93 (MQ=255)

                                                        aTTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACg                                 >  2:2277048/1‑93 (MQ=255)

                                                         ttCGTGTTCTTTGTTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGc                                <  1:3709603/93‑1 (MQ=255)

                                                         ttCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGc                                >  2:3123184/1‑93 (MQ=255)

                                                         ttCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGc                                <  1:3772834/93‑1 (MQ=255)

                                                          tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGca                               <  2:3580460/93‑1 (MQ=255)

                                                           cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGcac                              <  1:2858375/93‑1 (MQ=255)

                                                            gtgtTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGcaca                             >  1:712774/1‑93 (MQ=255)

                                                             tgtTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaa                            >  1:1589879/1‑93 (MQ=255)

                                                              gtTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaa                           >  1:4588815/1‑93 (MQ=255)

                                                              gtTCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaa                           >  1:1108016/1‑93 (MQ=255)

                                                               ttCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaaa                          >  1:3387875/1‑93 (MQ=255)

                                                                tCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaaaa                         <  1:1898919/93‑1 (MQ=255)

                                                                tCTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaaaa                         >  2:886618/1‑93 (MQ=255)

                                                                 cTTTGCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAAc                        <  1:318567/93‑1 (MQ=255)

                                                                     gCTGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGCCGGCTTTGCCGTTACGCACAAAAACAGGt                    <  2:2906631/93‑1 (MQ=255)

                                                                       tGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt                  >  1:1852203/1‑93 (MQ=255)

                                                                       tGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt                  <  1:3247485/93‑1 (MQ=255)

                                                                       tGTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt                  <  2:2614316/93‑1 (MQ=255)

                                                                        gTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGTCTTTGCCGTTACGCACAAAAACAGGtttt                 <  2:4669974/93‑1 (MQ=255)

                                                                        gTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGtttt                 >  2:4565762/1‑93 (MQ=255)

                                                                        gTATTTCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGtttt                 <  1:2525825/93‑1 (MQ=255)

                                                                         tATTTCAGTTCTGCTGGCACGTTGCTCATAAATCAATCTCCTAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACACAAACAGGTTTTg                <  2:4209398/93‑1 (MQ=255)

                                                                           tttCAGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTc              <  2:3161346/93‑1 (MQ=255)

                                                                            ttCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCa             >  2:3890416/1‑93 (MQ=255)

                                                                               aGTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACtt          >  2:3825670/1‑93 (MQ=255)

                                                                                gTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACttt         >  2:2176627/1‑93 (MQ=255)

                                                                                gTTCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACttt         >  1:1685171/1‑93 (MQ=255)

                                                                                 ttCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTa        <  2:2397147/93‑1 (MQ=255)

                                                                                 ttCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTa        <  2:2410985/93‑1 (MQ=255)

                                                                                  tCTGCTGGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTaa       <  1:1265859/93‑1 (MQ=255)

                                                                                   ctgctgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAAc      <  2:1194945/93‑1 (MQ=255)

                                                                                   ctgctgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAAc      >  1:4355493/1‑93 (MQ=255)

                                                                                    tgctgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGATTTGTCACTTTAAcc     >  2:1827188/1‑93 (MQ=255)

                                                                                     gctgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAACCg    <  1:2988983/93‑1 (MQ=255)

                                                                                       tgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAACCGGc  >  1:493743/1‑93 (MQ=255)

                                                                                       tgGCACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAACCGGc  <  1:4560915/93‑1 (MQ=255)

CTGAGCATGTTCGGTAATACCAACGGTGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTTTGTCACTTTAACCGGC  >  NC_000913/3047468‑3047647

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