BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,522,275	T→G	H607P (CAT→CCT)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,522,275	0	T	G	98.9%	607.5 / ‑3.0	178	H607P (CAT→CCT)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base G (87/88); total (90/88)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 178 total aligned reads displayed.

TAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCT  >  NC_000913/4522179‑4522371

 aaaTACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCt                                                                                              >  1:827671/1‑100 (MQ=255)

 aaaTACAACCGGTCCATGGCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCt                                                                                              <  2:3225727/100‑1 (MQ=255)

  aaTACAACCGGTCCATGTCCCCCCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTg                                                                                             >  1:3482414/1‑100 (MQ=255)

    tACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTggg                                                                                           >  2:1213640/1‑100 (MQ=255)

      cAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTc                                                                                         >  1:1669062/1‑100 (MQ=255)

        aCCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTg                                                                                       >  2:355895/1‑100 (MQ=255)

        aCCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGCGTCTg                                                                                       <  1:2339509/100‑1 (MQ=255)

         ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTgg                                                                                      <  1:1957050/100‑1 (MQ=255)

         ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTgg                                                                                      <  1:2185585/100‑1 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    <  2:2815506/100‑1 (MQ=255)

              ccATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCGTCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                 <  1:3887762/100‑1 (MQ=255)

                aTGTCCCACCGCTCACGGCGTGAAGAATAGCCCATTGGCAGGTATCATCTGGAACTTTAGGATCGGCAAGTAAATCCTGAGGGTTGGGTCTGGCATTTaa                                                                               <  2:798807/100‑1 (MQ=255)

                 tGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAat                                                                              >  1:785568/1‑100 (MQ=255)

                 tGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGATGGGTCTGGCATTTAat                                                                              >  1:166243/1‑100 (MQ=255)

                    cccACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGACGTAATTCCTGAGGGCTGGGTCTGGCATTTAatata                                                                           >  2:1559124/1‑99 (MQ=255)

                    cccACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                           >  1:1110361/1‑100 (MQ=255)

                     ncaccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTACTATTg                                                                          >  1:1889537/2‑100 (MQ=255)

                     ccaccaCTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTg                                                                          >  2:238505/1‑100 (MQ=255)

                      caccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGc                                                                         <  2:1139295/100‑1 (MQ=255)

                      caccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGc                                                                         <  2:3501180/100‑1 (MQ=255)

                       accacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCa                                                                        <  1:1401594/100‑1 (MQ=255)

                        ccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAg                                                                       <  1:1048479/100‑1 (MQ=255)

                        ccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAg                                                                       <  2:450499/100‑1 (MQ=255)

                          acTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGtt                                                                     <  2:2551595/100‑1 (MQ=255)

                          acTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCT‑GCATTTAATATTGCAGTTg                                                                    <  1:763890/100‑1 (MQ=255)

                           cTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTg                                                                    <  2:3373804/100‑1 (MQ=255)

                             cACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcc                                                                  <  2:3773487/100‑1 (MQ=255)

                             cACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcc                                                                  >  2:3426460/1‑100 (MQ=255)

                             cACGGCCTGAAGCATTGCCCATAACCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcc                                                                  >  2:3958296/1‑100 (MQ=255)

                              aCGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcct                                                                 >  1:1470084/1‑100 (MQ=255)

                                 gCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcct                                                              >  2:1757731/1‑100 (MQ=255)

                                  ccTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcctc                                                             >  1:2000216/1‑100 (MQ=255)

                                  ccTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcctc                                                             >  1:393215/1‑100 (MQ=255)

                                    tGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcctcct                                                           >  1:235864/1‑100 (MQ=255)

                                    tGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcctcct                                                           >  1:1766544/1‑100 (MQ=255)

                                        gCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGa                                                       <  1:2384760/100‑1 (MQ=255)

                                        gCATTGCACATAGCCGTGTATCATCAGGCAATTCAGGATCGGGAAGTAAACCCTTAGGGCTGTATCTGGCATTTAATCTTGCAGTTGCCTCCTCCGGTGa                                                       >  2:1271720/1‑100 (MQ=255)

                                         cATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGaa                                                      >  2:861572/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:2191847/1‑100 (MQ=255)

                                              cccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTaa                                                 >  1:360752/1‑100 (MQ=255)

                                               ccATAGCCGGGTATCATCTGGCAATTCCGGATCGGGAGGTAAATCCTGAGGCCCGGGTCTGGCATTTAATAGTGCAGTTGCCTCCTCCCGTGAAGGTAAt                                                >  2:1952156/1‑100 (MQ=255)

                                               ccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAAt                                                >  2:2497638/1‑100 (MQ=255)

                                                cATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAAtg                                               >  1:2016768/1‑99 (MQ=255)

                                                 aTAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTg                                              <  1:3507520/100‑1 (MQ=255)

                                                 aTAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTg                                              >  1:2629585/1‑100 (MQ=255)

                                                  tAGCGGGGTATCATCTGGCAATTCAGGATCGAGAAGTAAATCCTGAGTGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGt                                             <  2:707783/100‑1 (MQ=255)

                                                  tAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGt                                             <  1:2908449/100‑1 (MQ=255)

                                                    gCCGGGTATCATCTGGCAATTCAGGATCTGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTc                                           >  2:2438750/1‑100 (MQ=255)

                                                        ggTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCg                                       >  1:2858628/1‑100 (MQ=255)

                                                         gTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGc                                      <  2:3409888/100‑1 (MQ=255)

                                                            tcatcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTAc                                   <  2:769093/100‑1 (MQ=255)

                                                            tcatcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTAc                                   <  2:1365362/100‑1 (MQ=255)

                                                             catcTGGGAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACg                                  <  2:3447719/100‑1 (MQ=255)

                                                             catcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACg                                  >  2:2155731/1‑100 (MQ=255)

                                                             catcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATAGCTACg                                  >  2:3305377/1‑100 (MQ=255)

                                                              atcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACgg                                 <  2:3920577/100‑1 (MQ=255)

                                                              atcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACgg                                 >  1:1285806/1‑100 (MQ=255)

                                                                cTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATCGTTCATCGCTACGGGt                               >  1:1262654/1‑100 (MQ=255)

                                                                 tGGCAATTCAGGATCGGGAAGTAAATCCTGACGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGtt                              >  1:3271533/1‑100 (MQ=255)

                                                                  ggCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTc                             >  2:2281114/1‑100 (MQ=255)

                                                                  ggCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTc                             >  2:1495278/1‑100 (MQ=255)

                                                                    cAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCg                           >  1:1447886/1‑100 (MQ=255)

                                                                    cAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCg                           >  2:1741364/1‑100 (MQ=255)

                                                                    cAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCg                           >  1:1602579/1‑100 (MQ=255)

                                                                     aaTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGa                          >  2:3047657/1‑100 (MQ=255)

                                                                      aTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAg                         <  1:2700052/100‑1 (MQ=255)

                                                                      aTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAg                         >  1:3367134/1‑100 (MQ=255)

                                                                      aTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAg                         <  1:3457083/100‑1 (MQ=255)

                                                                       ttCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAgg                        <  2:2468124/100‑1 (MQ=255)

                                                                       ttCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAgg                        <  2:2972215/100‑1 (MQ=255)

                                                                       ttCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAgg                        <  1:4027320/100‑1 (MQ=255)

                                                                         caagATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTCATATTGCAGTTGCCTCCTCCTGTGAAGCTAATTGTTCGTCGTTACTAGTTCCGAggta                      >  2:2186591/4‑98 (MQ=255)

                                                                         cAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTACTTGTTCATCGCTACGGGTTCCGAGGaa                      >  1:4076654/1‑100 (MQ=255)

                                                                          aGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGaaa                     <  1:3702775/100‑1 (MQ=255)

                                                                           ggATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAAt                    >  1:88459/1‑100 (MQ=255)

                                                                            gATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAAtt                   >  1:49009/1‑100 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  1:206012/100‑1 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  1:2241088/100‑1 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGGTGGCTCCCCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  2:512552/100‑1 (MQ=255)

                                                                              tCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGa                 <  1:1989510/100‑1 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               >  2:2676112/1‑100 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               >  1:3425587/1‑100 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               <  2:108696/100‑1 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               <  2:158610/100‑1 (MQ=255)

                                                                                 ggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACtt              >  1:690974/1‑100 (MQ=255)

                                                                                 ggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACtt              >  1:2223903/1‑100 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          >  2:2036794/1‑100 (MQ=255)

                                                                                      tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCg         <  1:3093772/100‑1 (MQ=255)

                                                                                      tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCg         >  1:248882/1‑100 (MQ=255)

                                                                                      tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCg         <  2:3108005/100‑1 (MQ=255)

                                                                                      tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCg         >  1:2990365/1‑100 (MQ=255)

                                                                                       aaaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGt        <  1:1323307/100‑1 (MQ=255)

                                                                                        aaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTg       <  1:2717407/100‑1 (MQ=255)

                                                                                        aaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTg       <  1:3451162/100‑1 (MQ=255)

                                                                                        aaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTg       <  1:952221/100‑1 (MQ=255)

                                                                                           ccTGAGTGCTGCGTCTGGCATTTAATATTGCAGTTGCCTCCGCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGTAATTGACTTCGCCGTGAAg    <  2:3655026/100‑1 (MQ=255)

                                                                                           ccTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAGTTGTTCATCGCTCCGGGTTCCGAGGAAATTGACTTCGCCGTGAAg    >  2:3382746/1‑100 (MQ=255)

                                                                                           ccTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAg    <  2:3282433/100‑1 (MQ=255)

                                                                                            cTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGc   >  1:647025/1‑100 (MQ=255)

                                                                                             tGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCt  >  2:1010070/1‑100 (MQ=255)

TAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCT  >  NC_000913/4522179‑4522371

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: `‑`