BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,585,577	G→A	intergenic (+126/‑38)	narQ → / → acrD	NarQ/AcrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,585,577	0	G	A	99.4%	553.3 / ‑5.2	161	intergenic (+126/‑38)	narQ/acrD	NarQ/AcrD
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (73/87); total (73/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 169 total aligned reads displayed.

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCCT  >  NC_000913/2585487‑2585667

gTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCACCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAAcc                                                                                          <  1:3802043/93‑1 (MQ=255)

 ttCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCt                                                                                         >  1:67761/1‑93 (MQ=255)

 ttCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTACACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCt                                                                                         <  1:3326419/93‑1 (MQ=255)

  tCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTa                                                                                        <  1:1019300/93‑1 (MQ=255)

  tCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTa                                                                                        >  1:4146213/1‑93 (MQ=255)

   cAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTAc                                                                                       <  2:1647511/93‑1 (MQ=255)

    aGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACa                                                                                      >  2:1299034/1‑93 (MQ=255)

     gTATAAAAGGGCATGATAATTTACATTAGCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACaa                                                                                     <  1:1016952/93‑1 (MQ=255)

     gTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACaa                                                                                     <  2:104794/93‑1 (MQ=255)

      tataAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAAc                                                                                    <  2:722053/93‑1 (MQ=255)

      tataAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAAc                                                                                    >  1:959233/1‑93 (MQ=255)

        taAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGa                                                                                  <  2:4137666/93‑1 (MQ=255)

        taAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGa                                                                                  <  1:3961365/93‑1 (MQ=255)

          aaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATa                                                                                <  1:2425687/93‑1 (MQ=255)

          aaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCATGCCGCTACAGTGAAGCAAGACAAACCTAAAACGATa                                                                                >  2:3619666/1‑93 (MQ=255)

           aaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATAc                                                                               <  1:4431972/93‑1 (MQ=255)

           aaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATAc                                                                               <  2:4001622/93‑1 (MQ=255)

            aGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACg                                                                              <  1:4239097/93‑1 (MQ=255)

               gCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAg                                                                           >  2:186978/1‑93 (MQ=255)

                cATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGa                                                                          <  1:3900823/93‑1 (MQ=255)

                 aTGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaa                                                                         >  1:1714370/1‑93 (MQ=255)

                  tGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaaa                                                                        >  1:4355549/1‑93 (MQ=255)

                  tGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaaa                                                                        >  1:1419925/1‑93 (MQ=255)

                  tGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaaa                                                                        >  2:720335/1‑93 (MQ=255)

                   gATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAac                                                                       <  1:2913553/93‑1 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                     <  2:4559430/93‑1 (MQ=255)

                      aaTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACg                                                                    <  1:3148579/93‑1 (MQ=255)

                       aTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGa                                                                   <  2:1173550/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                  <  1:88283/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                  <  2:4167712/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                  <  1:3922543/93‑1 (MQ=255)

                          tACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGt                                                                >  2:643407/1‑93 (MQ=255)

                           aCATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTc                                                               <  2:247157/93‑1 (MQ=255)

                              ttAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCctc                                                            <  1:4495418/93‑1 (MQ=255)

                               tAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCctct                                                           <  2:4416029/93‑1 (MQ=255)

                               tAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCctct                                                           >  2:3570185/1‑93 (MQ=255)

                               tAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCctct                                                           >  1:2817413/1‑93 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                          >  1:2343144/1‑93 (MQ=255)

                                 aCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCttt                                                         <  1:215267/93‑1 (MQ=255)

                                 aCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCttt                                                         <  1:1538559/93‑1 (MQ=255)

                                  cTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtttt                                                        >  2:2717102/1‑93 (MQ=255)

                                  cTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGGCCTCtttt                                                        >  2:1266467/1‑93 (MQ=255)

                                   tCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCTCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTa                                                       <  1:1983798/93‑1 (MQ=255)

                                   tCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTa                                                       >  2:2299080/1‑93 (MQ=255)

                                   tCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAAACCGAGGTCCTCTTTTa                                                       >  1:94057/1‑93 (MQ=255)

                                    ccTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTaa                                                      <  2:1462020/93‑1 (MQ=255)

                                    ccTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTaa                                                      <  2:2336026/93‑1 (MQ=255)

                                      tttttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATg                                                    >  2:3273589/1‑93 (MQ=255)

                                       ttttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATgg                                                   >  2:354456/1‑93 (MQ=255)

                                         ttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCg                                                 >  2:3007337/1‑93 (MQ=255)

                                         ttttCTCCACGATTGGCTCGTAACTTGACGCTACAGTGAAGCAAGTCTAACCTACAACGATACGCAGAAACACGGGGACCTCTTTTCATGGCg                                                 >  2:2043646/1‑93 (MQ=255)

                                          tttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGa                                                <  2:4599971/93‑1 (MQ=255)

                                           ttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGaa                                               <  2:2100256/93‑1 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGCTACGCAGAAACACGAGGTCCTATTTTAATGGCGACttt                                            >  1:3003845/1‑93 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAttt                                            >  1:3284499/1‑93 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAttt                                            >  1:370773/1‑93 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAttt                                            >  1:1444078/1‑93 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAttt                                            >  1:1470533/1‑93 (MQ=255)

                                              tcCACGATTGGCTCGTACCCTGCCCCTAGATCGAAGCACGTCAAACCTAGAACGATACGCAGAAACACGAGCGCCTCTTTTAGTGGCGAAttt                                            >  2:2224511/1‑93 (MQ=255)

                                               ccACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTCATGGCGAATTTc                                           <  1:2530601/93‑1 (MQ=255)

                                               ccACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTc                                           <  2:599308/93‑1 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                          >  2:1959464/1‑93 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                          >  2:1459465/1‑93 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                          >  2:4518024/1‑93 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                          >  1:4566378/1‑93 (MQ=255)

                                                 ccGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                         <  1:4137033/92‑1 (MQ=255)

                                                 aCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                         >  1:1576020/1‑93 (MQ=255)

                                                 aCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCCTTTAATGGCGAATTTCtt                                         <  2:508592/93‑1 (MQ=255)

                                                  cGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCttt                                        >  1:3155624/1‑93 (MQ=255)

                                                  cGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCttt                                        >  1:613734/1‑93 (MQ=255)

                                                  cGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCttt                                        <  1:2122621/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                      >  2:3947066/1‑93 (MQ=255)

                                                      tGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTg                                    <  1:1903442/93‑1 (MQ=255)

                                                      tGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTg                                    <  1:4308936/93‑1 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGa                                   >  2:2534489/1‑93 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGa                                   <  2:3721940/93‑1 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGa                                   <  2:2097847/93‑1 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                  >  1:2382795/1‑93 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                  <  1:1839938/93‑1 (MQ=255)

                                                              aCCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGcccc                            >  1:405284/1‑93 (MQ=255)

                                                              aCCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGGCCTCTTTTAATGGCGAATTTCTTTATTGATCGcccc                            >  1:500338/1‑93 (MQ=255)

                                                                 ttGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAtt                         >  1:3828104/1‑93 (MQ=255)

                                                                   gCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAtttt                       <  1:1863987/93‑1 (MQ=255)

                                                                      gCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGc                    >  1:3543814/1‑93 (MQ=255)

                                                                      gCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGc                    <  2:2522394/93‑1 (MQ=255)

                                                                        tACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCt                  <  2:2851093/93‑1 (MQ=255)

                                                                        tACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCt                  <  2:2674913/93‑1 (MQ=255)

                                                                        tACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCt                  <  2:3503195/93‑1 (MQ=255)

                                                                         aCAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTg                 <  2:867846/93‑1 (MQ=255)

                                                                          cAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTgg                >  1:2682034/1‑93 (MQ=255)

                                                                           aGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTggg               <  2:4314494/93‑1 (MQ=255)

                                                                                aGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTg          <  2:1979481/93‑1 (MQ=255)

                                                                                aGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTg          >  2:3390510/1‑93 (MQ=255)

                                                                                aGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTg          <  2:2758907/93‑1 (MQ=255)

                                                                                    aGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCaa      >  2:3601925/1‑93 (MQ=255)

                                                                                    aGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCaa      >  1:3329161/1‑93 (MQ=255)

                                                                                       cAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATcc   <  1:3699349/93‑1 (MQ=255)

                                                                                       cAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATcc   <  2:973177/93‑1 (MQ=255)

                                                                                        aaGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCCt  <  2:4353988/93‑1 (MQ=255)

                                                                                        aaGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCCt  <  2:296092/93‑1 (MQ=255)

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCCT  >  NC_000913/2585487‑2585667

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