BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,595,233	A→T	5.2%	I102F (ATC→TTC) ‡	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,595,233	0	A	T	5.2%	944.8 / 26.3	287	I102F (ATC→TTC) ‡	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (157/113); new base T (7/8); total (166/121)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 295 total aligned reads displayed.

GCTCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTA  >  NC_002947/3595084‑3595382

gCTCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTa                                                                                                                                                       <  1:915737/150‑1 (MQ=255)

    ggCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGa                                                                                                                                                   >  1:3506657/1‑150 (MQ=255)

         ccACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcctgg                                                                                                                                              <  1:55697/150‑1 (MQ=255)

         ccACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctgg                                                                                                                                              <  1:186846/150‑1 (MQ=255)

         ccACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctgg                                                                                                                                              <  1:3107762/150‑1 (MQ=255)

          cACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctggg                                                                                                                                             <  1:6371027/150‑1 (MQ=255)

           aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcctgggt                                                                                                                                            <  1:4279040/150‑2 (MQ=255)

           aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcctgggc                                                                                                                                            >  2:6081506/1‑150 (MQ=255)

           aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctgggt                                                                                                                                            <  1:9070786/150‑2 (MQ=255)

            cGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcctgggcc                                                                                                                                           <  2:3064481/150‑1 (MQ=255)

                 cTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACctgggcctgggc                                                                                                                                      <  2:7741891/150‑1 (MQ=255)

                  tGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCtgggcctgggcc                                                                                                                                     >  2:8581063/1‑150 (MQ=255)

                  tGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCtgggcctgggcc                                                                                                                                     >  2:7342444/1‑150 (MQ=255)

                   ggTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGACGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCg                                                                                                                                    >  1:5519578/1‑150 (MQ=255)

                         cAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGcc                                                                                                                              >  1:6682097/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  1:791418/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  2:4126872/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  1:121439/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             <  2:1350110/150‑1 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  2:6826179/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  1:3902641/1‑150 (MQ=255)

                          aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCa                                                                                                                             >  1:3780193/1‑150 (MQ=255)

                           aCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAg                                                                                                                            >  1:1432226/1‑150 (MQ=255)

                           aCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAg                                                                                                                            >  1:7125411/1‑150 (MQ=255)

                            cGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAgg                                                                                                                           <  2:5364306/150‑1 (MQ=255)

                              aaaTCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCa                                                                                                                         >  1:367662/1‑150 (MQ=255)

                              aaaTCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCa                                                                                                                         >  2:5825948/1‑150 (MQ=255)

                              aaaTCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGGCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGcc                                                                                                                         >  1:4939629/1‑149 (MQ=255)

                                 tCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGcccag                                                                                                                      <  1:6941616/150‑5 (MQ=255)

                                 tCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGt                                                                                                                      >  1:2502748/1‑150 (MQ=255)

                                  cGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGGTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTc                                                                                                                     >  2:7243483/1‑150 (MQ=255)

                                      aaCGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCggg                                                                                                                 >  1:4547794/1‑150 (MQ=255)

                                           cACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGAc                                                                                                            >  1:6594370/1‑150 (MQ=255)

                                            aCCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGAcc                                                                                                           >  1:5589461/1‑150 (MQ=255)

                                               cGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGtt                                                                                                        <  2:5086152/150‑1 (MQ=255)

                                               cGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGtt                                                                                                        <  2:8504990/150‑1 (MQ=255)

                                               cGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGtt                                                                                                        >  1:391309/1‑150 (MQ=255)

                                                  cacaTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCgg                                                                                                     <  1:6148489/150‑1 (MQ=255)

                                                   acaTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCggg                                                                                                    >  2:8415269/1‑150 (MQ=255)

                                                     aTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCa                                                                                                  >  2:826525/1‑150 (MQ=255)

                                                     aTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCa                                                                                                  >  1:7676126/1‑150 (MQ=255)

                                                       tCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAtt                                                                                                <  2:9029008/150‑1 (MQ=255)

                                                       tCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAtt                                                                                                >  2:2924618/1‑150 (MQ=255)

                                                       gcGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAtt                                                                                                <  2:6962401/149‑1 (MQ=255)

                                                         gtgtGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTgc                                                                                              >  1:4150425/1‑150 (MQ=255)

                                                         gtgtGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTgc                                                                                              >  1:460249/1‑150 (MQ=255)

                                                            tgGGTGCTGGCGGTATCATGCTGCCCAACCTTGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGcgcg                                                                                           >  2:5492546/1‑150 (MQ=255)

                                                                   tGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgc                                                                                    >  2:7929936/1‑150 (MQ=255)

                                                                     gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgcgc                                                                                  <  2:4351960/150‑1 (MQ=255)

                                                                      cGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgcgcc                                                                                 <  2:4497216/150‑1 (MQ=255)

                                                                        gTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCag                                                                               >  1:2845396/1‑150 (MQ=255)

                                                                         tATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCaga                                                                              <  1:2753834/150‑1 (MQ=255)

                                                                           tCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGc                                                                            <  2:16199/150‑1 (MQ=255)

                                                                             aTGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCgg                                                                          <  1:1620314/150‑1 (MQ=255)

                                                                              tgctgcCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCggg                                                                         >  2:3633889/1‑150 (MQ=255)

                                                                               gctgcCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGa                                                                        >  1:8395094/1‑150 (MQ=255)

                                                                               gctgcCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGa                                                                        <  1:6443966/150‑1 (MQ=255)

                                                                               gctgcCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGa                                                                        <  1:7565449/150‑1 (MQ=255)

                                                                                   cccAACCCTGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCCGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                    <  2:4861390/150‑1 (MQ=255)

                                                                                   cccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                    >  2:1734166/1‑150 (MQ=255)

                                                                                   cccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                    >  1:5693912/1‑150 (MQ=255)

                                                                                   cccAACCATGCGCCGTTGGTGGTGGCCGAACAGGTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                    >  1:5390097/1‑150 (MQ=255)

                                                                                   ccaaaCAATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCTACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                    <  1:8197101/147‑1 (MQ=255)

                                                                                    ccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCACTGCGCGGCGCCGCGCCAGAGCGGGATTTCg                                                                   <  1:7985907/150‑1 (MQ=255)

                                                                                       aCCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGcc                                                                <  1:5264090/150‑1 (MQ=255)

                                                                                        ccATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCa                                                               <  1:8927455/150‑1 (MQ=255)

                                                                                             cgcCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACa                                                          <  1:7754868/150‑1 (MQ=255)

                                                                                              gcCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGGTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACAt                                                         >  2:2628401/1‑150 (MQ=255)

                                                                                                  ttGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGcc                                                     <  2:6030071/150‑1 (MQ=255)

                                                                                                  gtggtggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGcc                                                     <  2:9046394/149‑1 (MQ=255)

                                                                                                     gtggtggTCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAg                                                  >  1:2088175/1‑150 (MQ=255)

                                                                                                         gggCCGAACAGTTTGGCACCCTCAATACCTTGCCCCCAGGTCGGACCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGc                                              <  1:2250109/149‑1 (MQ=255)

                                                                                                           gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCTTTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                            <  1:5937405/150‑1 (MQ=255)

                                                                                                                 cAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAc                                      <  2:5195861/150‑1 (MQ=255)

                                                                                                                 cAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAc                                      <  2:6594370/150‑1 (MQ=255)

                                                                                                                  aGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAcc                                     >  1:4765479/1‑150 (MQ=255)

                                                                                                                    tttGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGacaaccgc                                   >  2:572410/1‑143 (MQ=255)

                                                                                                                    tttGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTg                                   >  1:2928463/1‑150 (MQ=255)

                                                                                                                      tGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGc                                 <  2:8808798/150‑1 (MQ=255)

                                                                                                                       ggCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGcc                                <  2:3243072/150‑1 (MQ=255)

                                                                                                                         cACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGa                              >  2:220286/1‑150 (MQ=255)

                                                                                                                             cTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAAc                          >  1:4477960/1‑150 (MQ=255)

                                                                                                                              tCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGa                                          <  2:5553596/133‑1 (MQ=255)

                                                                                                                              tCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg                         <  2:5953185/150‑1 (MQ=255)

                                                                                                                                aaTACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGc                       >  1:1687965/1‑150 (MQ=255)

                                                                                                                                aaTACATTGCACCCAGGTCGTAGCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGc                       <  1:8895273/150‑1 (MQ=255)

                                                                                                                                  tACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCaa                     <  1:8105225/150‑1 (MQ=255)

                                                                                                                                  tACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCaa                     >  2:4663081/1‑150 (MQ=255)

                                                                                                                                      ttGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGc                 >  1:9119290/1‑150 (MQ=255)

                                                                                                                                       tGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGcc                >  1:7092102/1‑150 (MQ=255)

                                                                                                                                        gCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCg               <  2:480239/150‑1 (MQ=255)

                                                                                                                                         cACCCTGGTCGTATCGACGTGGGCCTGGGCCGTGCGCCAGGCCCGTCGGGGCCGACCGTTGGGGCTTTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCgg              <  1:17829/150‑1 (MQ=255)

                                                                                                                                         cACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCgg              <  2:758976/150‑1 (MQ=255)

                                                                                                                                           cccAGGTCGATACGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt            >  2:7172552/1‑150 (MQ=255)

                                                                                                                                            ccAGGTCGATACGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTg           >  1:3457771/1‑150 (MQ=255)

                                                                                                                                                 tCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcgg      >  2:8984233/1‑150 (MQ=255)

                                                                                                                                                  cGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggc     <  2:8002620/150‑1 (MQ=255)

                                                                                                                                                    tATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGa                                                   <  1:7223629/102‑1 (MQ=255)

                                                                                                                                                     aTCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTa  >  1:4020625/1‑150 (MQ=255)

                                                                                                                                                     aTCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTa  >  2:1851744/1‑150 (MQ=255)

GCTCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTA  >  NC_002947/3595084‑3595382

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