BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,496	+C	coding (301/408 nt)	phnQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,495	1	.	C	98.7%	1816.0 / 9.9	460	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (268/186); total (270/190)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 521 total aligned reads displayed.

GATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGA  >  NC_000913/4311400‑4311585

gATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCa                                                                                             <  2:7989496/96‑1 (MQ=255)

 gtGAGTAATGATTCATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑c                                                                                           >  1:11383606/2‑96 (MQ=255)

 aTGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑c                                                                                           <  1:10373846/96‑1 (MQ=255)

  tGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCCCGCTGT‑CCA‑cc                                                                                          <  1:7635755/96‑2 (MQ=255)

  tGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑cc                                                                                          >  2:12135519/1‑95 (MQ=255)

  tGAGTAATGATTAATGTGGATAGAGCTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑cc                                                                                          >  1:9165319/1‑95 (MQ=255)

   gAGTAATGATTAATGTGGATAGAGTTTCTTTTTTAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccg                                                                                         >  2:5377253/1‑94 (MQ=255)

   gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCACGCTGT‑CCA‑ccg                                                                                         <  1:2515951/96‑3 (MQ=255)

     gTAATTATTAATGTGGAGAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccgtt                                                                                       <  1:3231019/96‑5 (MQ=255)

     gTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccgtt                                                                                       <  1:8919941/96‑5 (MQ=255)

        aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCg                                                                                     <  2:5604150/96‑1 (MQ=255)

         tGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCgg                                                                                    <  2:7750002/96‑1 (MQ=255)

          gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCACGCTGT‑CCACCGTTCCGGc                                                                                   <  2:4110144/96‑1 (MQ=255)

             tAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCa                                                                                <  1:10916507/96‑1 (MQ=255)

              aaTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAg                                                                               >  1:8120597/1‑96 (MQ=255)

                tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGc                                                                             >  1:11247832/1‑96 (MQ=255)

                 gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGcc                                                                            <  2:8161954/96‑1 (MQ=255)

                 gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGcc                                                                            >  1:8266695/1‑96 (MQ=255)

                 gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGcc                                                                            <  1:655685/96‑1 (MQ=255)

                      tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGatat                                                                       >  1:3684268/1‑96 (MQ=255)

                      tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGatat                                                                       >  2:3461784/1‑96 (MQ=255)

                       agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCCCCGTTCCGGCCCAGGCCGAtatt                                                                      <  2:10519599/96‑1 (MQ=255)

                       agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtatt                                                                      >  2:1900758/1‑96 (MQ=255)

                          gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATc                                                                   <  1:2217753/96‑1 (MQ=255)

                          gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATc                                                                   >  2:13420141/1‑96 (MQ=255)

                           tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcg                                                                  >  1:4477662/1‑96 (MQ=255)

                            ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcga                                                                 >  2:8890832/1‑96 (MQ=255)

                            ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcga                                                                 >  2:13306674/1‑96 (MQ=255)

                            ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcga                                                                 >  2:580049/1‑96 (MQ=255)

                              cTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGG‑CCACCGTTCCGGCCCAGGCCGATATTATcggcg                                                               <  1:13926622/96‑1 (MQ=255)

                               tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgacga                                                              <  2:3110471/96‑1 (MQ=255)

                                ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCgacgac                                                             >  1:5176050/1‑96 (MQ=255)

                                ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCgacgac                                                             >  2:11651248/1‑96 (MQ=255)

                                 tttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGaca                                                            <  2:366185/96‑1 (MQ=255)

                                  ttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgacgcaac                                                           >  1:6807304/1‑92 (MQ=255)

                                   tGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACAcc                                                          <  2:1418354/96‑1 (MQ=255)

                                    gAGGTTGGCTAATGAACCGCTTTCCCCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCa                                                         <  2:1648585/96‑1 (MQ=255)

                                       gTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                      <  2:8414338/96‑1 (MQ=255)

                                        ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAGCCTGCTCACGCTGT‑CCACCGTGCCGCCCCAGGCCGGTATTATCGACGACACCATTgg                                                     >  1:2268620/1‑96 (MQ=255)

                                        ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                     <  2:9484428/96‑1 (MQ=255)

                                        ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                     >  2:6575222/1‑96 (MQ=255)

                                         tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                                    <  1:2685631/96‑1 (MQ=255)

                                          ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                   >  1:2881685/1‑96 (MQ=255)

                                          ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                   >  1:14123283/1‑96 (MQ=255)

                                          cgCTAATGAAACTCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                   <  1:10759180/95‑1 (MQ=255)

                                            cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                                 >  1:13669129/1‑96 (MQ=255)

                                            cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                                 >  2:10802290/1‑96 (MQ=255)

                                              aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                               <  1:14196859/96‑1 (MQ=255)

                                              aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                               >  1:9481923/1‑96 (MQ=255)

                                                taaaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTGCGGCCCAGGCCGATATTATGGACGACACCATCGGCAATATTc                                             >  1:11114159/3‑96 (MQ=255)

                                                tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                             >  1:7240103/1‑96 (MQ=255)

                                                tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                             <  1:1855994/96‑1 (MQ=255)

                                                 gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTca                                            <  1:8514478/96‑1 (MQ=255)

                                                  aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                           >  2:9601786/1‑96 (MQ=255)

                                                   aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                          <  1:6410087/96‑1 (MQ=255)

                                                   aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                          <  2:11494970/96‑1 (MQ=255)

                                                      gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                       >  2:11102354/1‑96 (MQ=255)

                                                       cTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                                      <  2:5119254/96‑1 (MQ=255)

                                                        tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                                     >  1:479960/1‑96 (MQ=255)

                                                           ccGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                                  <  2:5497167/96‑1 (MQ=255)

                                                             gCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                                >  1:3460117/1‑96 (MQ=255)

                                                               tATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                              >  2:2275513/1‑96 (MQ=255)

                                                                aTTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGGTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                             >  1:1703277/1‑96 (MQ=255)

                                                                 tttCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                            >  2:2675243/1‑96 (MQ=255)

                                                                    cTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                         <  1:7938778/96‑1 (MQ=255)

                                                                      tttATTCCCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                       <  1:3033524/96‑1 (MQ=255)

                                                                       ttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAGTATTCAGCAGGCGATTAATGACGCTtat                      >  1:13261087/1‑96 (MQ=255)

                                                                       ttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                      >  1:9651090/1‑96 (MQ=255)

                                                                       ttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                      >  2:14348684/1‑96 (MQ=255)

                                                                        tATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCCACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtata                     >  2:6069647/1‑96 (MQ=255)

                                                                         aTTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                    >  1:9786063/1‑96 (MQ=255)

                                                                         aTTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                    >  1:4129927/1‑96 (MQ=255)

                                                                          ttCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                   >  2:11272731/1‑96 (MQ=255)

                                                                          ttCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                   >  2:7028767/1‑96 (MQ=255)

                                                                          ttCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                   >  2:126967/1‑96 (MQ=255)

                                                                           tCACCCTGCTCACGTTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc                  <  1:6219408/96‑1 (MQ=255)

                                                                           tCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc                  >  2:646968/1‑96 (MQ=255)

                                                                            cACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc                 <  2:4477662/96‑1 (MQ=255)

                                                                            cACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc                 >  2:13604623/1‑96 (MQ=255)

                                                                             aCCTTGCTCACGCTGTCCCC‑CGTCCAGGTCCAGGCTGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGTCGTTTATACCCCt                <  1:11564525/96‑1 (MQ=255)

                                                                             aCCCTGCTCCCGCTGTCCCC‑CGTTCCGGCCCAGGCCGATATTATCGACGACCCCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt                <  1:5491549/96‑1 (MQ=255)

                                                                             aCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt                >  2:10289732/1‑96 (MQ=255)

                                                                              cccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg               <  2:10995326/96‑1 (MQ=255)

                                                                              cccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg               <  1:12036671/96‑1 (MQ=255)

                                                                              cccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg               >  2:1237988/1‑96 (MQ=255)

                                                                                cTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc             >  2:5715541/1‑96 (MQ=255)

                                                                                cTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc             <  2:7345948/96‑1 (MQ=255)

                                                                                  gCTCCCGCTGT‑CCACCGTTCCGCCCCAGGCCGATATTATCGACGACACCATGGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGActg           <  1:13049717/96‑3 (MQ=255)

                                                                                    tCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg         <  2:9165319/96‑1 (MQ=255)

                                                                                    tCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg         <  1:10032774/96‑1 (MQ=255)

                                                                                     cACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCgggg        >  1:13437136/1‑96 (MQ=255)

                                                                                     cACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAACCTGACCGTgg        >  2:4848064/1‑96 (MQ=255)

                                                                                      aCGCTGTCCCC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCACTATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg       <  2:6790900/96‑1 (MQ=255)

                                                                                      aCGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg       >  2:3925738/1‑96 (MQ=255)

                                                                                      aCGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg       >  1:5018501/1‑96 (MQ=255)

                                                                                      aCGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgag       >  1:12878301/1‑94 (MQ=255)

                                                                                        gCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCg     >  2:14026347/1‑96 (MQ=255)

                                                                                        gCTGT‑CCACCGTTCCGGCCCAGGACGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACACTGACCGTGGGCg     >  1:8689246/1‑96 (MQ=255)

                                                                                          tGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCgag   >  1:5959211/1‑96 (MQ=255)

                                                                                           gT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCgaga  >  1:12373852/1‑96 (MQ=255)

GATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGA  >  NC_000913/4311400‑4311585

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