BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,401,532	T→G	8.1%	R108R (CGA→CGC) ‡	PP_3878 ←	minor capsid protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,401,532	0	T	G	8.1%	476.0 / 25.1	149	R108R (CGA→CGC) ‡	PP_3878	minor capsid protein C
Reads supporting (aligned to +/- strand): ref base T (72/65); new base G (7/5); total (79/70)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 155 total aligned reads displayed.

TTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGCCTGGTA  >  NC_002947/4401383‑4401681

ttAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGgctgct                                                                                                                                                       >  2:2882406/1‑150 (MQ=255)

ttAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGgctgct                                                                                                                                                       >  1:711587/1‑150 (MQ=255)

     cAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTAAGGGCTGCGAGCAg                                                                                                                                                  >  2:881817/1‑150 (MQ=255)

         ccGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTAAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCt                                                                                                                                              <  1:3568727/150‑1 (MQ=255)

          cGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCtt                                                                                                                                             <  2:2868817/150‑1 (MQ=255)

              gATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATa                                                                                                                                         >  1:541229/1‑150 (MQ=255)

               aTTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTTAGGGCTGCGAGCAGCCCtgaataa                                                                                                                                        <  1:3341724/150‑7 (MQ=255)

                    tttGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGt                                                                                                                                   >  2:247493/1‑150 (MQ=255)

                    tttGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGt                                                                                                                                   >  2:2012553/1‑150 (MQ=255)

                    tttGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGt                                                                                                                                   >  2:4003043/1‑150 (MQ=255)

                          cTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGt                                                                                                                             <  1:615202/150‑1 (MQ=255)

                           ttGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGtc                                                                                                                            <  2:698338/150‑1 (MQ=255)

                           ttGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGtc                                                                                                                            <  2:4301889/150‑1 (MQ=255)

                             ggTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGtctc                                                                                                                          <  2:173052/150‑1 (MQ=255)

                               tCGCGGCCCGACTGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGc                                                                                                                        <  1:943344/150‑1 (MQ=255)

                                     cccGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGt                                                                                                                  <  1:1743868/150‑1 (MQ=255)

                                      ccGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGtt                                                                                                                 >  1:2602881/1‑150 (MQ=255)

                                      ccGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCGGGTCATGGTCTCGCAGGCGtt                                                                                                                 >  2:4554519/1‑150 (MQ=255)

                                       cGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCa                                                                                                                                           <  1:2532225/123‑1 (MQ=255)

                                       cGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCa                                                                                                                                           >  2:2532225/1‑123 (MQ=255)

                                       cGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTc                                                                                                                <  2:4161412/150‑1 (MQ=255)

                                         aaTGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCaa                                                                                                              >  1:2235477/1‑150 (MQ=255)

                                            gTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTg                                                                                                           <  2:2529971/150‑1 (MQ=255)

                                             tCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTgg                                                                                                          <  2:1084433/150‑1 (MQ=255)

                                             tCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTgg                                                                                                          <  1:5209905/150‑1 (MQ=255)

                                              ccgccgCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGc                                                                                                         >  2:3941647/1‑150 (MQ=255)

                                               cgccgcCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCt                                                                                                        <  2:3028491/150‑1 (MQ=255)

                                                gccgccAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTg                                                                                                       >  1:4826967/1‑150 (MQ=255)

                                                   gccAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcg                                                                                                    >  2:2018914/1‑150 (MQ=255)

                                                   gccAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcg                                                                                                    >  1:1941502/1‑150 (MQ=255)

                                                      aGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGa                                                                                                 <  1:3688738/150‑1 (MQ=255)

                                                      aGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGa                                                                                                 <  2:3966634/150‑1 (MQ=255)

                                                      aGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCAGTCAAAGTCGAGCACAATGTGTTCGACCAGGGGGTCAGTGATGGCCTTATTCAGGGCAGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGTCGTTCAAGTGGCTGCCGCGa                                                                                                 >  2:3635925/1‑150 (MQ=255)

                                                          cAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGTCTGCTCGCAGCCCTTCATAGCTGTTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTca                                                                                             <  2:781901/150‑1 (MQ=255)

                                                           aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTcac                                                                                            >  1:3097330/1‑150 (MQ=255)

                                                           aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTcac                                                                                            >  2:4877205/1‑150 (MQ=255)

                                                           aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCGCGCAGGCGTTCAAGTGGCTGCCGCGACTcac                                                                                            >  1:4976234/1‑150 (MQ=255)

                                                             aaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACaa                                                                                          <  2:1511016/150‑1 (MQ=255)

                                                                  ccGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACg                                                                                     >  1:1369332/1‑150 (MQ=255)

                                                                            tGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGGTGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACgg                                                                           <  2:1753007/150‑1 (MQ=255)

                                                                            tGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACgg                                                                           <  2:4297457/150‑1 (MQ=255)

                                                                             gccgccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGaa                                                                                       >  1:3059133/1‑137 (MQ=255)

                                                                             gccgccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGaa                                                                                       <  2:3059133/137‑1 (MQ=255)

                                                                              ccgccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACgggg                                                                         <  2:2534951/150‑1 (MQ=255)

                                                                              ccgccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACgggg                                                                         <  2:2924158/150‑1 (MQ=255)

                                                                               caccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTAAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGa                                                                        <  2:2192789/148‑1 (MQ=255)

                                                                                gccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATg                                                                                                                               <  1:1424373/94‑1 (MQ=255)

                                                                                gccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATg                                                                                                                               >  2:1424373/1‑94 (MQ=255)

                                                                                gccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGAt                                                                       >  2:4620937/1‑150 (MQ=255)

                                                                                   tGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCAc                                                                    >  1:929119/1‑150 (MQ=255)

                                                                                   tGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCAc                                                                    <  1:1900159/150‑1 (MQ=255)

                                                                                   tGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCAc                                                                    <  1:3155608/150‑1 (MQ=255)

                                                                                    gggCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCAcc                                                                   <  2:4165708/150‑1 (MQ=255)

                                                                                     ggCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCt                                                                  >  2:5290469/1‑150 (MQ=255)

                                                                                      gCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTc                                                                 <  2:5244126/150‑1 (MQ=255)

                                                                                        tGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCca                                                               >  1:1625611/1‑150 (MQ=255)

                                                                                        tGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCca                                                               >  1:1625593/1‑150 (MQ=255)

                                                                                          tCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACACGAACGCCATGGACGGGGATCACCTCcaca                                                             >  2:2914691/1‑150 (MQ=255)

                                                                                           cAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCTCCTCCacac                                                            <  2:2602881/150‑1 (MQ=255)

                                                                                            aaTGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACAcc                                                           >  2:1132854/1‑150 (MQ=255)

                                                                                             aTGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACAccc                                                          <  2:1670949/150‑1 (MQ=255)

                                                                                                  gAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTg                                                     >  2:2971947/1‑150 (MQ=255)

                                                                                                  gAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTg                                                     >  1:4636561/1‑150 (MQ=255)

                                                                                                      acaATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCg                                                 >  2:2492760/1‑150 (MQ=255)

                                                                                                       caATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa                                                <  1:4361769/150‑1 (MQ=255)

                                                                                                       caATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa                                                >  2:1043312/1‑150 (MQ=255)

                                                                                                        aaTGTGTTCGACCATGGGGTCAGTTATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGAt                                               >  1:905723/1‑150 (MQ=255)

                                                                                                         aTGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATg                                              <  1:3355694/150‑1 (MQ=255)

                                                                                                           gtgtTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGa                                                                                        <  1:1518976/106‑1 (MQ=255)

                                                                                                           gtgtTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGa                                                                                        >  2:1518976/1‑106 (MQ=255)

                                                                                                              ttCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGATCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATggcggc                                         >  1:4292228/1‑150 (MQ=255)

                                                                                                                   ccATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGgcggcgcggc                                    <  1:1132854/150‑1 (MQ=255)

                                                                                                                     aTGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCggcgcggcgc                                  <  1:4620937/150‑1 (MQ=255)

                                                                                                                       ggggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTg                                <  1:2545911/150‑1 (MQ=255)

                                                                                                                        gggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcgc                                                                                                   <  2:3891017/82‑1 (MQ=255)

                                                                                                                        gggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcgc                                                                                                   >  1:3891017/1‑82 (MQ=255)

                                                                                                                        gggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGt                               >  1:1442102/1‑150 (MQ=255)

                                                                                                                        gggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGt                               >  1:4189182/1‑150 (MQ=255)

                                                                                                                           tCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGttct                            >  2:694627/1‑150 (MQ=255)

                                                                                                                             aGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGttcttc                          >  2:1184931/1‑150 (MQ=255)

                                                                                                                             aGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGttcttc                          >  2:5269879/1‑150 (MQ=255)

                                                                                                                             aGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGttcttc                          <  1:247493/150‑1 (MQ=255)

                                                                                                                              gTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcg                         >  2:5014376/1‑150 (MQ=255)

                                                                                                                               ggagggCCTTATGCAGGGCTGCTCGCAGCCCGGCATAGCTGGGCATGGTCTCGCTGGCGTTCACGTGGCTGCCGCGACCCCCCAGACCGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGCTGGCGGCGCGGCGCTGTTCTTcgg                        >  2:3409862/5‑149 (MQ=255)

                                                                                                                                gATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcgcg                       <  2:3738586/150‑1 (MQ=255)

                                                                                                                                     ccTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGAt                  >  1:2313152/1‑150 (MQ=255)

                                                                                                                                      cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATc                 >  1:3245720/1‑150 (MQ=255)

                                                                                                                                      cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATc                 >  1:2538587/1‑150 (MQ=255)

                                                                                                                                      cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATc                 <  2:1023071/150‑1 (MQ=255)

                                                                                                                                      cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATc                 <  1:1068843/150‑1 (MQ=255)

                                                                                                                                          ttCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGc             >  2:1515925/1‑150 (MQ=255)

                                                                                                                                            aaGGGCTGCGAGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTc           >  2:3895679/2‑150 (MQ=255)

                                                                                                                                            aaGGGCTGCGAGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTc           >  2:5426119/2‑150 (MQ=255)

                                                                                                                                             aGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCa          >  2:383707/1‑150 (MQ=255)

                                                                                                                                             aGGGCTGCGAGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCa          <  1:3895679/150‑1 (MQ=255)

                                                                                                                                             aGGGCTGCGAGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCa          >  1:523665/1‑150 (MQ=255)

                                                                                                                                               ggCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGc        <  2:93717/150‑1 (MQ=255)

                                                                                                                                                gctgctCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGcc       <  1:1516767/150‑1 (MQ=255)

                                                                                                                                                gctgctCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGcc       <  1:130891/150‑1 (MQ=255)

                                                                                                                                                     tCGCAGCCCTTCCGAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCCCCAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGCCTGGTa  <  2:4720120/150‑1 (MQ=255)

TTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTCAGCCTGGTA  >  NC_002947/4401383‑4401681

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