BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,724,581	T→C	Q137R (CAG→CGG)	yiaA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,724,581	0	T	C	99.8%	1425.6 / ‑5.5	452	Q137R (CAG→CGG)	yiaA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (235/216); total (235/217)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 471 total aligned reads displayed.

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

aGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCcgt                                                                                         <  1:453102/93‑1 (MQ=255)

 gAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCcgtg                                                                                        >  1:344730/1‑93 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       <  2:1804452/93‑1 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       <  2:1166144/93‑1 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       >  2:6353211/1‑93 (MQ=255)

   aaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGtt                                                                                      <  1:4924277/93‑1 (MQ=255)

     cAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTc                                                                                    <  1:206796/93‑1 (MQ=255)

     cAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTc                                                                                    >  2:1464058/1‑93 (MQ=255)

       gATTAGTGTGTTACGTAAAAGCCCGGTCACATTGGACGGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCCGTGTTTCtt                                                                                  <  1:3349823/93‑1 (MQ=255)

       gATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtt                                                                                  >  1:1623648/1‑93 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                >  1:2846889/1‑93 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                >  1:1485384/1‑93 (MQ=255)

           aGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtgg                                                                              <  2:3632883/93‑1 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             >  2:423261/1‑93 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             >  2:3995781/1‑93 (MQ=255)

             tGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtg                                                                            >  1:2591380/1‑93 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           <  2:6735232/93‑1 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           <  2:5460813/93‑1 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           <  1:6589578/93‑1 (MQ=255)

               aGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTggg                                                                          >  1:6313533/1‑93 (MQ=255)

                gtgaaGTAAAAGCCGGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         <  1:5264048/90‑1 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGTTGGGt                                                                         >  1:4761777/1‑93 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         >  1:7930312/1‑93 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         >  1:817944/1‑93 (MQ=255)

                 ttAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGtt                                                                        >  1:5467178/1‑93 (MQ=255)

                   aaGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTa                                                                      >  2:2679500/1‑93 (MQ=255)

                      tAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTc                                                                   <  1:5404902/93‑1 (MQ=255)

                       aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTcc                                                                  >  1:1747891/1‑93 (MQ=255)

                        aaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCg                                                                 >  2:4859930/1‑93 (MQ=255)

                        aaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCg                                                                 >  2:5337619/1‑93 (MQ=255)

                        aaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCg                                                                 >  2:7332955/1‑93 (MQ=255)

                         aaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCgg                                                                <  1:1994534/93‑1 (MQ=255)

                         aaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCgg                                                                >  2:2370615/1‑93 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              >  1:5570498/1‑93 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             >  2:2227313/1‑93 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             <  2:888405/93‑1 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             <  1:4145981/93‑1 (MQ=255)

                             ccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATc                                                            >  1:4410375/1‑93 (MQ=255)

                              cGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCa                                                           >  1:510599/1‑93 (MQ=255)

                               ggTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCAc                                                          >  1:7379701/1‑93 (MQ=255)

                                gTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                                         >  1:2306004/1‑93 (MQ=255)

                                  cacaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGaa                                                       <  2:7128586/93‑1 (MQ=255)

                                  cacaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGaa                                                       <  2:7631100/93‑1 (MQ=255)

                                  cacaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGaa                                                       >  1:733651/1‑93 (MQ=255)

                                   tcaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAat                                                      <  2:4162304/92‑1 (MQ=255)

                                   acaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAat                                                      >  1:4279504/1‑93 (MQ=255)

                                   acaTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAat                                                      <  2:1274525/93‑1 (MQ=255)

                                    caTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAata                                                     <  1:5825522/93‑1 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  1:1139684/1‑93 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    <  1:6892312/93‑1 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  1:6010686/1‑93 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  1:241087/1‑93 (MQ=255)

                                      ttGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtt                                                   >  2:3724538/1‑93 (MQ=255)

                                      ttGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtt                                                   >  2:7516176/1‑93 (MQ=255)

                                        ggACTGACCGGGCTTACGTGATTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtttt                                                 >  1:1380611/1‑93 (MQ=255)

                                        ggACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtttt                                                 <  2:4449005/93‑1 (MQ=255)

                                             gACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGc                                            >  2:4389227/1‑93 (MQ=255)

                                              aCCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCa                                           <  1:1825584/93‑1 (MQ=255)

                                                  ggCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGc                                       >  1:764523/1‑93 (MQ=255)

                                                  ggCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGc                                       >  2:1998166/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    >  2:6878417/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    >  1:2749052/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    >  2:3296999/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    <  2:4177950/93‑1 (MQ=255)

                                                      tACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACt                                   >  2:7805612/1‑93 (MQ=255)

                                                        cGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGc                                 >  2:1324672/1‑93 (MQ=255)

                                                        cGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGc                                 >  1:5149202/1‑93 (MQ=255)

                                                         gTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGGTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCa                                >  1:3814414/1‑93 (MQ=255)

                                                          tGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCAc                               >  2:6143485/1‑93 (MQ=255)

                                                             gTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaa                            >  1:463558/1‑93 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           >  2:2167679/1‑93 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           >  1:7487109/1‑93 (MQ=255)

                                                               tATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAAc                          >  2:6003022/1‑93 (MQ=255)

                                                                aTTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACa                         >  1:7836766/1‑93 (MQ=255)

                                                                 ttCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACaa                        >  1:1328589/1‑93 (MQ=255)

                                                                 ttCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACaa                        >  2:4651666/1‑93 (MQ=255)

                                                                    gCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCt                     <  2:1739325/93‑1 (MQ=255)

                                                                    gCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCAGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCt                     >  2:5584453/1‑93 (MQ=255)

                                                                    gCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTAGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCt                     <  1:6364067/93‑1 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   >  1:7142551/1‑93 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   >  1:5655307/1‑93 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   >  2:5873232/1‑93 (MQ=255)

                                                                        tATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaag                 <  2:2864902/93‑1 (MQ=255)

                                                                         aTTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaaga                <  2:4542043/93‑1 (MQ=255)

                                                                          ttcttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaagaa               <  1:3207916/93‑1 (MQ=255)

                                                                            cttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAagaaga             <  1:7901450/93‑1 (MQ=255)

                                                                             ttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAgaagaa            >  1:1041489/1‑93 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  1:2704265/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           >  1:4037155/1‑93 (MQ=255)

                                                                               ccTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAg          <  2:5158237/93‑1 (MQ=255)

                                                                               ccTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAg          <  1:6176650/93‑1 (MQ=255)

                                                                                cTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc         <  2:7685455/93‑1 (MQ=255)

                                                                                cTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc         >  1:3945281/1‑93 (MQ=255)

                                                                                 tGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGcc        >  2:5701154/1‑93 (MQ=255)

                                                                                 tGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGcc        >  1:4349081/1‑93 (MQ=255)

                                                                                   ggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAg      >  2:7161582/1‑93 (MQ=255)

                                                                                     tAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTc    <  1:3641515/93‑1 (MQ=255)

                                                                                     tAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTc    <  2:3242310/93‑1 (MQ=255)

                                                                                      aaCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTcc   <  2:1670805/93‑1 (MQ=255)

                                                                                       aCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCa  <  2:2994433/93‑1 (MQ=255)

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

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