BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,380,790	T→C	5.2%	A137A (GCT→GCC)	oprF →	porin F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,380,790	0	T	C	5.2%	285.9 / 3.1	96	A137A (GCT→GCC)	oprF	porin F
Reads supporting (aligned to +/- strand): ref base T (46/44); new base C (2/3); total (48/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 105 total aligned reads displayed.

TCGAACACTGCCCTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAGCTGCTCC  >  NC_002947/2380644‑2380939

tCGAACACTGCCCTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTtggcgctggc                                                                                                                                                     <  2:1811219/150‑1 (MQ=255)

tCGAACACTGCCCTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTtggcgctggc                                                                                                                                                     >  1:286308/1‑150 (MQ=255)

     cacTGCCCTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCCAGCGCCaa                                                                                                                                                <  1:497748/150‑1 (MQ=255)

           ccTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTa                                                                                                                                          <  2:1536722/150‑1 (MQ=255)

           ccTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTa                                                                                                                                          <  1:1657903/150‑1 (MQ=255)

            cTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAc                                                                                                                                         >  1:2137485/1‑150 (MQ=255)

            cTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAc                                                                                                                                         <  2:1485235/150‑1 (MQ=255)

               gACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGt                                                                                                                                           <  1:2282133/145‑1 (MQ=255)

                aCGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCa                                                                                                                                    >  2:1079116/1‑150 (MQ=255)

                 cGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATccc                                                                                                                                   <  2:2242810/150‑3 (MQ=255)

                 cGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCAc                                                                                                                                   <  1:2370323/150‑1 (MQ=255)

                 cGCCGTTTACCACCTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCAc                                                                                                                                   <  1:111768/150‑1 (MQ=255)

                  gCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGt                                                                                                                                           >  2:174086/1‑142 (MQ=255)

                  gCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGt                                                                                                                                           <  1:174086/142‑1 (MQ=255)

                  gCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCAcc                                                                                                                                  >  2:1385916/1‑150 (MQ=255)

                  gCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCAcc                                                                                                                                  <  1:1385916/150‑1 (MQ=255)

                    cGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGa                                                                                                                                >  2:40057/1‑150 (MQ=255)

                     gTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGAc                                                                                                                               <  2:2137485/150‑1 (MQ=255)

                         aCCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGt                                                                                                                           >  1:1448844/1‑150 (MQ=255)

                          ccACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCTGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGtt                                                                                                                          >  1:912195/1‑150 (MQ=255)

                           cACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCTCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTc                                                                                                                         >  1:855402/1‑150 (MQ=255)

                            aCTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCt                                                                                                                        >  1:1457956/1‑150 (MQ=255)

                             cTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTg                                                                                                                                             >  1:2149061/1‑129 (MQ=255)

                             cTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTg                                                                                                                                             <  2:2149061/129‑1 (MQ=255)

                              ttCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGa                                                                                                                                >  2:1286777/1‑140 (MQ=255)

                              ttCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGa                                                                                                                                <  1:1286777/140‑1 (MQ=255)

                                 aacaacCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGcc                                                                                                                   >  1:817363/1‑150 (MQ=255)

                                 aacaacCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGcc                                                                                                                   >  2:1340540/1‑150 (MQ=255)

                                   caacCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCg                                                                                                                 >  1:1701997/1‑150 (MQ=255)

                                      cccGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTa                                                                                                                                          <  2:1618055/123‑1 (MQ=255)

                                      cccGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCAACAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTa                                                                                                                                          >  1:1618055/1‑123 (MQ=255)

                                       ccGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCCAGCGCCaacgttgg                                                                                                                                          <  1:1810784/122‑7 (MQ=255)

                                       ccGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCCAGCGCCaacgttgg                                                                                                                                          >  2:1810784/1‑116 (MQ=255)

                                        cGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCg                                                                                                            >  1:442823/1‑150 (MQ=255)

                                                   aTCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTGCGGTCGCTG‑CATGTTCTATGCCCGTGCCGGCGTAGAAGCt                                                                                                 >  2:565223/1‑150 (MQ=255)

                                                   aTCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCt                                                                                                 >  1:392344/1‑150 (MQ=255)

                                                   aTCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCt                                                                                                 >  1:1771230/1‑150 (MQ=255)

                                                   aTCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCt                                                                                                 <  2:60334/150‑1 (MQ=255)

                                                     ccGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCa                                                                                               <  2:1701997/150‑1 (MQ=255)

                                                     ccGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCa                                                                                               <  2:2190148/150‑1 (MQ=255)

                                                          cATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTaca                                                                                          >  2:1554648/1‑150 (MQ=255)

                                                           aTACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAg                                                                                                      >  1:1428848/1‑137 (MQ=255)

                                                           aTACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAg                                                                                                      <  2:1428848/137‑1 (MQ=255)

                                                            tACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTacaac                                                                                        <  2:1905207/150‑1 (MQ=255)

                                                            tACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTacaac                                                                                        >  1:922186/1‑150 (MQ=255)

                                                            tACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTacaac                                                                                        <  2:688279/150‑1 (MQ=255)

                                                                  tCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGAc                                                                                  <  2:2149525/150‑1 (MQ=255)

                                                                       tGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAggg                                                                             >  1:2401862/1‑150 (MQ=255)

                                                                       tGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCCAGCGCCaacgt                                                                                                                                             <  1:345942/87‑4 (MQ=255)

                                                                       tGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCCAGCGCCaacgt                                                                                                                                             >  2:345942/1‑84 (MQ=255)

                                                                        ggTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGc                                                                            <  2:2073123/150‑1 (MQ=255)

                                                                         gTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCCTGACCACTCCCCCTTCGCCAACGATGGCACTGGCGCCAATTGGTAC‑ATCCCCGACATGTTCTATCCCCGTACAGGCGTAGAAGCTCAGTACAACATCGACCAGGGCg                                                                           >  2:1049554/1‑150 (MQ=255)

                                                                            tctcGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGaca                                                                        <  1:449515/150‑1 (MQ=255)

                                                                            ctgtgCACCAATCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGaca                                                                        >  2:297968/5‑150 (MQ=255)

                                                                                gCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGa                                                                    >  1:319940/1‑150 (MQ=255)

                                                                                   ccAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTg                                                                 <  2:77839/150‑1 (MQ=255)

                                                                                   ccAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCAACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTg                                                                 <  1:818447/150‑1 (MQ=255)

                                                                                   acAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTg                                                                 >  1:1219098/2‑150 (MQ=255)

                                                                                    cAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTgg                                                                <  1:1018648/150‑1 (MQ=255)

                                                                                       ccGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCCCCCCCCCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGcc                                                             <  1:725434/149‑1 (MQ=255)

                                                                                             ttccATACCGGCCGTGGCGGTCGTGACCCCTCCACCTTCGCCCACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGc                                                       <  2:392344/148‑1 (MQ=255)

                                                                                              gCCAGACCGGCCGTGGCGATCGTGACCACTCCACCTTCGCCGACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCCGGGCGACACCGAGCGGGCCCCgctcg                                                      >  2:630040/1‑146 (MQ=255)

                                                                                               ccAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGt                                                     >  1:980086/1‑150 (MQ=255)

                                                                                                 aGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCACCGCGCGTTg                                                   <  2:1913968/150‑1 (MQ=255)

                                                                                                     cGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTCGGCCCCGCGCGTTGgtgt                                               >  2:369187/1‑150 (MQ=255)

                                                                                                       gCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGCTg                                             <  1:1340540/150‑1 (MQ=255)

                                                                                                         cGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGgtgt                                               >  2:80658/1‑146 (MQ=255)

                                                                                                         cGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGgtgt                                               <  1:80658/146‑1 (MQ=255)

                                                                                                           tGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCt                                         <  1:1554648/150‑1 (MQ=255)

                                                                                                           tGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCt                                         >  2:1968362/1‑150 (MQ=255)

                                                                                                           tGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCt                                         <  2:286308/150‑1 (MQ=255)

                                                                                                                  cGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTc                                  >  2:383211/1‑150 (MQ=255)

                                                                                                                   gTGACCACTCCACCTTCGCCAACGTTGGCGCGGGGGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTcg                                 <  1:635889/150‑1 (MQ=255)

                                                                                                                      aCCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTcggcg                              >  2:346482/1‑150 (MQ=255)

                                                                                                                        cACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGt                            <  2:2186998/150‑1 (MQ=255)

                                                                                                                           tCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGt                                                                  >  1:268875/1‑109 (MQ=255)

                                                                                                                           tCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGt                                                                  <  2:268875/109‑1 (MQ=255)

                                                                                                                            ccACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTAAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTcgg                                <  2:270517/142‑1 (MQ=255)

                                                                                                                            ccACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTAAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTcgg                                >  1:270517/1‑142 (MQ=255)

                                                                                                                              aCCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCg                      >  2:363345/1‑150 (MQ=255)

                                                                                                                                ctccGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTTCAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    <  2:670957/147‑1 (MQ=255)

                                                                                                                                cTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    <  1:1049554/150‑1 (MQ=255)

                                                                                                                                cTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    >  1:1908164/1‑150 (MQ=255)

                                                                                                                                 ttCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGGCTGAACTTCGGCGGTAGCCCGaag                   >  2:1211333/1‑150 (MQ=255)

                                                                                                                                 ttCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    >  2:2099861/1‑149 (MQ=255)

                                                                                                                                 ttCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    <  1:2099861/149‑1 (MQ=255)

                                                                                                                                  tCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaagc                  <  2:1771230/150‑1 (MQ=255)

                                                                                                                                     ccAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaa                    >  1:834954/1‑145 (MQ=255)

                                                                                                                                      cAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGaagcaagc              >  2:1198316/1‑150 (MQ=255)

                                                                                                                                       aaCGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCaagct             <  1:2092484/150‑1 (MQ=255)

                                                                                                                                         cGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCaagctga           <  2:699763/150‑1 (MQ=255)

                                                                                                                                           ttGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGccc                                                                                                                  <  2:2094869/45‑1 (MQ=255)

                                                                                                                                           ttGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGccc                                                                                                                  >  1:2094869/1‑45 (MQ=255)

                                                                                                                                               cgctggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAgctgc     <  1:264634/150‑1 (MQ=255)

                                                                                                                                               cgctggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGACCCTACCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAgctgc     <  1:1356412/150‑1 (MQ=255)

                                                                                                                                                gctggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAgctgct    <  2:1562908/150‑1 (MQ=255)

                                                                                                                                                 ctggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTg                                                                 >  2:1740704/1‑88 (MQ=255)

                                                                                                                                                 ctggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTg                                                                 <  1:1740704/88‑1 (MQ=255)

                                                                                                                                                  tggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAGctgctcc  >  2:2071449/1‑150 (MQ=255)

                                                                                                                                                  tggcgcCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAGctgctcc  <  2:319940/150‑1 (MQ=255)

TCGAACACTGCCCTGGACGCCGTTTACCACTTCAACAACCCGTACGACGCTATCCGTCCATACGTTTCCGCTGGTTTCTCGCACCAGTCGCTGGGCCAGACCGGCCGTGGCGGTCGTGACCACTCCACCTTCGCCAACGTTGGCGCTGGCGCCAAGTGGTAC‑ATCACCGACATGTTCTATGCCCGTGCCGGCGTAGAAGCTCAGTACAACATCGACCAGGGCGACACCGAGTGGGCCCCGAGCGTTGGTGTTGGCCTGAACTTCGGCGGTAGCCCGAAGCAAGCTGAAGCTGCTCC  >  NC_002947/2380644‑2380939

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