BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,030,839	T→A	5.1%	D379E (GAT→GAA)	PP_4432 →	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,030,839	0	T	A	5.1%	1001.9 / 28.0	294	D379E (GAT→GAA)	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand): ref base T (145/134); new base A (9/6); total (154/140)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 299 total aligned reads displayed.

AGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTT  >  NC_002947/5030690‑5030988

  cTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATAc                                                                                                                                                     >  2:3039819/1‑150 (MQ=255)

   tATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAAca                                                                                                                                                    >  1:5107315/1‑149 (MQ=255)

    aTCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTg                                                                                                                                                   >  2:544094/1‑150 (MQ=255)

     tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                                                                  >  1:6695905/1‑150 (MQ=255)

      ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                                                                 <  2:5604057/150‑1 (MQ=255)

        accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCAACCTGATCACTGATACTGGGGt                                                                                                                                               >  2:2398548/1‑150 (MQ=255)

          cggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCg                                                                                                                                             >  1:1566088/6‑150 (MQ=255)

           ggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCga                                                                                                                                            <  1:8071018/146‑1 (MQ=255)

           aGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCga                                                                                                                                            <  1:8436342/150‑1 (MQ=255)

              ctGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatg                                                                                                                                         <  1:3942309/149‑4 (MQ=255)

                 gggCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGtt                                                                                                                                      <  1:8649828/150‑1 (MQ=255)

                  ggCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTc                                                                                                                                     <  2:6784157/150‑1 (MQ=255)

                   gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                                                                    >  1:548101/1‑150 (MQ=255)

                   gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                                                                    >  1:7213270/1‑150 (MQ=255)

                   gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                                                                    <  2:7316611/150‑1 (MQ=255)

                       gcgcACCCTGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTAGGACCTTCCCCTTCATGCCGGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTc                                                                                                                                <  2:1155778/150‑1 (MQ=255)

                       gcgcACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTc                                                                                                                                <  2:393667/150‑1 (MQ=255)

                          cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                                                                             <  2:8854860/150‑1 (MQ=255)

                           aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                                                                            >  1:8214616/1‑150 (MQ=255)

                           aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                                                                            >  2:9098080/1‑150 (MQ=255)

                             cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCGAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                                                          >  2:4676379/1‑150 (MQ=255)

                             cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTCTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                                                          >  1:6907256/1‑150 (MQ=255)

                             cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                                                          >  1:5518840/1‑150 (MQ=255)

                             cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                                                          >  1:260301/1‑150 (MQ=255)

                                   ccTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCaa                                                                                                                    >  1:7395194/1‑150 (MQ=255)

                                   ccTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCaa                                                                                                                    <  2:6407774/150‑1 (MQ=255)

                                    cTGCGCCCGAGCGACAACAGCGTGCTAAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAg                                                                                                                   >  2:6943189/1‑150 (MQ=255)

                                     tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGc                                                                                                                  <  2:2636366/150‑1 (MQ=255)

                                          ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                                                             >  2:3782811/1‑150 (MQ=255)

                                          ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                                                             <  2:5094827/150‑1 (MQ=255)

                                          ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                                                             <  2:2139831/150‑1 (MQ=255)

                                               cgcctacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATTGATGGATGATTGGGGGCTCGAGATTTCCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                                                        <  1:2821426/145‑1 (MQ=255)

                                               cGACAACAGCGTGCTCACACCTGGTATGACCTTCCACTTCATGCCTGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                                                        <  2:9081429/150‑1 (MQ=255)

                                               cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                                                        <  2:4356175/150‑1 (MQ=255)

                                                gACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAg                                                                                                       >  1:7369364/1‑150 (MQ=255)

                                                      aGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGa                                                                                                 <  1:891958/150‑1 (MQ=255)

                                                         gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                                                                              <  1:4211957/150‑1 (MQ=255)

                                                             tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                                                                          >  2:3321255/1‑150 (MQ=255)

                                                               aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                                                        >  2:8634492/1‑150 (MQ=255)

                                                               aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                                                        <  2:5004261/150‑1 (MQ=255)

                                                               aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                                                        >  1:7126856/1‑150 (MQ=255)

                                                                aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                                                                       <  1:2944172/150‑1 (MQ=255)

                                                                  ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGc                                                                                     >  2:4390807/1‑150 (MQ=255)

                                                                  ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGc                                                                                     >  2:177869/1‑150 (MQ=255)

                                                                  ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGc                                                                                     <  1:4624887/150‑1 (MQ=255)

                                                                     ggTATGACCTTCCACTTCATGTCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTTGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCGAg                                                                                  <  2:1299220/150‑1 (MQ=37)

                                                                     ggTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTTGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCGAg                                                                                  <  2:1092556/150‑1 (MQ=38)

                                                                      gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGt                                                                                 >  2:5799428/1‑150 (MQ=255)

                                                                      gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCAATGATACTGGGGTCTAGACGATCTGCTCGACTCAACGCAAGAATTTCGTCAAGGACTGAAGAAAAGTTAGAAAGt                                                                                 >  1:1966295/1‑150 (MQ=255)

                                                                        aTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGTTACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCg                                                                               >  2:6556738/1‑150 (MQ=255)

                                                                           aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                                                                            <  1:957914/150‑1 (MQ=255)

                                                                             cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaa                                                                                            <  2:6671096/132‑1 (MQ=255)

                                                                             cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                                                                          <  2:4052870/150‑1 (MQ=255)

                                                                             cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                                                                          >  1:6602247/1‑150 (MQ=255)

                                                                                ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGtt                                                                       >  2:232283/1‑150 (MQ=255)

                                                                                 cACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTta                                                                      >  2:2339898/1‑150 (MQ=255)

                                                                                   cgacATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCGGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACGGACGAAAAGTTTGCAAGTCGCTGGAGTtata                                                                    <  2:5955869/147‑1 (MQ=255)

                                                                                   cTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTtata                                                                    <  1:5885918/150‑1 (MQ=255)

                                                                                     tCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGt                                                                  >  1:8175295/1‑150 (MQ=255)

                                                                                        tGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGcc                                                               >  1:2750348/1‑150 (MQ=255)

                                                                                          ccTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCAc                                                             >  1:5333301/1‑150 (MQ=255)

                                                                                            tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg                                                           >  1:6361459/1‑150 (MQ=255)

                                                                                            tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg                                                           >  1:1761254/1‑150 (MQ=255)

                                                                                             ggCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGc                                                          <  1:5799428/150‑1 (MQ=255)

                                                                                                cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTAt                                                       >  2:6190136/1‑150 (MQ=255)

                                                                                                 tATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATc                                                      >  1:7609992/1‑150 (MQ=255)

                                                                                                     gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa                                                  <  1:3356342/150‑1 (MQ=255)

                                                                                                       tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc                                                <  2:3069486/150‑1 (MQ=255)

                                                                                                       tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc                                                >  1:4486136/1‑150 (MQ=255)

                                                                                                         gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTCTCTTCAGCGa                                              >  2:7941339/1‑150 (MQ=255)

                                                                                                         gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGa                                              <  2:8410497/150‑1 (MQ=255)

                                                                                                               tGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt                                        <  1:6252423/150‑1 (MQ=255)

                                                                                                               tGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt                                        <  1:6864098/150‑1 (MQ=255)

                                                                                                                gggggCCCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgtg                                       >  1:6423586/1‑150 (MQ=255)

                                                                                                                 ggggCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGa                                      >  2:672259/1‑150 (MQ=255)

                                                                                                                      tCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTg                                 >  1:7794862/1‑150 (MQ=255)

                                                                                                                      tCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTg                                 >  2:1715027/1‑150 (MQ=255)

                                                                                                                            ttACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaa                           <  2:8204182/150‑1 (MQ=255)

                                                                                                                             tACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                          <  1:3732098/150‑1 (MQ=255)

                                                                                                                             tACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                          <  2:7693526/150‑1 (MQ=255)

                                                                                                                             tACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                          <  2:8073940/150‑1 (MQ=255)

                                                                                                                               ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGc                        <  2:574595/150‑1 (MQ=255)

                                                                                                                               ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGc                        <  1:672259/150‑1 (MQ=255)

                                                                                                                                  agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAcc                     >  2:8032539/1‑150 (MQ=255)

                                                                                                                                  agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAcc                     <  2:158186/150‑1 (MQ=255)

                                                                                                                                     gCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATc                  <  1:7148551/150‑1 (MQ=255)

                                                                                                                                      cATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCa                 <  1:8026473/150‑1 (MQ=255)

                                                                                                                                       aTCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAc                >  2:6228603/1‑150 (MQ=255)

                                                                                                                                        tCCTGATCACTGATACTGGGGTCGAGACGTTCTTCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACt               >  1:2319177/1‑150 (MQ=255)

                                                                                                                                          cTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACttt             <  1:9102486/150‑3 (MQ=255)

                                                                                                                                           tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa            <  2:6907256/150‑1 (MQ=255)

                                                                                                                                           tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa            <  1:7923516/150‑1 (MQ=255)

                                                                                                                                            gATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTc                                                   <  2:6077693/110‑1 (MQ=255)

                                                                                                                                            gATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTc                                                   >  1:6077693/1‑110 (MQ=255)

                                                                                                                                                aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc       >  1:3477861/1‑150 (MQ=255)

                                                                                                                                                aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc       >  1:301709/1‑150 (MQ=255)

                                                                                                                                                 cTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACt      >  1:5220301/1‑150 (MQ=255)

                                                                                                                                                  tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc     >  2:3412554/1‑150 (MQ=255)

                                                                                                                                                  tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc     <  1:4400639/150‑1 (MQ=255)

                                                                                                                                                  tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc     <  1:1081991/150‑1 (MQ=255)

AGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTT  >  NC_002947/5030690‑5030988

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