BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,587,225	A→T	5.8%	M95L (ATG→TTG)	yhhA →	DUF2756 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,225	0	A	T	5.8%	472.4 / 15.5	137	M95L (ATG→TTG)	yhhA	DUF2756 family protein
Reads supporting (aligned to +/- strand): ref base A (59/70); new base T (3/5); total (62/75)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

Only 100 of 140 total aligned reads displayed.

AGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTA  >  NC_000913/3587079‑3587372

  cAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgct                                                                                                                                                 <  3:222303/149‑1 (MQ=255)

   aaaTCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                <  4:832625/149‑1 (MQ=255)

      tCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtg                                                                                                                                                   >  4:142499/1‑143 (MQ=255)

       ccAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccgg                                                                                                                                            <  3:693425/149‑3 (MQ=255)

        cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCaa                                                                                                                                           >  2:1009958/1‑149 (MQ=255)

        cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCaa                                                                                                                                           >  1:579818/1‑149 (MQ=255)

         aGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccgggc                                                                                                                                          <  2:481517/149‑5 (MQ=255)

         aGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAac                                                                                                                                          <  2:708824/149‑1 (MQ=255)

         aGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAac                                                                                                                                          <  1:56589/149‑1 (MQ=255)

          gCAAAAAGGGATGCTGAATCAGTAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGCGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAaca                                                                                                                                         <  1:945817/149‑1 (MQ=255)

          gCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAaca                                                                                                                                         <  3:168106/149‑1 (MQ=255)

            aaaaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACAcc                                                                                                                                       >  1:501871/1‑149 (MQ=255)

            aaaaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACAcc                                                                                                                                       >  2:807103/1‑149 (MQ=255)

             aaaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGAAAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCa                                                                                                                                      >  1:809241/1‑149 (MQ=255)

              aaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCaa                                                                                                                                     >  3:322638/1‑149 (MQ=255)

               aaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAAc                                                                                                                                    <  2:501871/149‑1 (MQ=255)

                     tGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCggg                                                                                                                              <  3:129686/149‑1 (MQ=255)

                         gAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGt                                                                                                                          <  3:118357/149‑1 (MQ=255)

                         gAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGt                                                                                                                          <  4:322638/149‑1 (MQ=255)

                          aaTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGtt                                                                                                                         <  1:1009958/149‑1 (MQ=255)

                          aaTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGtt                                                                                                                         >  3:433795/1‑149 (MQ=255)

                            tCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaa                                                                                                                       >  2:436123/1‑149 (MQ=255)

                             cAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaaa                                                                                                                      <  4:810028/149‑1 (MQ=255)

                             cAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaaa                                                                                                                      <  3:181598/149‑1 (MQ=255)

                              aGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccg                                                                                                                                             <  2:14522/125‑2 (MQ=255)

                              aGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccg                                                                                                                                             >  1:14522/1‑124 (MQ=255)

                              aGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAAc                                                                                                                     <  2:328205/149‑1 (MQ=255)

                               gCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACa                                                                                                                    <  1:511407/149‑1 (MQ=255)

                                 aaCTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGAt                                                                                                                            <  2:212543/139‑1 (MQ=255)

                                 aaCTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGc                                                                                                                  >  3:453470/1‑149 (MQ=255)

                                    tGAAAACGCAAACTCAGTTGCATCAGCAGCACTGAGAAAACCAGATAAACAATAATTCTCATCGAGTGTCGCAATCAGAGCCGGGGAAGCGAAATCCCACCCGCCAGCAAATGCTGCCCAAAACCAACGTCGGGATGTTAaacaacaac                                                                                                               >  3:973403/1‑149 (MQ=255)

                                    tGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAAc                                                                                                               <  1:225965/149‑1 (MQ=255)

                                    tGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAAc                                                                                                               >  4:900222/1‑149 (MQ=255)

                                     gAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAAcc                                                                                                              <  2:1029570/149‑1 (MQ=255)

                                     gAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAAcc                                                                                                              <  4:234142/149‑1 (MQ=255)

                                      aaaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCg                                                                                                             >  1:677067/1‑149 (MQ=255)

                                        aaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           >  3:1058784/1‑149 (MQ=255)

                                        aaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           <  1:312840/149‑1 (MQ=255)

                                         aCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTaa                                                                                                          >  2:259870/1‑149 (MQ=255)

                                            cAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATcc                                                                                                       <  3:375270/149‑1 (MQ=255)

                                             aaaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCg                                                                                                      >  2:720129/1‑149 (MQ=255)

                                              aaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGccc                                                                                                                                             >  3:49968/1‑109 (MQ=255)

                                              aaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGccc                                                                                                                                             <  4:49968/109‑1 (MQ=255)

                                               aCTCAGTTGCAACAGCAGCATTTAGAGAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    <  1:1069318/149‑1 (MQ=255)

                                               aCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    <  1:679915/149‑1 (MQ=255)

                                               aCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    <  4:402264/149‑1 (MQ=255)

                                                cTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGAt                                                                                                   <  3:376586/149‑1 (MQ=255)

                                                     ttGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTc                                                                                              <  4:1071721/149‑1 (MQ=255)

                                                       gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCTCAGCCGGGGGAGCGTAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAACCCGGATAGTTCgc                                                                                            <  4:936315/149‑2 (MQ=255)

                                                       gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaaa                                                                                                                      >  3:21125/1‑123 (MQ=255)

                                                       gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaaa                                                                                                                      <  4:21125/123‑1 (MQ=255)

                                                       gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGt                                                                                            <  3:435439/149‑1 (MQ=255)

                                                        cAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGtt                                                                                           <  1:265806/149‑1 (MQ=255)

                                                                 caTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGtt                                                                                               <  4:207350/136‑1 (MQ=255)

                                                                 caTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                                  <  1:259870/149‑1 (MQ=255)

                                                                 caTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                                  >  1:1070656/1‑149 (MQ=255)

                                                                  aTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagcag                                                                                 >  4:246667/1‑149 (MQ=255)

                                                                   tttAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCaa                                                                                                                                     >  4:261878/1‑96 (MQ=255)

                                                                   tttAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGaa                                                                                        >  2:677197/1‑141 (MQ=255)

                                                                          aaCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTc                                                                                              >  1:658738/1‑128 (MQ=255)

                                                                          aaCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgct                                                                         >  1:325060/1‑149 (MQ=255)

                                                                             cAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGcc                                                                      <  3:891982/149‑1 (MQ=255)

                                                                               gATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACATCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCgg                                                                    >  4:640184/1‑149 (MQ=255)

                                                                               gATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCgg                                                                    >  3:354226/1‑149 (MQ=255)

                                                                                  aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 >  1:756443/1‑149 (MQ=255)

                                                                                  aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 <  1:918618/149‑1 (MQ=255)

                                                                                  aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAAGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 <  3:650123/149‑1 (MQ=255)

                                                                                     cAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagga                                                              <  3:640184/149‑1 (MQ=255)

                                                                                     cAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagga                                                              <  1:464676/149‑1 (MQ=255)

                                                                                      aataatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagaggag                                                             >  2:344055/1‑149 (MQ=255)

                                                                                      aataatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagaggag                                                             >  2:102323/1‑149 (MQ=255)

                                                                                        taatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGaa                                                           >  4:985035/1‑149 (MQ=255)

                                                                                        taatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGaa                                                           >  2:218879/1‑149 (MQ=255)

                                                                                         aatTCCCAGCGGGTGTTGCAGACACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCACAACACCAACGGCGGGATGTTAATCAGCAACCGTAATCCGCATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGaaa                                                          <  4:339279/149‑1 (MQ=255)

                                                                                             ctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCTGGAGAGGAGAAACGGc                                                      <  4:756465/149‑1 (MQ=255)

                                                                                             ctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGc                                                      <  2:11600/149‑1 (MQ=255)

                                                                                              tcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCg                                                     >  3:993247/1‑149 (MQ=255)

                                                                                               cAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGa                                                    >  3:332640/1‑149 (MQ=255)

                                                                                               cAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGa                                                    <  4:354226/149‑1 (MQ=255)

                                                                                                 gCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGAAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 >  3:1021778/1‑134 (MQ=255)

                                                                                                     gtgtTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCt                                              <  2:981647/149‑1 (MQ=255)

                                                                                                     gtgtTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCt                                              <  2:756443/149‑1 (MQ=255)

                                                                                                         tGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGTGACATGCTGAAt                                          <  4:964839/149‑1 (MQ=255)

                                                                                                         tGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAAt                                          <  1:203998/149‑1 (MQ=255)

                                                                                                         tGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAAt                                          <  2:1029124/149‑1 (MQ=255)

                                                                                                          gCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATc                                         >  2:38452/1‑149 (MQ=255)

                                                                                                          gCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATc                                         >  2:603817/1‑149 (MQ=255)

                                                                                                          gCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGGTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATc                                         >  3:1007968/1‑149 (MQ=255)

                                                                                                              tCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGcc                                     <  1:492142/149‑1 (MQ=255)

                                                                                                                 caGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAg                                  <  3:158958/149‑1 (MQ=255)

                                                                                                                  aGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGc                                 >  3:948350/1‑149 (MQ=255)

                                                                                                                               aaaTCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGAc                                                   >  1:83148/1‑118 (MQ=255)

                                                                                                                               aaaTCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTg                    >  1:800435/1‑149 (MQ=255)

                                                                                                                                  tCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACa                 >  3:769441/1‑149 (MQ=255)

                                                                                                                                    ccgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGTAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACAtt               >  1:991719/1‑149 (MQ=255)

                                                                                                                                         gggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          >  4:799839/2‑149 (MQ=255)

                                                                                                                                              gcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGaaaa     >  3:754180/1‑149 (MQ=255)

                                                                                                                                                aaaTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACt   >  2:8350/1‑149 (MQ=255)

                                                                                                                                                aaaTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACt   >  2:197017/1‑149 (MQ=255)

                                                                                                                                                 aaTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTa  >  4:424062/1‑149 (MQ=255)

AGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTA  >  NC_000913/3587079‑3587372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`