BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,422,211	T→C	A53A (GCA→GCG)	yjfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,422,211	0	T	C	99.2%	849.9 / ‑5.7	265	A53A (GCA→GCG)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (131/132); total (132/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

cATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACgcgc                                                                                         <  1:4869253/93‑1 (MQ=255)

 atatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCt                                                                                        <  2:4660144/93‑1 (MQ=255)

 atatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCt                                                                                        >  1:198119/1‑93 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  2:6190737/1‑93 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  2:3108105/1‑93 (MQ=255)

    ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCg                                                                                     >  1:1721183/1‑93 (MQ=255)

     tGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGc                                                                                    <  2:732753/93‑1 (MQ=255)

      gCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCt                                                                                   <  1:2620008/93‑1 (MQ=255)

      gCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCt                                                                                   <  1:4358280/93‑1 (MQ=255)

      gCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCt                                                                                   <  1:3026399/93‑1 (MQ=255)

       ccGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGcgcttcgctt                                                                                  <  1:69154/93‑1 (MQ=255)

       ccGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGcgcttcgctt                                                                                  >  2:2806453/1‑93 (MQ=255)

         gTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTc                                                                                <  1:1032484/93‑1 (MQ=255)

          ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCtttta                                                                               <  2:5009883/93‑3 (MQ=255)

          ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               >  2:3761832/1‑93 (MQ=255)

          ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               <  1:516501/93‑1 (MQ=255)

          ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               >  2:1611746/1‑93 (MQ=255)

            aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  2:1866680/93‑1 (MQ=255)

            aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  1:499485/93‑1 (MQ=255)

             cTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTg                                                                            <  2:4020355/93‑1 (MQ=255)

              tCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTgg                                                                           <  1:2887969/93‑1 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     >  2:3892820/1‑93 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     >  1:1142359/1‑93 (MQ=255)

                     ttCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAtt                                                                    >  2:1440295/1‑93 (MQ=255)

                     ttCCCACATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAtt                                                                    <  2:6131859/93‑1 (MQ=255)

                      tCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTa                                                                   <  2:1053113/93‑1 (MQ=255)

                      tCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTa                                                                   >  1:4635903/1‑93 (MQ=255)

                      tCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTa                                                                   <  1:4246130/93‑1 (MQ=255)

                       cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  <  2:432783/93‑1 (MQ=255)

                        ccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtg                                                                 >  2:1662601/1‑93 (MQ=255)

                         cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                >  2:2564754/1‑93 (MQ=255)

                         cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                <  2:5901797/93‑1 (MQ=255)

                         cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                <  2:5600109/93‑1 (MQ=255)

                          gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgat                                                               >  2:3543701/1‑93 (MQ=255)

                          gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgat                                                               >  1:3224058/1‑93 (MQ=255)

                          gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgat                                                               >  1:928740/1‑93 (MQ=255)

                            aTCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattg                                                             >  1:2766902/1‑93 (MQ=255)

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            >  1:5844733/1‑93 (MQ=255)

                              cAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattgat                                                           <  1:4468750/93‑1 (MQ=255)

                              cAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattgat                                                           <  1:5774045/93‑1 (MQ=255)

                                gAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGatat                                                         <  2:4195665/93‑1 (MQ=255)

                                gAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGatat                                                         >  1:1593958/1‑93 (MQ=255)

                                 aaTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGAtata                                                        <  1:5183900/93‑1 (MQ=255)

                                 aaTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGAtata                                                        >  2:4140056/1‑93 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  2:4909014/93‑1 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  2:1742243/93‑1 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  1:2112762/93‑1 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  2:1354873/93‑1 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  1:1928992/93‑1 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   >  2:459506/1‑93 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  1:5362904/93‑1 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   >  1:5921641/1‑93 (MQ=255)

                                       gTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAc                                                  >  2:1211194/1‑93 (MQ=255)

                                       gTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAc                                                  <  1:4760073/93‑1 (MQ=255)

                                        tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 <  1:6174080/93‑1 (MQ=255)

                                         ggTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc                                                <  1:3717714/93‑1 (MQ=255)

                                         ggTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc                                                >  1:4042564/1‑93 (MQ=255)

                                          gtagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGcc                                               >  1:3321915/1‑93 (MQ=255)

                                          gtagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGcc                                               <  2:1519390/93‑1 (MQ=255)

                                          gtagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGcc                                               <  2:2928272/93‑1 (MQ=255)

                                            agtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTa                                             >  2:950581/1‑93 (MQ=255)

                                            agtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTa                                             <  1:839910/93‑1 (MQ=255)

                                              taCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAg                                           <  2:244781/93‑1 (MQ=255)

                                              taCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAg                                           >  2:2748368/1‑93 (MQ=255)

                                                 gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        <  2:4446939/93‑1 (MQ=255)

                                                 gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        <  1:345878/93‑1 (MQ=255)

                                                 gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        >  2:5748274/1‑93 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       >  1:1902129/1‑93 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       >  1:6217443/1‑93 (MQ=255)

                                                   tGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTg                                      <  2:5865/93‑1 (MQ=255)

                                                       ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  >  2:2302236/1‑93 (MQ=255)

                                                       ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  <  1:4910834/93‑1 (MQ=255)

                                                       ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  <  2:2369073/93‑1 (MQ=255)

                                                        cTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACAt                                 <  2:5534730/93‑1 (MQ=255)

                                                         ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                <  1:2287736/93‑1 (MQ=255)

                                                           gcgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGt                              <  1:2691363/93‑1 (MQ=255)

                                                           gcgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGt                              <  2:5054589/93‑1 (MQ=255)

                                                             gcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTcc                            >  2:3354066/1‑93 (MQ=255)

                                                             gcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTcc                            <  1:5652760/93‑1 (MQ=255)

                                                              cATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCa                           <  2:5426667/93‑1 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCata                         >  2:4792671/1‑93 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCata                         <  2:4403178/93‑1 (MQ=255)

                                                                 cTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCatat                        >  1:6136066/1‑93 (MQ=255)

                                                                 cTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCatat                        <  1:603642/93‑1 (MQ=255)

                                                                    tttCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTc                     <  2:3944449/93‑1 (MQ=255)

                                                                       cttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTgg                  >  1:5440114/1‑93 (MQ=255)

                                                                        ttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGc                 >  2:5110194/1‑93 (MQ=255)

                                                                         tcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGcc                >  2:2680390/1‑93 (MQ=255)

                                                                         tcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGcc                >  1:2523805/1‑93 (MQ=255)

                                                                          cATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCt               >  2:62707/1‑93 (MQ=255)

                                                                           aTTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtg              <  2:1962470/93‑1 (MQ=255)

                                                                            tttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtgg             <  1:1261215/93‑1 (MQ=255)

                                                                             ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            <  1:3857845/93‑1 (MQ=255)

                                                                              tAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtg           <  1:5264540/93‑1 (MQ=255)

                                                                                  gCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtta       >  2:3196838/1‑93 (MQ=255)

                                                                                  gCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtta       >  1:1579401/1‑93 (MQ=255)

                                                                                  gCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtta       >  1:995458/1‑93 (MQ=255)

                                                                                    cTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatt     >  1:2200475/1‑93 (MQ=255)

                                                                                      gACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTaa   >  2:3276328/1‑93 (MQ=255)

                                                                                      gACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTaa   <  1:3572525/93‑1 (MQ=255)

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

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