BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,193,818	A→G	6.1%	Q907R (CAG→CGG)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,193,818	0	A	G	6.1%	463.8 / 12.5	147	Q907R (CAG→CGG)	PP_1937	helicase
Reads supporting (aligned to +/- strand): ref base A (74/64); new base G (3/6); total (77/70)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 150 total aligned reads displayed.

ATGTGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCA  >  NC_002947/2193671‑2193967

aTGTGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGc                                                                                                                                                     <  1:1455800/150‑1 (MQ=255)

aTGTGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGc                                                                                                                                                     >  1:4627280/1‑150 (MQ=255)

     tATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGtc                                                                                                                                                >  1:3832303/1‑150 (MQ=255)

      aTGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGtca                                                                                                                                               <  1:3352845/150‑1 (MQ=255)

           cTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatcaga                                                                                                                                          <  1:3704430/150‑3 (MQ=255)

           cTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatcaga                                                                                                                                          <  1:825341/150‑3 (MQ=255)

              tATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAAtctc                                                                                                                                       >  2:601175/1‑150 (MQ=255)

               aTCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTcc                                                                                                                                      <  1:2505530/150‑1 (MQ=255)

                  gaaacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCcgcg                                                                                                                                   >  1:4912832/1‑150 (MQ=255)

                   aaacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCga                                                                                                                                  >  1:539636/1‑150 (MQ=255)

                         aacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCGCGGAGATTga                                                                                                                            <  1:1929678/150‑2 (MQ=255)

                         aacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGc                                                                                                                            <  2:3764273/150‑1 (MQ=255)

                          acTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTgat                                                                                                                           <  2:4447311/150‑3 (MQ=255)

                              aCAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGc                                                                                                                       >  1:1857386/1‑150 (MQ=255)

                               cAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCa                                                                                                                      <  2:3809221/150‑1 (MQ=255)

                                aaGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTTGATCTTCCTGTCCGTCTCTGATGACAAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAg                                                                                                                     >  1:5135950/1‑150 (MQ=255)

                                aaGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAg                                                                                                                     >  1:2111468/1‑150 (MQ=255)

                                    tCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAAGGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGa                                                                                                                 <  1:2903864/150‑1 (MQ=255)

                                     cATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAg                                                                                                                <  2:2082027/150‑1 (MQ=255)

                                     cATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAg                                                                                                                <  2:1775522/150‑1 (MQ=255)

                                          cgcgTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCAc                                                                                                           >  1:1171379/1‑150 (MQ=255)

                                          cgcgTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCAc                                                                                                           >  1:1171440/1‑150 (MQ=255)

                                            cgTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCGTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACgg                                                                                                         >  2:3776541/1‑150 (MQ=255)

                                             gtgtGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCCTCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGt                                                                                                        >  1:5158350/1‑150 (MQ=255)

                                              tgtgAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGtt                                                                                                       <  1:3787995/150‑1 (MQ=255)

                                              tgtgAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGtt                                                                                                       >  2:3787995/1‑150 (MQ=255)

                                                tgAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTAt                                                                                                     <  2:966034/150‑1 (MQ=255)

                                                 gAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCCCGGTGATc                                                                                                    >  1:4700661/1‑150 (MQ=255)

                                                 gAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATc                                                                                                    <  2:4940719/150‑1 (MQ=255)

                                                 gAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCGGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATc                                                                                                    <  2:3985379/150‑1 (MQ=255)

                                                  aGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCa                                                                                                   <  1:3183688/150‑1 (MQ=255)

                                                     gTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCga                                                                                                                                  >  1:2730782/1‑116 (MQ=255)

                                                     gTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCga                                                                                                                                  <  2:2730782/116‑1 (MQ=255)

                                                           ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                          >  2:1213783/1‑150 (MQ=255)

                                                           ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                          <  2:3643921/150‑1 (MQ=255)

                                                           ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                          >  1:2783093/1‑150 (MQ=255)

                                                           ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                          >  1:2783159/1‑150 (MQ=255)

                                                            cTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGc                                                                                         >  2:2602308/1‑150 (MQ=255)

                                                             ttCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCa                                                                                        >  1:3419937/1‑150 (MQ=255)

                                                                cAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGaa                                                                                     <  2:944634/150‑1 (MQ=255)

                                                                 aGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAAc                                                                                    <  2:235781/150‑1 (MQ=255)

                                                                      aCCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCaa                                                                               >  2:627459/1‑150 (MQ=255)

                                                                           gggTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGa                                                                          <  2:2907696/150‑1 (MQ=255)

                                                                               gtTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAg                                                                      <  1:4047908/150‑1 (MQ=255)

                                                                               gtTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAg                                                                      <  1:3673629/150‑1 (MQ=255)

                                                                                tttCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAgg                                                                     >  1:3083242/1‑150 (MQ=255)

                                                                                 ttCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAggg                                                                    >  1:3713001/1‑150 (MQ=255)

                                                                                        cTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAAt                                                             >  1:1849937/1‑150 (MQ=255)

                                                                                        cTTGGCCGGTGAGAAGGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAAt                                                             <  1:5014750/150‑1 (MQ=255)

                                                                                          tGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATgg                                                           <  1:5036556/150‑1 (MQ=255)

                                                                                              cGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAt                                                       >  1:551244/1‑150 (MQ=255)

                                                                                               ggTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCACCCCGACAAGGGCTGCAATGGGCAtg                                                      >  2:2294777/1‑150 (MQ=255)

                                                                                               ggTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtg                                                      <  2:3318972/150‑1 (MQ=255)

                                                                                                 tGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtgct                                                    <  1:4377247/150‑1 (MQ=255)

                                                                                                 tGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtgct                                                    >  1:3792240/1‑150 (MQ=255)

                                                                                                 tGAGAAGGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtgct                                                    <  1:1415588/150‑1 (MQ=255)

                                                                                                      aTGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAAt                                               <  2:3183648/150‑1 (MQ=255)

                                                                                                      aTGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAAt                                               >  1:1302895/1‑150 (MQ=255)

                                                                                                      aTGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAAt                                               <  1:1729693/150‑1 (MQ=255)

                                                                                                       tGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGa                                                                                                                 <  1:1057557/83‑1 (MQ=255)

                                                                                                        gcgcgcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCg                                             >  1:3460690/1‑150 (MQ=255)

                                                                                                            gcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTc                                         <  2:672300/150‑1 (MQ=255)

                                                                                                            gcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTc                                         <  2:539636/150‑1 (MQ=255)

                                                                                                             cgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCt                                        <  1:2428528/150‑1 (MQ=255)

                                                                                                             cgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGCCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCt                                        <  2:4330785/150‑1 (MQ=255)

                                                                                                              gcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTg                                       <  1:4699203/150‑1 (MQ=255)

                                                                                                              gcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTg                                       <  1:4742821/150‑1 (MQ=255)

                                                                                                              gcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTg                                       <  2:3665713/150‑1 (MQ=255)

                                                                                                                aCCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGc                                     >  1:1576625/1‑150 (MQ=255)

                                                                                                                aCCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGc                                     >  1:1576691/1‑150 (MQ=255)

                                                                                                                aCCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGc                                     <  1:4863957/150‑1 (MQ=255)

                                                                                                                aCCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGc                                     <  1:2705359/150‑1 (MQ=255)

                                                                                                                   aaGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAg                                  >  2:3979118/1‑150 (MQ=255)

                                                                                                                   aaGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAg                                  >  1:968384/1‑150 (MQ=255)

                                                                                                                        ggATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTgg                             <  2:5317383/150‑1 (MQ=255)

                                                                                                                        ggATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTgg                             <  2:4808829/150‑1 (MQ=255)

                                                                                                                              tCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCgaga                       <  2:691191/150‑1 (MQ=255)

                                                                                                                              tCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCgaga                       <  2:1662014/150‑1 (MQ=255)

                                                                                                                                cTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCTCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATg                     <  1:1816644/150‑1 (MQ=255)

                                                                                                                                cTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATg                     <  2:1576691/150‑1 (MQ=255)

                                                                                                                                cTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATg                     <  2:1576625/150‑1 (MQ=255)

                                                                                                                                  gTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGga                   >  2:3294832/1‑150 (MQ=255)

                                                                                                                                     ggagattgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaa                >  1:4671046/7‑150 (MQ=255)

                                                                                                                                     cGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaa                >  1:2658210/1‑150 (MQ=255)

                                                                                                                                      ggCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaa               <  1:4271724/150‑1 (MQ=255)

                                                                                                                                      ggCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaa               >  1:2861023/1‑150 (MQ=255)

                                                                                                                                      ggCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaa               >  1:1074012/1‑150 (MQ=255)

                                                                                                                                      ggCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaa               >  2:5346550/1‑150 (MQ=255)

                                                                                                                                       gCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGTCGATTTGCAGAACCTCCACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTTGGCGAGATGgagtata              <  2:4002247/150‑5 (MQ=255)

                                                                                                                                        gattgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAAt             >  2:4604893/4‑150 (MQ=255)

                                                                                                                                          ctGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGc           <  2:4912832/150‑1 (MQ=255)

                                                                                                                                           tgatgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCt          <  2:1345959/150‑1 (MQ=255)

                                                                                                                                             atgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGgaga                 <  1:1984800/141‑1 (MQ=255)

                                                                                                                                             atgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGgaga                 >  2:1984800/1‑141 (MQ=255)

                                                                                                                                             atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTAc        >  2:584699/1‑150 (MQ=255)

                                                                                                                                             atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTAc        >  1:136683/1‑150 (MQ=255)

                                                                                                                                                aCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGc     <  1:601147/150‑1 (MQ=255)

                                                                                                                                                aCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGc     <  1:601175/150‑1 (MQ=255)

                                                                                                                                                 ccAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGcc    >  1:3639105/1‑150 (MQ=255)

                                                                                                                                                   aGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTGTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCa  >  2:5119736/1‑150 (MQ=255)

ATGTGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCA  >  NC_002947/2193671‑2193967

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