BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,765,259	T→C	D484G (GAC→GGC)	ydiJ ←	putative oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,765,259	0	T	C	87.1%	563.3 / 63.1	209	D484G (GAC→GGC)	ydiJ	putative oxidase
Reads supporting (aligned to +/- strand): ref base T (11/16); new base C (96/86); total (107/102)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 226 total aligned reads displayed.

CACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACGTGAC  >  NC_000913/1765167‑1765349

   cGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGccat                                                                                         <  1:1009945/93‑1 (MQ=255)

    gAAAAACGCCTGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCAACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGccatc                                                                                        <  1:2817293/93‑1 (MQ=255)

    gAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACGACCGTATTTCGCAGTCAGCGCCACCACGccatc                                                                                        <  1:6647020/93‑1 (MQ=255)

    gAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGccatc                                                                                        <  2:2389026/93‑1 (MQ=255)

       aaaCGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCaga                                                                                     >  1:1724507/1‑93 (MQ=255)

        aaCGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCagag                                                                                    <  1:667204/93‑1 (MQ=255)

         aCGCCGGGCTGTATTCAGCGCGCAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGACATCAgaga                                                                                   >  1:1544365/1‑93 (MQ=255)

         aCGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAgaga                                                                                   <  2:357878/93‑1 (MQ=255)

          cGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAt                                                                                  <  1:5242935/93‑1 (MQ=255)

           gCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAtt                                                                                 >  2:2118278/1‑93 (MQ=255)

           gCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAtt                                                                                 >  2:1282482/1‑93 (MQ=255)

            ccGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAttt                                                                                >  2:3706272/1‑93 (MQ=255)

            ccGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAttt                                                                                <  1:2956670/93‑1 (MQ=255)

            ccGGGCTGTATTCAGCGCGAAAACCCTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGAttt                                                                                <  2:1620229/93‑1 (MQ=255)

              gggCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGc                                                                              >  1:4988290/1‑93 (MQ=255)

              gggCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGc                                                                              >  2:5456971/1‑93 (MQ=255)

              gggCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGc                                                                              >  1:3509737/1‑93 (MQ=255)

               ggCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACGAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCt                                                                             >  1:4135873/1‑93 (MQ=255)

                 cTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTtc                                                                           <  1:3339426/93‑1 (MQ=255)

                  tGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTtca                                                                          >  1:2171511/1‑93 (MQ=255)

                   gTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTtcat                                                                         <  1:2112781/93‑1 (MQ=255)

                   gTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTtcat                                                                         >  2:2558436/1‑93 (MQ=255)

                   ctATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTtcat                                                                         <  2:1222245/92‑1 (MQ=255)

                    tATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTtcatc                                                                        >  1:5736935/1‑93 (MQ=255)

                       tCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCaaa                                                                     >  2:5984200/1‑93 (MQ=255)

                       tCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTAGTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCaaa                                                                     <  2:1080545/93‑1 (MQ=255)

                         aGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAAt                                                                   >  2:3762072/1‑93 (MQ=255)

                          gcgcGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAAtc                                                                  >  1:6518407/1‑93 (MQ=255)

                          gcgcGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAAtc                                                                  <  1:4177038/93‑1 (MQ=255)

                          gcgcGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAAtc                                                                  >  2:5825082/1‑93 (MQ=255)

                             cgAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATctct                                                               <  1:5547012/93‑1 (MQ=255)

                              gAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCtt                                                              <  2:3324690/93‑1 (MQ=255)

                              gAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCtt                                                              >  2:1001880/1‑93 (MQ=255)

                               aaaaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttg                                                             >  2:6114759/1‑93 (MQ=255)

                                aaaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCttgt                                                            <  2:2528667/93‑1 (MQ=255)

                                aaaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgt                                                            <  1:4229086/93‑1 (MQ=255)

                                aaaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgt                                                            <  2:1864969/93‑1 (MQ=255)

                                aaaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgt                                                            >  1:4314603/1‑93 (MQ=255)

                                 aaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgtt                                                           >  1:2281794/1‑93 (MQ=255)

                                 aaCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgtt                                                           <  1:2867027/93‑1 (MQ=255)

                                  aCCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCttgttg                                                          >  2:1645377/1‑93 (MQ=255)

                                   ccTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGa                                                         >  2:345444/1‑93 (MQ=255)

                                    ctctGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAg                                                        <  1:6315382/90‑1 (MQ=255)

                                    cTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAg                                                        <  1:2044170/93‑1 (MQ=255)

                                    cTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAg                                                        <  1:5896463/93‑1 (MQ=255)

                                     tttGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAgg                                                       >  1:5339733/1‑93 (MQ=255)

                                      ttGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGa                                                      <  2:1696829/93‑1 (MQ=255)

                                       tGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGAt                                                     <  2:4609470/93‑1 (MQ=255)

                                         ccGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCg                                                   >  2:4382050/1‑93 (MQ=255)

                                           gTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGca                                                 <  2:3729779/93‑1 (MQ=255)

                                           gTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGca                                                 <  2:5952524/93‑1 (MQ=255)

                                            tGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGcac                                                >  2:724327/1‑93 (MQ=255)

                                             gCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGcaca                                               <  2:3149674/93‑1 (MQ=255)

                                                 gCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATc                                           <  1:3088473/93‑1 (MQ=255)

                                                 gCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATc                                           >  1:2346713/1‑93 (MQ=255)

                                                  ccccACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATcc                                          <  2:233591/93‑1 (MQ=255)

                                                    ccACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAg                                        <  2:6088426/93‑1 (MQ=255)

                                                    ccACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAg                                        >  1:2439424/1‑93 (MQ=255)

                                                    ccACAACAAACCACCGCATTTCGCAGGCAGCGACACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAAGAGCGCACATATcctg                                        >  1:3261614/1‑91 (MQ=255)

                                                     cacaACAAACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAgc                                       <  2:2861070/93‑1 (MQ=255)

                                                      ccaacaaACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAgcg                                      <  2:4482865/92‑1 (MQ=255)

                                                        aacaaACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCt                                    <  1:5420933/93‑1 (MQ=255)

                                                        aacaaACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCt                                    <  1:3654819/93‑1 (MQ=255)

                                                         acaaACCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTg                                   <  1:1865908/93‑1 (MQ=255)

                                                          caaACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTgg                                  <  1:1349462/93‑1 (MQ=255)

                                                          aaaaCCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTgg                                  <  2:6316130/92‑1 (MQ=255)

                                                            aaCCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggac                                <  1:399583/93‑1 (MQ=255)

                                                            aaCCACCGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATACAGCGCTggac                                <  2:3282159/93‑1 (MQ=255)

                                                             accaccGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacg                               >  2:4175992/1‑93 (MQ=255)

                                                             accaccGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacg                               <  1:2695308/93‑1 (MQ=255)

                                                             accaccGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacg                               <  1:3798601/93‑1 (MQ=255)

                                                              ccaccGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacgg                              >  2:4788993/1‑93 (MQ=255)

                                                               caccGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacgga                             <  1:1102855/93‑1 (MQ=255)

                                                               caccGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTggacgga                             <  1:560801/93‑1 (MQ=255)

                                                                 ccGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGgacggacg                           <  1:6158697/93‑1 (MQ=255)

                                                                  cGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGt                          <  2:84608/93‑1 (MQ=255)

                                                                  cGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGt                          >  2:1744023/1‑93 (MQ=255)

                                                                  cGTATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATACAGCGCTGGACGGACGt                          <  2:909746/93‑1 (MQ=255)

                                                                    tATTTCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGc                        >  1:3712405/1‑93 (MQ=255)

                                                                      gttCGCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCaa                      >  1:6382237/2‑93 (MQ=255)

                                                                         cGCAGTCACCCCCACCACGCCATCAGATGTTTGCCTCATCAAAATCTCGTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGAGCAAGac                   <  2:1718580/93‑1 (MQ=255)

                                                                          gCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGaca                  <  2:4937208/93‑1 (MQ=255)

                                                                          gCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGaca                  <  1:11573/93‑1 (MQ=255)

                                                                          gCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGaca                  >  2:6415113/1‑93 (MQ=255)

                                                                          gCAGTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGaca                  <  1:4066863/93‑1 (MQ=255)

                                                                           cAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGacac                 <  1:2500848/93‑1 (MQ=255)

                                                                             gTCAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCt               >  2:3092285/1‑93 (MQ=255)

                                                                               cAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGc             >  1:6113314/1‑93 (MQ=255)

                                                                               cAGCGCCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGc             >  2:3604047/1‑93 (MQ=255)

                                                                                 gcgcCACCACGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGt           >  1:1686408/1‑93 (MQ=255)

                                                                                  cgcCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTc          <  2:4955809/93‑1 (MQ=255)

                                                                                  cgcCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTc          <  2:3984010/93‑1 (MQ=255)

                                                                                  cgcCACCACGCCAACAGATATTTGCTTCATCAAAACCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCCAGACACCTGCGTc          <  1:4568511/93‑1 (MQ=255)

                                                                                    ccaccaCGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCgc        >  1:179283/1‑92 (MQ=255)

                                                                                    ccaccaCGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCAGCGTCGa        <  1:5988908/93‑1 (MQ=255)

                                                                                     caccacGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGAc       >  2:2215556/1‑93 (MQ=255)

                                                                                      accacGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACg      >  1:4383610/1‑93 (MQ=255)

                                                                                      accacGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACg      >  2:1972739/1‑93 (MQ=255)

                                                                                       ccacGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACGt     <  1:6208771/93‑1 (MQ=255)

                                                                                         acGCCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACGTGa   >  2:5118609/1‑93 (MQ=255)

CACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAAATCTCTTGTTGAGGATCGCACATATCCAGCGCTGGACGGACGTGCAAGACACCTGCGTCGACGTGAC  >  NC_000913/1765167‑1765349

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