BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,965,001	T→C	M486V (ATG→GTG)	ptsP ←	PtsP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,965,001	0	T	C	99.6%	814.7 / ‑4.7	253	M486V (ATG→GTG)	ptsP	PtsP
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (119/133); total (120/133)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 268 total aligned reads displayed.

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

tttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         <  1:5596252/92‑1 (MQ=255)

ctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         >  1:4094400/1‑93 (MQ=255)

 ttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAcccc                                                                                        <  1:1728459/93‑1 (MQ=255)

 ttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAcccc                                                                                        <  1:2737462/93‑1 (MQ=255)

   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  2:3584387/93‑1 (MQ=255)

   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  1:2185717/93‑1 (MQ=255)

      tGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGAt                                                                                   >  1:3734203/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  <  1:2646122/93‑1 (MQ=255)

        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATat                                                                                 >  1:5959358/1‑92 (MQ=255)

        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATAc                                                                                 >  1:4582871/1‑93 (MQ=255)

        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATAc                                                                                 >  1:5558395/1‑93 (MQ=255)

          cGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTg                                                                               >  1:1511828/1‑93 (MQ=255)

          cGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTg                                                                               <  2:6447731/93‑1 (MQ=255)

           ggCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgc                                                                              >  2:2123973/1‑93 (MQ=255)

           ggCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgc                                                                              >  2:5175042/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  1:921495/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          <  1:6714509/93‑1 (MQ=255)

                gTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGCTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAc                                                                         >  1:6129582/1‑93 (MQ=255)

                 tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCCAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAcc                                                                        >  1:3413671/1‑93 (MQ=255)

                 tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCGCCCCCTGATACTGCGCcccc                                                                        <  1:1782272/93‑1 (MQ=255)

                  aaGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCt                                                                       >  1:2544731/1‑93 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  2:3489581/1‑93 (MQ=255)

                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     <  2:3235633/93‑1 (MQ=255)

                        gctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttctt                                                                 >  1:6273157/1‑93 (MQ=255)

                        gctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttctt                                                                 >  2:844754/1‑93 (MQ=255)

                         ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttc                                                                <  2:6594855/93‑1 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               <  2:3663304/93‑1 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               <  1:4340229/93‑1 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:5232290/1‑93 (MQ=255)

                            ctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttcttc                                                             <  1:2666504/93‑1 (MQ=255)

                            ctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttcttc                                                             >  1:3778037/1‑93 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  2:104838/93‑1 (MQ=255)

                              tATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCCCCCCCTGATACTGCGCCACCTGTTCTTCTTCCg                                                           <  1:1476955/93‑1 (MQ=255)

                              tATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCg                                                           >  1:5944493/1‑93 (MQ=255)

                               aTCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGa                                                          >  1:5820375/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         >  1:4411351/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         <  2:6518128/93‑1 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         >  2:5384602/1‑93 (MQ=255)

                                  tGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACgg                                                       <  2:273374/93‑1 (MQ=255)

                                   gCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACggg                                                      >  1:1308265/1‑93 (MQ=255)

                                    cTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGa                                                     <  1:3214682/93‑1 (MQ=255)

                                     tCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGaa                                                    >  1:4630024/1‑93 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCCCCTGTTCTTCTTCCGACGGGAAAcc                                                 <  1:351263/93‑1 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 >  2:3073868/1‑93 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 >  1:4005273/1‑93 (MQ=255)

                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCGGGCGGGAAACCa                                                >  2:1098334/1‑93 (MQ=255)

                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                >  2:3047840/1‑93 (MQ=255)

                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                >  2:5518353/1‑93 (MQ=255)

                                          cATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCAc                                               >  2:4057819/1‑93 (MQ=255)

                                          cATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCAc                                               <  1:5162418/93‑1 (MQ=255)

                                          cATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCAc                                               <  1:1137458/93‑1 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              <  1:170954/93‑1 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:839836/1‑93 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             <  1:2001993/93‑1 (MQ=255)

                                             ccAGCGTACGCAAGGTGACGTGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACttt                                            <  1:2050338/93‑1 (MQ=255)

                                               aGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgc                                          >  1:3893420/1‑93 (MQ=255)

                                               aGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgc                                          <  1:3965536/93‑1 (MQ=255)

                                                   tACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCGCCTGTTCTTCTTCCGACGGGAAACCACTTTgcagca                                      <  2:1568745/93‑1 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     <  2:383009/93‑1 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     <  2:2847482/93‑1 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  >  1:1902392/1‑93 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  2:5450013/93‑1 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  2:2189455/93‑1 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  2:1642438/93‑1 (MQ=255)

                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 <  2:2275106/93‑1 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                >  2:2439308/1‑93 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                >  1:2343596/1‑93 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                <  2:4882698/93‑1 (MQ=255)

                                                             gACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAt                            >  2:3354501/1‑93 (MQ=255)

                                                              aCGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                           <  1:5483333/93‑1 (MQ=255)

                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  2:1559452/93‑1 (MQ=255)

                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  2:172055/93‑1 (MQ=255)

                                                                 ggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCa                        <  2:811362/93‑1 (MQ=255)

                                                                 ggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCa                        >  1:2376894/1‑93 (MQ=255)

                                                                  gTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAg                       <  1:6424820/93‑1 (MQ=255)

                                                                   tttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      <  2:2287415/93‑1 (MQ=255)

                                                                    ttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTg                     <  1:4726996/93‑1 (MQ=255)

                                                                     tATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGc                    >  2:6177955/1‑93 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  <  1:1914162/93‑1 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  >  2:950136/1‑93 (MQ=255)

                                                                         aTTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATa                <  2:3380581/93‑1 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  1:2126926/1‑93 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  2:6084905/1‑93 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  2:2471359/1‑93 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               <  1:1092821/93‑1 (MQ=255)

                                                                           tGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaaa              <  2:6696405/93‑1 (MQ=255)

                                                                           tGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaaa              <  1:6201064/93‑1 (MQ=255)

                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             >  2:764617/1‑93 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            <  2:5604316/93‑1 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            >  2:4219564/1‑93 (MQ=255)

                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:4412224/93‑1 (MQ=255)

                                                                               cATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAcc          >  1:4281680/1‑93 (MQ=255)

                                                                                aTTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCg         <  2:6334567/93‑1 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        >  2:1424808/1‑93 (MQ=255)

                                                                                  ttGCAGCACCCCCTGATACTGCGCCGCCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAGAAAGACCGAt       <  1:3158231/93‑1 (MQ=255)

                                                                                    gcagcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGc     <  2:5294406/93‑1 (MQ=255)

                                                                                     cagcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGcc    >  2:2579378/1‑93 (MQ=255)

                                                                                      agcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   >  2:5025658/1‑93 (MQ=255)

                                                                                      agcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   <  2:1151015/93‑1 (MQ=255)

                                                                                       gcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAt  >  1:1630883/1‑93 (MQ=255)

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

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