BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,995,357	+G	coding (562/717 nt)	yigB →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,353	1	.	G	99.1%	774.3 / ‑1.1	225	coding (558/717 nt)	yigB	putative phosphatase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base G (108/115); total (108/117)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 265 total aligned reads displayed.

GCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑GGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATG  >  NC_000913/3995261‑3995440

gcTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGtt                                                                                          >  1:4858907/1‑93 (MQ=255)

   ggCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑ggg                                                                                      >  2:4410888/1‑93 (MQ=255)

    gCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggg                                                                                     <  2:92551/93‑1 (MQ=255)

    gCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggg                                                                                     <  1:1598014/93‑1 (MQ=255)

    gCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggg                                                                                     <  1:3600182/93‑1 (MQ=255)

     cgcGCACGGGCGGTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑ggggg                                                                                    <  2:1433016/91‑2 (MQ=255)

     cccGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑ggggg                                                                                    <  1:3773301/93‑2 (MQ=255)

       cGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggggac                                                                                  >  2:687824/1‑90 (MQ=255)

       cGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggggac                                                                                  <  2:4965427/93‑4 (MQ=255)

       cGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggggac                                                                                  <  2:3281872/93‑4 (MQ=255)

        gCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑gggggccg                                                                                 >  2:2819945/1‑89 (MQ=255)

          aCGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGAt                                                                                <  2:609553/93‑1 (MQ=255)

          aCGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGAt                                                                                <  1:3100972/93‑1 (MQ=255)

              gcgcTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTca                                                                            >  1:4491735/1‑93 (MQ=255)

                gcTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTcacc                                                                          >  1:3708492/1‑93 (MQ=255)

                 cTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTcacca                                                                         >  1:4328034/1‑93 (MQ=255)

                  tCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTcaccac                                                                        >  2:3774933/1‑93 (MQ=255)

                    aaaaCCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTg                                                                      >  2:2112442/1‑93 (MQ=255)

                     aaaCCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGa                                                                     >  2:1580797/1‑93 (MQ=255)

                     aaaCCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGa                                                                     <  1:857099/93‑1 (MQ=255)

                     aaaCCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGa                                                                     <  1:4636347/93‑1 (MQ=255)

                       aCCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACg                                                                   >  1:5409042/1‑93 (MQ=255)

                        ccGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgt                                                                  <  1:5415123/93‑1 (MQ=255)

                          gtacAGCGATATGTACTTTTTGCCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgg                                                                <  1:5309762/90‑1 (MQ=255)

                          gTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgg                                                                >  2:1686892/1‑93 (MQ=255)

                          gTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgg                                                                >  2:1805916/1‑93 (MQ=255)

                             cAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgggtg                                                             >  1:6217710/1‑93 (MQ=255)

                             cAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgggtg                                                             >  2:2308075/1‑93 (MQ=255)

                              aGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACgtgggtgg                                                            <  1:1055616/93‑1 (MQ=255)

                               gCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGtgggtggg                                                           >  2:1358899/1‑93 (MQ=255)

                               gCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGtggggggg                                                           >  1:1926855/1‑93 (MQ=255)

                                 gATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGc                                                         >  1:2433020/1‑93 (MQ=255)

                                  atatGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCa                                                        >  1:3637405/1‑93 (MQ=255)

                                    atGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAAt                                                      >  1:3443788/1‑93 (MQ=255)

                                     tGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAAtt                                                     <  1:3352133/93‑1 (MQ=255)

                                     tGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAAtt                                                     >  2:46132/1‑93 (MQ=255)

                                      gTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTc                                                    <  1:6438126/93‑1 (MQ=255)

                                       tACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCg                                                   >  2:5983664/1‑93 (MQ=255)

                                         cTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCa                                                 <  1:2224857/93‑1 (MQ=255)

                                          tttttGGCTGCGGAAAAACTCAGCGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAg                                                <  1:5534876/93‑1 (MQ=255)

                                          tttttGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAg                                                >  1:1765599/1‑93 (MQ=255)

                                           ttttGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGc                                               >  2:3895492/1‑93 (MQ=255)

                                           ttttGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGc                                               >  2:3117999/1‑93 (MQ=255)

                                           ttttGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGc                                               <  1:494948/93‑1 (MQ=255)

                                               ggCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGaa                                           >  1:2142575/1‑93 (MQ=255)

                                                gCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAAt                                          <  1:6131027/93‑1 (MQ=255)

                                                gCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAAt                                          <  1:1677896/93‑1 (MQ=255)

                                                  tGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGc                                        <  2:4018050/93‑1 (MQ=255)

                                                  tGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGc                                        <  1:536937/93‑1 (MQ=255)

                                                  tGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGc                                        >  1:4745462/1‑93 (MQ=255)

                                                  cgCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGc                                        <  2:4702753/92‑1 (MQ=255)

                                                   gCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCa                                       <  2:27272/93‑1 (MQ=255)

                                                   gCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCa                                       >  1:4691571/1‑93 (MQ=255)

                                                    cGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAg                                      >  2:3557455/1‑93 (MQ=255)

                                                    cGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAg                                      <  1:1097529/93‑1 (MQ=255)

                                                    cGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAg                                      >  1:4456499/1‑93 (MQ=255)

                                                      gAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGc                                    <  2:2461247/93‑1 (MQ=255)

                                                             tCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTgg                             >  2:5119184/1‑93 (MQ=255)

                                                              cAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGa                            <  1:3025097/93‑1 (MQ=255)

                                                              cAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGa                            >  2:107939/1‑93 (MQ=255)

                                                              cAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGa                            <  1:5028526/93‑1 (MQ=255)

                                                              cAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGa                            <  2:4166373/93‑1 (MQ=255)

                                                              cAACGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGa                            <  2:3685384/93‑1 (MQ=255)

                                                               aaCGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGCAATGCAGGCTTGTTGGAt                           <  1:568048/93‑1 (MQ=255)

                                                                aCGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATc                          >  1:1380242/1‑93 (MQ=255)

                                                                aCGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATc                          >  1:939762/1‑93 (MQ=255)

                                                                 cGTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCa                         >  2:5665572/1‑93 (MQ=255)

                                                                  gTGCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGGTCAg                        >  2:3423037/1‑93 (MQ=255)

                                                                   tGCGGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGa                       <  1:4585444/93‑1 (MQ=255)

                                                                    gCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGAc                      <  2:5293352/93‑1 (MQ=255)

                                                                    gCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGAc                      >  2:5584843/1‑93 (MQ=255)

                                                                    gCCGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGAc                      >  2:6450053/1‑93 (MQ=255)

                                                                      cGATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCg                    <  1:5961538/93‑1 (MQ=255)

                                                                       gATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCgg                   >  1:6367824/1‑93 (MQ=255)

                                                                       gATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCgg                   >  2:3383628/1‑93 (MQ=255)

                                                                       gATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCgg                   <  1:4282113/93‑1 (MQ=255)

                                                                       gATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCgg                   <  1:5192049/93‑1 (MQ=255)

                                                                       gATCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCgg                   >  2:2525160/1‑93 (MQ=255)

                                                                        aTCGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGa                  >  2:1920952/1‑93 (MQ=255)

                                                                         tCGGCGAGATCTTACACGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGaa                 >  2:2925161/1‑93 (MQ=255)

                                                                          cGGCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGaaa                <  1:128252/93‑1 (MQ=255)

                                                                           ggCGAGATCTTACATGTTGGGGGACGATCTCCCCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGaaaa               <  2:846447/93‑1 (MQ=255)

                                                                           ggCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGaaaa               >  1:4056352/1‑93 (MQ=255)

                                                                           ggCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGAACAGACCGGgaaa               >  1:4683188/1‑93 (MQ=255)

                                                                            gCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGaaaag              >  2:4569762/1‑92 (MQ=255)

                                                                            gCGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAAt              <  1:2289481/93‑1 (MQ=255)

                                                                             cGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATg             <  1:4128157/93‑1 (MQ=255)

                                                                             cGAGATCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATg             <  2:4092926/93‑1 (MQ=255)

                                                                              gagaTCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATgg            <  1:6276143/93‑1 (MQ=255)

                                                                               agaTCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGc           >  2:369575/1‑93 (MQ=255)

                                                                               agaTCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGc           >  2:5811591/1‑93 (MQ=255)

                                                                               agaTCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGc           >  1:3988580/1‑93 (MQ=255)

                                                                                 aTCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGa         <  1:1004784/93‑1 (MQ=255)

                                                                                  tCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGAt        <  2:1245832/93‑1 (MQ=255)

                                                                                  tCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGAt        >  2:3212145/1‑93 (MQ=255)

                                                                                  tCTTACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGAt        <  2:3029274/93‑1 (MQ=255)

                                                                                    ttACATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCt      <  1:74873/93‑1 (MQ=255)

                                                                                      aCATGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGa    <  1:885223/93‑1 (MQ=255)

                                                                                        aTGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATg  >  1:1952693/1‑93 (MQ=255)

                                                                                        aTGTTGGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATg  <  2:2644856/93‑1 (MQ=255)

GCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTT‑GGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATG  >  NC_000913/3995261‑3995440

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