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.6%	809.9 / NA	232	intergenic (+126/‑38)	narQ/acrD	NarQ/AcrD
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (108/123); minor base C (1/0); total (109/123)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCC  >  NC_000913/2585487‑2585666

gTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAAcc                                                                                         >  2:5607098/1‑93 (MQ=255)

gTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAAcc                                                                                         <  2:1405552/93‑1 (MQ=255)

    aGTATAAAAGGGCATGATAATTTACTTTAACTCCTTTTTTTCTCCACGATTGGCTCGTAACTTGCAGCTTCCGTGAAGTCAGTCACACCTACa                                                                                     >  1:7459872/1‑93 (MQ=255)

     gTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACaa                                                                                    >  2:4784967/1‑93 (MQ=255)

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

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

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

         aaaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGAt                                                                                <  2:6150446/93‑1 (MQ=255)

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

          aaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATa                                                                               >  1:5299048/1‑93 (MQ=255)

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

              ggCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCa                                                                           >  1:867176/1‑93 (MQ=255)

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

                 aTGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaa                                                                        <  1:3782271/93‑1 (MQ=255)

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

                  tGATAATTTACATTAACTCCTTTGTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaaa                                                                       <  1:960488/93‑1 (MQ=255)

                   gATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAac                                                                      >  2:5652768/1‑93 (MQ=255)

                    aTAATTTACATTAACTACTTTTTTTCTCCACGACTGGCTCGTGCCTTGCAGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAaca                                                                     <  2:3962567/93‑1 (MQ=255)

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

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    >  2:2056946/1‑93 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    >  1:4253948/1‑93 (MQ=255)

                      aaTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACg                                                                   >  2:2976457/1‑93 (MQ=255)

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

                       aTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGa                                                                  >  2:1866730/1‑93 (MQ=255)

                       aTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGa                                                                  >  1:2084635/1‑93 (MQ=255)

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

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGGCGCTACAATGCAGCAAGCCAAACCTACAACGATACGCAGAAACACgcg                                                                 >  1:220533/1‑91 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                 >  2:820395/1‑93 (MQ=255)

                         ttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAgg                                                                >  1:6361025/1‑93 (MQ=255)

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

                           aCATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTc                                                              >  1:2222554/1‑93 (MQ=255)

                           aCATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTc                                                              >  1:6946549/1‑93 (MQ=255)

                            cATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTcc                                                             <  2:6664495/93‑1 (MQ=255)

                             aTTAACTTCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCct                                                            <  1:5695791/93‑1 (MQ=255)

                             aTTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGCAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCct                                                            >  1:1176210/1‑93 (MQ=255)

                              ttAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCctc                                                           >  2:6721812/1‑93 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                         <  1:4549853/93‑1 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                         <  1:4193410/93‑1 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                         <  1:1655045/93‑1 (MQ=255)

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

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

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

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

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

                                    ccTTTTTTGCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTaa                                                     <  1:6909017/93‑1 (MQ=255)

                                     cTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAAt                                                    >  1:5717646/1‑93 (MQ=255)

                                      tttttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATg                                                   <  2:4619880/93‑1 (MQ=255)

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

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

                                       ttttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATgg                                                  <  2:1618489/93‑1 (MQ=255)

                                        tttttCTCCACGATTGGCTCGTACCTTGCTGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGc                                                 >  2:6712923/1‑93 (MQ=255)

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

                                           ttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCCAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGaa                                              >  1:1069029/1‑93 (MQ=255)

                                            tctcCACGATTGGCTCTTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAt                                             <  2:1115447/93‑1 (MQ=255)

                                             ctcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAtt                                            <  2:1859874/93‑1 (MQ=255)

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

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

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

                                                 aCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                        <  2:5948304/93‑1 (MQ=255)

                                                   gATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtttt                                      >  2:6422563/1‑92 (MQ=255)

                                                   gATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTa                                      >  1:4002038/1‑93 (MQ=255)

                                                   gATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTa                                      >  1:4956886/1‑93 (MQ=255)

                                                   gATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTa                                      <  2:1735190/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                     <  2:2222908/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                     <  1:1063006/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                     <  1:6330581/93‑1 (MQ=255)

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

                                                     ttGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAtt                                    >  2:5721117/1‑93 (MQ=255)

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

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

                                                      tGGCTCGTACCTTGCCGCTACAGTGAAGCAAGGCAAACCTACAACGATACGCAGAAACACGAGGGCCTCTTTTAATGGCGAATTTCTTTATTg                                   >  2:2301953/1‑93 (MQ=255)

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

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

                                                         cTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATc                                <  1:859481/93‑1 (MQ=255)

                                                             tACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGccc                            >  1:3561853/1‑93 (MQ=255)

                                                             tACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGccc                            >  1:5576530/1‑93 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          <  1:1745532/93‑1 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          >  2:7503825/1‑93 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          >  2:3038932/1‑93 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          >  2:426787/1‑93 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          >  2:4371660/1‑93 (MQ=255)

                                                                  tGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAttt                       <  1:369117/93‑1 (MQ=255)

                                                                    ccGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAttttt                     <  1:6480902/93‑1 (MQ=255)

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

                                                                       cTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGcc                  <  2:5154843/93‑1 (MQ=255)

                                                                         aCAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTg                >  2:1195616/1‑93 (MQ=255)

                                                                          cAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTgg               <  2:4160423/93‑1 (MQ=255)

                                                                           aGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTggg              >  1:5979312/1‑93 (MQ=255)

                                                                            gTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGt             <  1:7163729/93‑1 (MQ=255)

                                                                              gAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGc           <  2:6324826/93‑1 (MQ=255)

                                                                               aaGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATAGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCt          >  1:6008089/1‑93 (MQ=255)

                                                                                 gCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGAGTGCTgg        >  1:1449261/1‑93 (MQ=255)

                                                                                   aaGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCa      >  2:7652834/1‑93 (MQ=255)

                                                                                   aaGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCa      >  2:3788289/1‑93 (MQ=255)

                                                                                    aGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCaa     <  1:2428596/93‑1 (MQ=255)

                                                                                    aGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCaa     <  2:1965932/93‑1 (MQ=255)

                                                                                     gTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAAt    <  2:73058/93‑1 (MQ=255)

                                                                                      tCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATc   <  1:3227183/93‑1 (MQ=255)

                                                                                      tCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATc   <  1:4465778/93‑1 (MQ=255)

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

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCC  >  NC_000913/2585487‑2585666

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