BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,750,775	C→G	7.4%	A31G (GCC→GGC)	PP_1551 →	replication protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,750,775	0	C	G	7.4%	672.1 / 32.6	201	A31G (GCC→GGC)	PP_1551	replication protein O
Reads supporting (aligned to +/- strand): ref base C (108/77); new base G (11/4); total (120/81)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCC  >  NC_002947/1750628‑1750924

gAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggccc                                                                                                                                                     >  2:1750076/1‑150 (MQ=255)

gAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggccc                                                                                                                                                     >  2:3193966/1‑150 (MQ=255)

 aaaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccg                                                                                                                                                    >  1:103041/1‑150 (MQ=255)

   aaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCggg                                                                                                                                                  >  1:3268213/1‑150 (MQ=255)

   aaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCggg                                                                                                                                                  >  2:4284969/1‑150 (MQ=255)

   aaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCggg                                                                                                                                                  >  1:4437988/1‑150 (MQ=255)

    aaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGga                                                                                                                                                 <  2:3999017/150‑1 (MQ=255)

    aaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGga                                                                                                                                                 <  2:1356359/150‑1 (MQ=255)

     aaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgag                                                                                                                                                <  1:348403/150‑1 (MQ=255)

     aaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgag                                                                                                                                                >  2:287250/1‑150 (MQ=255)

      aaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAACGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgaga                                                                                                                                               >  2:2696110/1‑150 (MQ=255)

      aaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgaga                                                                                                                                               >  2:3866523/1‑150 (MQ=255)

      aaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgaga                                                                                                                                               <  2:1277853/150‑1 (MQ=255)

      aaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgaga                                                                                                                                               <  2:3538331/150‑1 (MQ=255)

          cGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGaa                                                                                                                                           <  1:2543372/150‑1 (MQ=255)

           gCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAg                                                                                                                                          <  1:334969/150‑1 (MQ=255)

             cTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAt                                                                                                                                        >  2:708856/1‑150 (MQ=255)

              tGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                       >  1:2651015/1‑150 (MQ=255)

                 cAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                    <  1:2162843/150‑1 (MQ=255)

                 cAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                    <  2:1285279/150‑1 (MQ=255)

                 cAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                    <  1:1164962/150‑1 (MQ=255)

                  aGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGc                                                                                                                                   >  1:2514551/1‑150 (MQ=255)

                         cTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGt                                                                                                                            <  1:1527685/150‑1 (MQ=255)

                         cTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGt                                                                                                                            <  2:2386231/150‑1 (MQ=255)

                         cTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGt                                                                                                                            <  2:1073110/150‑1 (MQ=255)

                          tttCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                           >  1:2987351/1‑150 (MQ=255)

                              tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCa                                                                                                                       >  1:1428286/1‑150 (MQ=255)

                                tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                     <  1:50083/150‑1 (MQ=255)

                                tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                     >  1:4169137/1‑150 (MQ=255)

                                 aCAGCATTAGAACGGGTTAAAGCATGACAAACATCGACCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAgg                                                                                                                    >  2:2884222/1‑150 (MQ=255)

                                  cAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGc                                                                                                                   >  1:887099/1‑150 (MQ=255)

                                      aTTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacg                                                                                                               <  1:2118082/150‑2 (MQ=255)

                                      aTTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGccacc                                                                                                               >  2:1525398/1‑150 (MQ=255)

                                           caCGGGTTACAGCATGACCAACCTCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCaa                                                                                                          <  1:4051947/149‑1 (MQ=255)

                                           aaCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCaa                                                                                                          >  1:3235475/1‑150 (MQ=255)

                                                gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                     <  1:3649038/150‑1 (MQ=255)

                                                gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                     <  2:3312817/150‑1 (MQ=255)

                                                gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                     <  2:1248639/150‑1 (MQ=255)

                                                gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                     >  2:1036365/1‑150 (MQ=255)

                                                 ttAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCgg                                                                                                    <  1:3066204/150‑1 (MQ=255)

                                                  tAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGt                                                                                                   <  2:1689578/150‑1 (MQ=255)

                                                   aaaGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTg                                                                                                  >  2:4460594/1‑150 (MQ=255)

                                                      gCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCgg                                                                                               <  2:370410/150‑1 (MQ=255)

                                                      gCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCACCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCgg                                                                                               <  2:369304/150‑1 (MQ=255)

                                                           aCAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCt                                                                                          >  2:2135113/1‑150 (MQ=255)

                                                           aCAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCt                                                                                          >  1:1211838/1‑150 (MQ=255)

                                                            cAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTc                                                                                         >  1:1074869/1‑150 (MQ=255)

                                                             aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                           <  2:2873737/115‑1 (MQ=255)

                                                             aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                           >  1:2873737/1‑115 (MQ=255)

                                                             aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCa                                                                                        <  2:3799703/150‑1 (MQ=255)

                                                              aaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCaa                                                                                       <  2:2332686/150‑1 (MQ=255)

                                                               aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                      <  2:734000/150‑1 (MQ=255)

                                                                      cccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATccc                                                                               <  2:3437055/150‑1 (MQ=255)

                                                                          cTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCt                                                                           >  1:190119/1‑150 (MQ=255)

                                                                                aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                                     >  1:3345491/1‑150 (MQ=255)

                                                                                aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                                     >  1:4137166/1‑150 (MQ=255)

                                                                                aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                                     >  2:1648682/1‑150 (MQ=255)

                                                                                     cAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                                <  2:532828/150‑1 (MQ=255)

                                                                                     cAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                                <  2:4159332/150‑1 (MQ=255)

                                                                                        ggggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaa                                                             <  1:737295/150‑1 (MQ=255)

                                                                                         gggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                            <  2:2599556/150‑1 (MQ=255)

                                                                                           gTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGc                                                          <  1:3172980/150‑1 (MQ=255)

                                                                                                cccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                     >  1:3407019/1‑86 (MQ=255)

                                                                                                cccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATcc                                                     >  2:1828005/1‑150 (MQ=255)

                                                                                                 ccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCa                                                    >  1:858681/1‑150 (MQ=255)

                                                                                                    gATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGc                                                 >  1:600186/1‑150 (MQ=255)

                                                                                                       ggACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCa                                              <  2:2987351/150‑1 (MQ=255)

                                                                                                         aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                            >  2:3933636/1‑150 (MQ=255)

                                                                                                          cAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAccc                                           <  1:3975239/150‑1 (MQ=255)

                                                                                                           aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGACAGCCATCCACGCTCACCCt                                          <  2:1002913/150‑1 (MQ=255)

                                                                                                           aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAACGCCATCCACGCTCACCCt                                          >  2:3023656/1‑150 (MQ=255)

                                                                                                                  tCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCg                                   >  2:1937999/1‑150 (MQ=255)

                                                                                                                  tCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCg                                   >  2:4206541/1‑150 (MQ=255)

                                                                                                                    atggatggCCTGCTGGCTATAGATCTCCCGGACCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                                 >  1:354363/1‑150 (MQ=255)

                                                                                                                      ggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTc                               >  2:2666864/1‑150 (MQ=255)

                                                                                                                      ggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTc                               <  2:2730979/150‑1 (MQ=255)

                                                                                                                            ccTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                         >  1:3830139/1‑150 (MQ=255)

                                                                                                                             ctgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCg                        >  2:2186532/1‑150 (MQ=255)

                                                                                                                              tgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCgg                       >  2:204910/1‑150 (MQ=255)

                                                                                                                               gctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGt                      <  2:3288195/150‑1 (MQ=255)

                                                                                                                               gctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGt                      >  1:2517960/1‑150 (MQ=255)

                                                                                                                                ctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGtt                     <  1:1857005/150‑1 (MQ=255)

                                                                                                                                  ggCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                   >  2:4117843/1‑150 (MQ=255)

                                                                                                                                    cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCa                 <  1:512047/150‑1 (MQ=255)

                                                                                                                                     tATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAg                >  2:813211/1‑150 (MQ=255)

                                                                                                                                      aTCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGc               <  1:4460594/150‑1 (MQ=255)

                                                                                                                                       ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc              >  1:1931291/4‑150 (MQ=255)

                                                                                                                                       ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc              >  1:1254230/4‑150 (MQ=255)

                                                                                                                                       ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc              >  1:1254285/4‑150 (MQ=255)

                                                                                                                                        tcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCct             >  2:2247036/3‑150 (MQ=255)

                                                                                                                                        tcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCct             >  1:490378/3‑150 (MQ=255)

                                                                                                                                          aTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgc           >  2:781817/1‑150 (MQ=255)

                                                                                                                                           tctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgct          >  2:3061259/1‑150 (MQ=255)

                                                                                                                                             tcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCGTCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCtgctgc        >  2:833180/1‑150 (MQ=255)

                                                                                                                                             tcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCtgctgc        >  2:4086521/1‑150 (MQ=255)

                                                                                                                                               ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGt      >  2:4252575/1‑150 (MQ=255)

                                                                                                                                                cggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTc     <  1:1330729/150‑1 (MQ=255)

                                                                                                                                                cggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTc     <  1:1461070/150‑1 (MQ=255)

                                                                                                                                                cggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTc     <  1:1750076/150‑1 (MQ=255)

                                                                                                                                                  gcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGc   <  1:287250/150‑1 (MQ=255)

GAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCC  >  NC_002947/1750628‑1750924

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