BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,929	T→G	L411R (CTG→CGG)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,929	0	T	G	99.3%	469.7 / ‑4.2	136	L411R (CTG→CGG)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (74/61); total (74/62)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 136 total aligned reads displayed.

CAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGG  >  NC_000913/1756832‑1757024

cAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGa                                                                                               >  1:406721/1‑100 (MQ=255)

cAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGAAAAAAAAAAAAAAAACGGCTCATCAGTTGGTACGGa                                                                                               <  2:3139035/100‑1 (MQ=255)

 aGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAg                                                                                              <  2:2055931/100‑1 (MQ=255)

  gggCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGc                                                                                             >  2:1547677/1‑100 (MQ=255)

  gggCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGc                                                                                             >  2:963162/1‑100 (MQ=255)

     cGGAAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCaaa                                                                                          >  2:3717369/1‑100 (MQ=255)

      ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAg                                                                                         >  1:1600540/1‑100 (MQ=255)

      ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAg                                                                                         <  1:2452384/100‑1 (MQ=255)

         aaaTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCg                                                                                      >  1:1415059/1‑100 (MQ=255)

           aTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGtt                                                                                    >  1:1393639/1‑100 (MQ=255)

            tCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTtg                                                                                   <  1:1867603/100‑1 (MQ=255)

            tCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTtg                                                                                   <  1:2574918/100‑1 (MQ=255)

            tCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTtg                                                                                   >  1:3585102/1‑100 (MQ=255)

            tCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTAAAAACAAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTtg                                                                                   <  1:3930327/100‑1 (MQ=255)

             cTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTtgt                                                                                  <  1:1538072/100‑1 (MQ=255)

               gCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAGAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGc                                                                                <  2:598939/100‑1 (MQ=255)

                cTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACCAACAAACGGCTCACCAGTTGATACGGAGCAAAGGCATTGTAcc                                                                               >  2:1164500/1‑100 (MQ=255)

                  cgcgCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGCGCAAAGGCGTTGTGCCGc                                                                             >  1:4033095/1‑100 (MQ=255)

                  cgcgCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGc                                                                             >  1:137862/1‑100 (MQ=255)

                   gcgcAGTACGTAAATAATGCCCGTATGCAACCCTCCTGGCACAGAACAAAAAAGAAAAAACGGCTCATAAGTTGGTACGGAGCAAAGGCGTTGTGCCGCa                                                                            <  2:387907/100‑1 (MQ=255)

                    agcAGTACGTAAATACTTCCCGGATGCCACCATGCTGGCACTGACCACCAACGAAAAAACGGCTCTTCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAg                                                                           <  2:2309745/99‑1 (MQ=255)

                        gTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTg                                                                       >  2:850516/1‑100 (MQ=255)

                         tACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAGACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGt                                                                      <  2:928912/100‑1 (MQ=255)

                         tACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGt                                                                      <  1:3496361/100‑1 (MQ=255)

                          aCGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGtt                                                                     >  1:3524424/1‑100 (MQ=255)

                           cGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTa                                                                    <  1:2052350/100‑1 (MQ=255)

                            gTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa                                                                   >  1:1127209/1‑100 (MQ=255)

                             tAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTaaa                                                                  <  1:2543570/100‑1 (MQ=255)

                             tAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTCTGCCGCAGCATGTTaaa                                                                  >  1:3757899/1‑100 (MQ=255)

                              aaaTACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCCGTTGGTACGGAGCAAAGGCGTTGTGTCGCAGCTTGTTAAag                                                                 >  2:3478279/1‑100 (MQ=255)

                              aaaTACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                 <  2:2272929/100‑1 (MQ=255)

                               aaTACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAaga                                                                >  2:678593/1‑100 (MQ=255)

                                atccTTCCCGGAGGCCACCAGCCGGCCACTGACCACCAACGAAAATACGTCTCATCAGTTTGTACGGAGCAAAGGCGTAGTGCCCCAGCTTGTTAAagag                                                               <  2:1980800/97‑1 (MQ=255)

                                 tACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAgaga                                                              >  1:2122791/1‑100 (MQ=255)

                                   cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATc                                                            <  2:1116639/100‑1 (MQ=255)

                                   cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATc                                                            <  1:2556801/100‑1 (MQ=255)

                                    ttCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCa                                                           <  2:138446/100‑1 (MQ=255)

                                    ttCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCa                                                           >  1:669988/1‑100 (MQ=255)

                                     tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                          <  2:1040319/100‑1 (MQ=255)

                                     tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                          <  1:4031489/100‑1 (MQ=255)

                                     tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                          >  1:1821482/1‑100 (MQ=255)

                                      cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACt                                                         >  1:3685730/1‑100 (MQ=255)

                                      cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACt                                                         >  1:885108/1‑100 (MQ=255)

                                      cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACt                                                         >  2:405198/1‑100 (MQ=255)

                                      cccGGATGCCACCATCCTGGCACTGAACAACAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACt                                                         <  2:3574/100‑1 (MQ=255)

                                       ccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACtt                                                        <  1:894426/100‑1 (MQ=255)

                                       ccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACtt                                                        <  1:569544/100‑1 (MQ=255)

                                       ccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACtt                                                        >  1:275173/1‑100 (MQ=255)

                                        cGGATGCCACCATCCTGGCACAGAAAAAAAAAAAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTc                                                       <  2:2487163/100‑1 (MQ=255)

                                        cGGATGCCACCATCCGGGGACTGACGACCAACGAAAGAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTc                                                       <  2:3430933/100‑1 (MQ=255)

                                         ggATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCt                                                      >  1:785773/1‑100 (MQ=255)

                                          gATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCCTCAGTTGGTACGGCGCAAAGGCGTTGTGCGGCAGCTTGTTAAAGCGATCACTTCTa                                                     >  1:3068157/1‑100 (MQ=255)

                                            tGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACt                                                   <  2:1832308/100‑1 (MQ=255)

                                              ccaccaTCCTGGCACTGACGACCAACGAAAAAACGGCTCCTCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtga                                                 >  2:3688822/1‑100 (MQ=255)

                                              ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtga                                                 >  1:3093338/1‑100 (MQ=255)

                                              ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtga                                                 >  1:2386530/1‑100 (MQ=255)

                                               caccaTCCTGGCACTGACCACCAACGAAACAACGGCTCATGAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTCCTACtgat                                                <  2:3838290/100‑1 (MQ=255)

                                               caccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtgat                                                >  2:3556700/1‑100 (MQ=255)

                                                accaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtgatg                                               <  2:3003369/100‑1 (MQ=255)

                                                accaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCACCTTGTTAAAGAGATCACTTCTACtgatg                                               >  2:1859722/1‑100 (MQ=255)

                                                      cTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCt                                         >  2:1014288/1‑100 (MQ=255)

                                                       tGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTa                                        >  1:1164233/1‑100 (MQ=255)

                                                          cACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCg                                     >  2:2614117/1‑100 (MQ=255)

                                                           aCTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGt                                    >  2:3954174/1‑100 (MQ=255)

                                                            cTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTc                                   >  1:2897396/1‑100 (MQ=255)

                                                             tGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCt                                  >  2:3001781/1‑100 (MQ=255)

                                                             tGACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAACAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCt                                  <  2:769302/100‑1 (MQ=255)

                                                              gACCACCAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTg                                 >  2:768770/1‑100 (MQ=255)

                                                               accaccAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTgg                                <  1:2322885/100‑1 (MQ=255)

                                                               accaccAACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTgg                                <  1:340660/100‑1 (MQ=255)

                                                                ccaccaACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTggg                               <  2:3689966/100‑1 (MQ=255)

                                                                ccaccaACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTggg                               >  1:1489890/1‑100 (MQ=255)

                                                                 caccaACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGt                              <  2:3392861/100‑1 (MQ=255)

                                                                   ccaACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTaa                            >  1:1567724/1‑100 (MQ=255)

                                                                    caACGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTaaa                           >  1:3059969/1‑100 (MQ=255)

                                                                     aaCGAGAACACGGCCCATCAGTTGCTACGGAGCAAAGGCGCTGCGCCGCTGCTTGTGCACGAGATCTCTTCTACTGATGATTTCTACCGTCTGGGTAAAg                          <  2:1610237/100‑1 (MQ=255)

                                                                       cGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGaa                        >  1:3867265/1‑100 (MQ=255)

                                                                        gAAAAAACGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                       <  1:2878829/100‑1 (MQ=255)

                                                                         aaaaaaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACt                      >  1:1153949/1‑100 (MQ=255)

                                                                            aaaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGc                   <  1:1696353/100‑1 (MQ=255)

                                                                            aaaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGc                   <  1:2751722/100‑1 (MQ=255)

                                                                            aaaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGc                   >  2:1503465/1‑100 (MQ=255)

                                                                            aaaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGc                   <  2:1248431/100‑1 (MQ=255)

                                                                              aCGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctc                 <  2:173530/100‑1 (MQ=255)

                                                                               cGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctct                <  2:325841/100‑1 (MQ=255)

                                                                               cGGCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctct                <  2:1672577/100‑1 (MQ=255)

                                                                                gtcTCATCAGGTGAAACGGAACAAAGCCGTTGGGCCGCAGCAAGTGAAATAGATAAAATCTACTGAAGATTTCTACCGTCTGGGTAAAGAACGGGCTCTg               <  2:2402499/98‑1 (MQ=255)

                                                                                ggCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTg               <  1:2133795/100‑1 (MQ=255)

                                                                                 gCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc              <  2:2627017/100‑1 (MQ=255)

                                                                                 gCTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc              >  1:2863816/1‑100 (MQ=255)

                                                                                 gCTCATCAGTTGGTAAGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGCTACATTCTACTGATGATTTCTACGGTATGGGTAAAGCAATGGCTATga              >  2:1998307/1‑99 (MQ=255)

                                                                                  cTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa             >  2:1282127/1‑100 (MQ=255)

                                                                                  cTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa             >  1:674656/1‑100 (MQ=255)

                                                                                  cTCATCAGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCGACTTATGATTTCTACGGTCTGGGTAAAGAACTGGCTCTGCa             >  2:1714607/1‑100 (MQ=255)

                                                                                   tcatcaGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCag            >  2:2483041/1‑100 (MQ=255)

                                                                                    catcaGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCaga           >  2:1081591/1‑100 (MQ=255)

                                                                                      tcaGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc         <  2:2232285/100‑1 (MQ=255)

                                                                                        aGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCgg       >  2:7297/1‑100 (MQ=255)

                                                                                            ggTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTg   >  2:67690/1‑100 (MQ=255)

                                                                                             gTACGGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTgg  <  1:1930471/100‑1 (MQ=255)

CAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGG  >  NC_000913/1756832‑1757024

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