BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,248	G→A	6.0%	P202P (CCG→CCA)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,248	0	G	A	6.0%	413.8 / 10.9	134	P202P (CCG→CCA)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base G (56/70); new base A (5/3); total (61/73)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 137 total aligned reads displayed.

AACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCC  >  NC_002947/1550099‑1550395

 aCTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCg                                                                                                                                                     >  1:4490689/1‑149 (MQ=255)

 aCTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCg                                                                                                                                                     <  2:4490689/149‑1 (MQ=255)

  cTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTa                                                                                                                                                   <  1:3360372/150‑1 (MQ=255)

   tCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTAc                                                                                                                                                  <  1:4196342/150‑1 (MQ=255)

      gACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTc                                                                                                                                               >  1:839711/1‑150 (MQ=255)

      gACAACTCCAGCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTc                                                                                                                                               <  2:226155/150‑1 (MQ=255)

         aaCTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTc                                                                                                                                            <  1:4924865/150‑1 (MQ=255)

           cTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaa                                                                                                                                          <  2:747816/150‑1 (MQ=255)

            tCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaac                                                                                                                                         <  1:4610996/150‑1 (MQ=255)

            tCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaac                                                                                                                                         <  1:20025/150‑1 (MQ=255)

             ccATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaaca                                                                                                                                        >  2:586285/1‑150 (MQ=255)

             ccATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaaca                                                                                                                                        <  1:4681393/150‑1 (MQ=255)

                    tGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGa                                                                                                                                 <  1:2649368/150‑1 (MQ=255)

                    tGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGa                                                                                                                                 <  2:1483914/150‑1 (MQ=255)

                        atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                             <  1:4638657/150‑1 (MQ=255)

                        atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                             <  1:4638923/150‑1 (MQ=255)

                        atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                             >  1:4257527/1‑150 (MQ=255)

                         tgatcGCTGAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGacaaga                                                                                                                            <  2:2804242/148‑4 (MQ=255)

                           atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                          >  2:3248376/1‑150 (MQ=255)

                           atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                          <  1:4957426/150‑1 (MQ=255)

                            tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                         >  1:2598555/1‑150 (MQ=255)

                            tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                         >  2:1581516/1‑150 (MQ=255)

                            tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                         >  1:4838727/1‑150 (MQ=255)

                            tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                         >  1:1405263/1‑150 (MQ=255)

                              gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                       <  1:708264/147‑1 (MQ=255)

                              gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                       >  2:1640083/4‑150 (MQ=255)

                              gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                       >  2:1640099/4‑150 (MQ=255)

                               ccGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTg                                                                                                                      >  2:4870210/1‑150 (MQ=255)

                               ccGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTg                                                                                                                      >  2:3878230/1‑150 (MQ=255)

                                tgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGTTGc                                                                                                                     <  2:1206031/149‑1 (MQ=255)

                                 gAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGcc                                                                                                                    <  2:357938/150‑1 (MQ=255)

                                  aaGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCg                                                                                                                   >  1:1799413/1‑150 (MQ=255)

                                     ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                                                                                >  1:4543610/1‑150 (MQ=255)

                                     ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                                                                                >  1:4543638/1‑150 (MQ=255)

                                     ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                                                                                <  2:2207475/150‑1 (MQ=255)

                                      cATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCt                                                                                                               >  1:2072622/1‑150 (MQ=255)

                                          gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                                                           >  2:372604/1‑150 (MQ=255)

                                          gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                                                           <  2:3365230/150‑1 (MQ=255)

                                           aaaaaGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTCTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGCTGCCGAGCTGGAAg                                                                                                          <  1:198000/150‑1 (MQ=255)

                                             aaaGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTGCCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGc                                                                                                        >  1:3266636/1‑150 (MQ=255)

                                               aGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGccc                                                                                                      <  2:2434435/150‑1 (MQ=255)

                                               aGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGccc                                                                                                      <  2:2357134/150‑1 (MQ=255)

                                                 tCGGTAAAGAAGGCGTGATTACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                                                                    >  1:5029348/1‑150 (MQ=255)

                                                 tCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                                                                    <  2:2719932/150‑1 (MQ=255)

                                                    gTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgctgc                                                                                                 >  2:5070327/1‑150 (MQ=255)

                                                       aagaagGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgctgctgc                                                                                              <  1:1042920/150‑1 (MQ=255)

                                                            ggCGAGATCCCCGTTGAAGAAGGCTCGGGCCGGGAAAACGAACTGTCCGACGTAGAAGGCAAGCAGTTCGACCGCGGCAACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGCTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                                                                         >  2:412156/1‑150 (MQ=255)

                                                                  aTCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGGACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                   <  2:411710/150‑1 (MQ=255)

                                                                    cACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACaa                                                                                 >  1:4585825/1‑150 (MQ=255)

                                                                     aCCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAg                                                                                >  1:2734205/1‑150 (MQ=255)

                                                                     aCCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAg                                                                                >  1:2734522/1‑150 (MQ=255)

                                                                      ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                                                                               <  2:2061347/150‑1 (MQ=255)

                                                                      ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                                                                               >  1:1020613/1‑150 (MQ=255)

                                                                        gTTGAAGTAGGCTAGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGTCATGCAGTTCGACCGCGGCTACCTGTCGACGTACTTCGTCATCCAGCCGGACTCCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCAGCTGGTCGACAAGaaa                                                                             <  2:841994/150‑1 (MQ=255)

                                                                        gTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGaaa                                                                             >  1:509875/1‑150 (MQ=255)

                                                                            aagaagGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACACGAAAAtct                                                                         >  2:1212196/1‑150 (MQ=255)

                                                                                  gCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACa                                                                   >  2:1517862/1‑150 (MQ=255)

                                                                                   cTCGGGTCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACAt                                                                  <  2:670751/150‑1 (MQ=255)

                                                                                    tCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTATCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATc                                                                 >  1:2278669/1‑150 (MQ=255)

                                                                                     cGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATcc                                                                <  2:1514774/150‑1 (MQ=255)

                                                                                       ggCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGt                                                              <  1:2593070/150‑1 (MQ=255)

                                                                                         ccTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGa                                                            <  2:3435918/150‑1 (MQ=255)

                                                                                           tGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgc                                                          <  2:3025222/150‑1 (MQ=255)

                                                                                           acgAAAACCAACTCTCCGTCCTAGAACGCAACCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCCGACACCATGCTAGCCGAGCTGGAAGGCCCGCTGCTGCTGCACGTCGACAACAAAATCTCCAACATCCGTGAgc                                                          <  2:3490645/148‑1 (MQ=255)

                                                                                               aaaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgctgct                                                      <  1:837723/150‑1 (MQ=255)

                                                                                                aacgaacgGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGctgctg                                                     <  1:3787619/150‑1 (MQ=255)

                                                                                                 aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAg                                                                                                                 >  1:1376258/1‑89 (MQ=255)

                                                                                                 aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAg                                                                                                                 <  2:1376258/89‑1 (MQ=255)

                                                                                                 aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCtgctgc                                                    <  1:3730789/150‑1 (MQ=255)

                                                                                                 aCGAACTGGCAGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGGTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCACCATCCGTGAGCtgctgc                                                    <  1:3730760/150‑1 (MQ=255)

                                                                                                     aCTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                                                                         >  2:4974126/1‑109 (MQ=255)

                                                                                                       tGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTc                                              >  1:817423/1‑150 (MQ=255)

                                                                                                        gTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCt                                             <  1:3978714/150‑1 (MQ=255)

                                                                                                        gTCCGTCGTAGAAGGCATGAAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCt                                             <  1:524131/150‑1 (MQ=255)

                                                                                                           cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                   <  2:2213476/109‑1 (MQ=255)

                                                                                                           cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                   >  1:2213476/1‑109 (MQ=255)

                                                                                                           cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGa                                          >  1:153331/1‑150 (MQ=255)

                                                                                                            gtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGaa                                         <  1:472586/150‑1 (MQ=255)

                                                                                                             tcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAg                                        >  2:884232/1‑150 (MQ=255)

                                                                                                              cgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGc                                       <  2:3591730/150‑1 (MQ=255)

                                                                                                              cgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGc                                       <  2:2956017/150‑1 (MQ=255)

                                                                                                               gtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGcc                                      <  1:3525968/150‑1 (MQ=255)

                                                                                                                tAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCg                                     >  2:1991931/1‑150 (MQ=255)

                                                                                                                tAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCg                                     >  2:5172331/1‑150 (MQ=255)

                                                                                                                 aGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGt                                    >  2:4628122/1‑150 (MQ=255)

                                                                                                                  gAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCAGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGtt                                   >  1:4239243/1‑150 (MQ=255)

                                                                                                                   aaGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTg                                  <  2:2972965/150‑1 (MQ=255)

                                                                                                                    aGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGc                                 >  1:2216835/1‑150 (MQ=255)

                                                                                                                      gCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCa                               >  2:466020/1‑150 (MQ=255)

                                                                                                                        aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                          <  2:2180961/57‑1 (MQ=255)

                                                                                                                        aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                          >  1:2180961/1‑57 (MQ=255)

                                                                                                                        aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAg                             >  1:5202996/1‑150 (MQ=255)

                                                                                                                        aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAg                             <  1:4628122/150‑1 (MQ=255)

                                                                                                                             gTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccgg                        <  2:1315697/150‑1 (MQ=255)

                                                                                                                               tCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccggcc                      <  1:466020/150‑1 (MQ=255)

                                                                                                                                   ccGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGccc                  >  1:3696362/1‑150 (MQ=255)

                                                                                                                                        gCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCActgc             <  2:1405263/150‑1 (MQ=255)

                                                                                                                                                 cGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTg    <  2:4194460/150‑1 (MQ=255)

                                                                                                                                                 cGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTg    <  2:5091704/150‑1 (MQ=255)

                                                                                                                                                  gCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGc   <  1:1212196/150‑1 (MQ=255)

AACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCC  >  NC_002947/1550099‑1550395

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