BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,602,918	G→C	5.0%	S94T (AGT→ACT)	garK →	glycerate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,602,918	0	G	C	5.0%	364.8 / 7.5	112	S94T (AGT→ACT)	garK	glycerate kinase
Reads supporting (aligned to +/- strand): ref base G (54/52); new base C (3/3); total (57/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GCGAAGGCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGAT  >  NC_002947/3602771‑3603057

gCGAAGGCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTg                                                                                                                                           >  1:974424/1‑150 (MQ=255)

  gAAGGCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGc                                                                                                                                         >  2:443848/1‑150 (MQ=255)

  gAAGGCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGc                                                                                                                                         <  2:1117253/150‑1 (MQ=255)

      gCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGc                                                                                                                                     >  2:850765/1‑150 (MQ=255)

       cACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCACTGGCCTggg                                                                                                                                    <  2:737876/150‑3 (MQ=255)

        aCCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCaa                                                                                                                                   <  1:850765/150‑1 (MQ=255)

        aCCATGGAGGTGATCCTGGCCGCCAGACACGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCaa                                                                                                                                   >  1:1142024/1‑150 (MQ=255)

          cATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACa                                                                                                                                 >  2:140242/1‑150 (MQ=255)

          cATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACa                                                                                                                                 >  2:1729235/1‑150 (MQ=255)

          cATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACa                                                                                                                                 >  1:1625754/1‑150 (MQ=255)

               aGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGc                                                                                                                            >  2:612820/1‑150 (MQ=255)

                ggTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAgg                                                                                                                               >  2:83284/1‑146 (MQ=255)

                ggTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAgg                                                                                                                               <  1:83284/146‑1 (MQ=255)

                ggTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGcc                                                                                                                           >  1:15988/1‑150 (MQ=255)

                 gtcgtCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCa                                                                                                                          >  2:1836887/5‑150 (MQ=255)

                   gATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAAc                                                                                                                        >  2:672494/1‑150 (MQ=255)

                   gATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAAc                                                                                                                        >  1:327654/1‑150 (MQ=255)

                    aTCCTGGCCGCCAGCCGCAGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACa                                                                                                                       <  1:391020/150‑1 (MQ=255)

                     tCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCTACAGGTGCCAACAg                                                                                                                      <  2:1142037/150‑1 (MQ=255)

                     tCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAg                                                                                                                      <  2:10496/150‑1 (MQ=255)

                     tCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAg                                                                                                                      >  1:1827780/1‑150 (MQ=255)

                       cTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGc                                                                                                                    <  1:405938/150‑1 (MQ=255)

                        tGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGcc                                                                                                                   >  1:1798458/1‑150 (MQ=255)

                                 gCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATg                                                                                                          <  1:782560/150‑1 (MQ=255)

                                  ccGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGc                                                                                                         <  2:1580220/150‑1 (MQ=255)

                                    gcgGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCtt                                                                                                       <  2:465350/150‑1 (MQ=255)

                                        tGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGCTGCTTGCCg                                                                                                   >  2:216445/1‑150 (MQ=255)

                                         gAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCgc                                                                                                  <  1:673587/150‑1 (MQ=255)

                                           aTTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCgcag                                                                                                <  2:1798458/150‑1 (MQ=255)

                                            ttGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCgcagc                                                                                               >  1:75806/1‑150 (MQ=255)

                                                    ccAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTg                                                                                       <  2:522436/150‑1 (MQ=255)

                                                    ccAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTg                                                                                       >  1:522436/1‑150 (MQ=255)

                                                    ccAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGCAGCCCCTCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTg                                                                                       <  1:612820/150‑1 (MQ=255)

                                                     cAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGGTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGc                                                                                                                                         >  1:1571941/1‑99 (MQ=255)

                                                     cAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGc                                                                                                                                         <  2:1571941/99‑1 (MQ=255)

                                                      aGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTggg                                                                                     <  2:1840343/150‑1 (MQ=255)

                                                      aGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCACACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTggg                                                                                     <  1:874953/150‑1 (MQ=255)

                                                       gACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGTTGCCAACAGGCCAACGCGATGCTTGCCgcagcagcatcgcgtt                                                                                    <  2:353914/150‑8 (MQ=255)

                                                              cgcgGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGt                                                                             <  1:1728485/150‑1 (MQ=255)

                                                                 gggCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGaa                                                                          >  2:395251/1‑150 (MQ=255)

                                                                 gggCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGaa                                                                          <  1:1669573/150‑1 (MQ=255)

                                                                   gCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACt                                                                        <  2:1048067/150‑1 (MQ=255)

                                                                   gCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACt                                                                        >  1:1482366/1‑150 (MQ=255)

                                                                     cACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGa                                                                      >  1:1473118/1‑150 (MQ=255)

                                                                     cACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAAAACTCGAAATGGCCCAGGAACCTGGTGGTAAGCTACGACCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGa                                                                      >  1:1407458/1‑150 (MQ=255)

                                                                      aCTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTgg                                                                                      <  2:1730112/133‑1 (MQ=255)

                                                                      aCTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTgg                                                                                      >  1:1730112/1‑133 (MQ=255)

                                                                      aCTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGAt                                                                     <  1:751732/150‑1 (MQ=255)

                                                                      aCTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGAt                                                                     >  1:605318/1‑150 (MQ=255)

                                                                       cTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTCGCCGCAGCAGCACCTGGGGTACCGGTGAACTGAtt                                                                    >  1:1727797/1‑150 (MQ=255)

                                                                        tGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGTAAAGTCGCACCGCAATCATCGAAATGGCCCATGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTg                                                                   <  1:1789867/150‑1 (MQ=255)

                                                                          gggCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGcc                                                                 <  2:1136289/150‑1 (MQ=255)

                                                                              aGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAg                                                             <  2:1155127/150‑1 (MQ=255)

                                                                                   gTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCtctc                                                        >  2:1433397/1‑150 (MQ=255)

                                                                                   gTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGACGCGGCCAGTGGCCTGCAGCAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCtctc                                                        >  1:1691512/1‑150 (MQ=255)

                                                                                    tCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCg                                                       >  1:1681317/1‑150 (MQ=255)

                                                                                      gAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGAc                                                     >  1:690380/1‑150 (MQ=255)

                                                                                          atggctggGGCTGGTTGCCGCAAAGTCGCACCGTAATAATCGTAATGGCCCAGGACAGTGGCATGCAACATATGTCAGTAGCTCAACGAGATGATTGAAGCAGCAGCTCCTGGGGTAACGGTGAACTAATTGCCGCAGCTCTTGACGCCg                                                 <  2:949745/149‑1 (MQ=255)

                                                                                           tggctggGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTAc                                                                                 >  1:1056775/1‑117 (MQ=255)

                                                                                           tggctggGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTAc                                                                                 <  2:1056775/117‑1 (MQ=255)

                                                                                           tggctggGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCgg                                                >  1:257539/1‑150 (MQ=255)

                                                                                                ggggctggTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGccc                                           >  1:1669817/1‑150 (MQ=255)

                                                                                                       gTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCgtgt                                    <  1:140242/150‑1 (MQ=255)

                                                                                                       gTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCgtgt                                    >  2:1488146/1‑150 (MQ=255)

                                                                                                        ttGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGa                                                      <  1:173098/131‑1 (MQ=255)

                                                                                                        ttGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGa                                                      >  2:173098/1‑131 (MQ=255)

                                                                                                        ttGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGtgtg                                   <  1:811364/150‑1 (MQ=255)

                                                                                                           ccGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTgtggtg                                >  2:703692/1‑150 (MQ=255)

                                                                                                            cGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGc                               >  1:1713275/1‑150 (MQ=255)

                                                                                                             gCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCt                              >  2:1659624/1‑150 (MQ=255)

                                                                                                             gCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTACACACCTGCGCCTGGCGTGTGGTGCt                              <  2:1142024/150‑1 (MQ=255)

                                                                                                              cAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTg                             <  2:1606699/150‑1 (MQ=255)

                                                                                                               aaaGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTgg                            >  1:1236089/1‑150 (MQ=255)

                                                                                                                     gCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcg                      <  1:1433397/150‑1 (MQ=255)

                                                                                                                     gCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGAAGTTGGCGACGCCGGCGCCCGGCGTGTGGTGCTGACCATcg                      >  1:1558278/1‑150 (MQ=255)

                                                                                                                       aCCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGACATGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcggc                    <  1:443848/150‑1 (MQ=255)

                                                                                                                       aCCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcggc                    >  2:1068185/1‑150 (MQ=255)

                                                                                                                       aCCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcggc                    <  2:1407458/150‑1 (MQ=255)

                                                                                                                        ccGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGa                                                      >  2:117091/1‑115 (MQ=255)

                                                                                                                        ccGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGa                                                      <  1:117091/115‑1 (MQ=255)

                                                                                                                         cGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGc                           <  2:1719331/141‑1 (MQ=255)

                                                                                                                         cGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGc                           >  1:1719331/1‑141 (MQ=255)

                                                                                                                         cGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcggcgg                  <  1:1729235/150‑1 (MQ=255)

                                                                                                                         cGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATcggcgg                  >  1:1915765/1‑150 (MQ=255)

                                                                                                                           cAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCa                <  2:1669817/150‑1 (MQ=255)

                                                                                                                            aaTCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAg               >  1:1562552/1‑150 (MQ=255)

                                                                                                                               catcGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAgcgc            <  1:196822/150‑1 (MQ=255)

                                                                                                                                 tcGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGcca          >  2:485805/1‑150 (MQ=255)

                                                                                                                                  cGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGccac         >  1:484808/1‑150 (MQ=255)

                                                                                                                                       tGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAAc     >  2:726621/1‑149 (MQ=255)

                                                                                                                                       tGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAAc     <  1:726621/149‑1 (MQ=255)

                                                                                                                                        gttgcAGGCCACTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAAc                                                                         <  2:216391/76‑1 (MQ=255)

                                                                                                                                        gttgcAGGCCACTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAAc                                                                         >  1:216391/5‑80 (MQ=255)

                                                                                                                                         gCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACATGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGAt  <  2:1625754/150‑1 (MQ=255)

                                                                                                                                           ccAGGCCAGTGACCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCtctc                                                        <  1:95805/94‑1 (MQ=255)

                                                                                                                                           ccAGGCCAGTGACCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCtctc                                                        >  2:95805/1‑94 (MQ=255)

                                                                                                                                               gCCACTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGc                                                                                                     <  2:576362/45‑1 (MQ=255)

                                                                                                                                               gCCACTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGc                                                                                                     >  1:576362/1‑45 (MQ=255)

                                                                                                                                                 cAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGAtt                                                                    >  2:313851/1‑76 (MQ=255)

                                                                                                                                                 cAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGAtt                                                                    <  1:313851/76‑1 (MQ=255)

GCGAAGGCACCATGGAGGTGATCCTGGCCGCCAGCCGCGGTGAATTGCGCCGCCAGACCGTGCGCGGGCCACTGGGGCAGTACGTCGAGGCTGGCTGGGGCTGGTTGCCGCAAAGTCGCACCGCAATCATCGAAATGGCCCAGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTCGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGAT  >  NC_002947/3602771‑3603057

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