BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,832	A→G	R146R (CGT→CGC)	yjjT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,832	0	A	G	99.6%	761.7 / ‑4.9	238	R146R (CGT→CGC)	yjjT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (124/113); total (125/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 249 total aligned reads displayed.

AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921

  acacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcgac                                                                                       <  2:5198112/93‑1 (MQ=255)

  acacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAgacgac                                                                                       <  1:1580589/93‑1 (MQ=255)

   cacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGAcc                                                                                      <  1:1053220/93‑1 (MQ=255)

   cacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGAcc                                                                                      <  2:3613326/93‑1 (MQ=255)

      gCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCaaa                                                                                   <  2:3204713/93‑1 (MQ=255)

          ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               <  1:1920148/93‑1 (MQ=255)

          ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               >  1:4140908/1‑93 (MQ=255)

            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACTAAAATagag                                                                             <  2:4423248/93‑1 (MQ=255)

            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             <  2:2589064/93‑1 (MQ=255)

            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             <  2:5317559/93‑1 (MQ=255)

            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             <  2:546303/93‑1 (MQ=255)

             tgCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            <  1:14182/92‑1 (MQ=255)

             aGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            >  1:5209516/1‑93 (MQ=255)

             aGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            >  1:2169390/1‑93 (MQ=255)

              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCACATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           >  1:285545/1‑93 (MQ=255)

              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           <  1:2394755/93‑1 (MQ=255)

              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           <  2:4524257/93‑1 (MQ=255)

              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           <  1:2465974/93‑1 (MQ=255)

                gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         >  1:2615902/1‑93 (MQ=255)

                 ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcac                                                                        >  1:2216547/1‑93 (MQ=255)

                     gACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcg                                                                    >  1:3883860/1‑93 (MQ=255)

                     gACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGACACAGcg                                                                    >  2:4582939/1‑93 (MQ=255)

                      aCCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcga                                                                   <  2:5863676/93‑1 (MQ=255)

                      aCCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcga                                                                   <  1:1195273/93‑1 (MQ=255)

                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  <  2:851899/93‑1 (MQ=255)

                         gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgacga                                                                <  2:2062533/93‑1 (MQ=255)

                          tCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAg                                                               <  1:5071838/93‑1 (MQ=255)

                             ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  2:3026508/93‑1 (MQ=255)

                                cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         <  1:5848811/93‑1 (MQ=255)

                                cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         <  1:893738/93‑1 (MQ=255)

                                cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         <  2:1611192/93‑1 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:4161785/1‑93 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:4000803/1‑93 (MQ=255)

                                  tCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGATGCCACAGCGACGAGCGCTGTCg                                                       >  2:3648934/1‑93 (MQ=255)

                                  tCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCg                                                       >  1:3086469/1‑93 (MQ=255)

                                   cGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCgc                                                      >  2:1918219/1‑92 (MQ=255)

                                   cGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGa                                                      <  2:1685302/93‑1 (MQ=255)

                                     aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  1:2161964/93‑1 (MQ=255)

                                      cGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACtt                                                   >  1:1260846/1‑93 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 <  1:4408117/93‑1 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAACAGAGGCCACAGCGACGAGCGCTGTCGACtttt                                                 >  2:1830486/1‑92 (MQ=255)

                                         tGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAt                                                >  1:2089540/1‑93 (MQ=255)

                                         tGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAt                                                <  2:135863/93‑1 (MQ=255)

                                          gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               >  1:3297503/1‑93 (MQ=255)

                                          gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               <  2:5353467/93‑1 (MQ=255)

                                           tATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTc                                              <  2:4102911/93‑1 (MQ=255)

                                           tATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTc                                              <  1:4225546/93‑1 (MQ=255)

                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             <  2:5266940/93‑1 (MQ=255)

                                             ttCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCaa                                            >  1:4791854/1‑93 (MQ=255)

                                             tccGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCaa                                            <  2:3098536/91‑1 (MQ=255)

                                              tCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAAc                                           <  1:4066626/93‑1 (MQ=255)

                                              tCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAAc                                           <  1:297688/93‑1 (MQ=255)

                                                gCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACgg                                         <  2:5664677/93‑1 (MQ=255)

                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCCACGgc                                        >  2:140183/1‑93 (MQ=255)

                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        >  1:1973717/1‑93 (MQ=255)

                                                   ccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgcgc                                      >  2:5116537/1‑93 (MQ=255)

                                                   ccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgcgc                                      >  2:4737789/1‑93 (MQ=255)

                                                   ccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgcgc                                      >  1:1082269/1‑93 (MQ=255)

                                                    cAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCa                                     >  2:6077232/1‑93 (MQ=255)

                                                     aGAATTTCTCCGCATCAAATACCGGCTGTTTTTTCAGGCGACCAAAATAGAGGCCAAAGCGACGAGCGCTGTCGACTTTATTCAACAGCGCAt                                    >  2:3543429/1‑93 (MQ=255)

                                                     aGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    >  2:176109/1‑93 (MQ=255)

                                                      gAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATa                                   <  2:1108676/93‑1 (MQ=255)

                                                      gAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATa                                   <  1:4648729/93‑1 (MQ=255)

                                                      gAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATa                                   <  1:3532331/93‑1 (MQ=255)

                                                       aaTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATaa                                  >  1:2438472/1‑93 (MQ=255)

                                                        aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 <  1:1021429/93‑1 (MQ=255)

                                                        aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 <  1:6012882/93‑1 (MQ=255)

                                                        aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 >  1:5570738/1‑93 (MQ=255)

                                                        aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 <  1:4871849/93‑1 (MQ=255)

                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  2:1137778/1‑93 (MQ=255)

                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                <  2:216659/93‑1 (MQ=255)

                                                             tcCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGACACAGCGACGAGCAGCGTCGACTTTATTCAACGGCGTATACTATGcc                            >  2:1211639/1‑93 (MQ=255)

                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  2:904718/93‑1 (MQ=255)

                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           >  1:5805866/1‑93 (MQ=255)

                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           >  2:5048931/1‑93 (MQ=255)

                                                               cGCATCAAATACCGGCTGTTTTTCCTGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAg                          <  1:6086488/93‑1 (MQ=255)

                                                                gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         >  2:2970557/1‑93 (MQ=255)

                                                                gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         <  1:5382782/93‑1 (MQ=255)

                                                                 cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        >  1:6032206/1‑93 (MQ=255)

                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:3339927/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:83798/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      <  1:6298281/93‑1 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCAAGCATc                      >  2:2694256/1‑93 (MQ=255)

                                                                       aTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCt                  >  1:95991/1‑93 (MQ=255)

                                                                       aTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCt                  <  2:2740396/93‑1 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  2:1552515/93‑1 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 >  2:4213255/1‑93 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  2:6099342/93‑1 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 >  1:5674443/1‑93 (MQ=255)

                                                                        tACCGGCTGTGTTGCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  1:729309/93‑1 (MQ=255)

                                                                           cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              <  2:3865552/93‑1 (MQ=255)

                                                                             gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            >  2:2914123/1‑93 (MQ=255)

                                                                             gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            >  2:3698928/1‑93 (MQ=255)

                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:5107288/1‑93 (MQ=255)

                                                                                gTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTgc         >  1:3456605/1‑93 (MQ=255)

                                                                                gTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTgc         >  2:4958849/1‑93 (MQ=255)

                                                                                  ttttCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTgcgc       >  2:2474940/1‑93 (MQ=255)

                                                                                    ttCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAc     <  1:2263280/93‑1 (MQ=255)

                                                                                      ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   <  1:1849149/93‑1 (MQ=255)

                                                                                      ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   >  2:4889264/1‑93 (MQ=255)

AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921

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