BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,928	T→A	42.0%	intergenic (‑1011/‑157)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,928	0	T	A	42.0%	73.6 / 151.6	127	intergenic (‑1011/‑157)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (34/39); new base A (35/18); total (70/57)
Fisher's exact test for biased strand distribution p-value = 4.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 128 total aligned reads displayed.

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATAC  >  CP000730/1929790‑1930060

tataTTTGCCGATGAAATAATAGGAGAGTTGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGtttttt                                                                                                                                            <  1:261145/140‑1 (MQ=255)

    ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                                       >  1:378478/7‑140 (MQ=255)

    ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                                       >  1:444181/7‑140 (MQ=255)

     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                      >  2:223806/6‑140 (MQ=255)

     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TGACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                      <  2:226637/135‑1 (MQ=255)

      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGATGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                     >  1:297137/5‑140 (MQ=255)

       acagatagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                    <  2:152124/132‑1 (MQ=255)

        cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAAAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGGTAACAACATAGGACCAGTCATCAGTGGCTGTGCTGTTTCGTATTATACTTATTgg                                                                                                                                   >  1:337951/3‑140 (MQ=255)

         agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                  <  2:36775/139‑1 (MQ=255)

         agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                  <  2:383188/139‑1 (MQ=255)

          gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                  <  1:437332/140‑1 (MQ=255)

          gatgatATAATAGTAAAGTCGCCTATCTCTCAGGCGTCAATTTATACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                  <  1:110705/140‑1 (MQ=255)

             ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  1:426444/138‑4 (MQ=255)

             ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTACTATACTTATTGGCTTAg                                                                                                                              <  2:335683/138‑1 (MQ=255)

              atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                              <  2:495638/140‑1 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                            >  2:558378/1‑140 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAA‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCATTCATCAGTGGCTTTTCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                            >  2:157312/1‑140 (MQ=255)

                ttaataGTAGAGCCGCTGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                            >  2:146383/2‑140 (MQ=255)

                cgatctatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTGACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                           >  2:6591/8‑139 (MQ=255)

                ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTCTTGGcggagaat                                                                                                                           >  2:135965/1‑133 (MQ=255)

                ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                           >  1:374131/1‑139 (MQ=255)

                 taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATAATAGCACCAGTTATCAGTGGCTGTGCCAATGCATTTTTTTCTTATTGGCTAAGTAg                                                                                                                           <  1:28062/140‑1 (MQ=255)

                  aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                         <  1:540592/140‑4 (MQ=255)

                  aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                         >  1:667840/1‑140 (MQ=255)

                   ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCT‑TTA‑G‑T‑TCTTAAAGCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                        <  1:154182/140‑1 (MQ=255)

                   ataCTAGAGTAGCCTATCTCTCAGGCAAAAATT‑‑GAAGAAGAGAGGAGGTGAATTAA‑T‑T‑TGAA‑AAGCATTTTCAGTTGCATCAGATCACCAGTCATCAGTGGCTGTGGTGTTGCGTATGATACTTATTGGCTTAGTAAACg                                                                                                                        <  2:657606/140‑1 (MQ=255)

                     agtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTACTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                       <  1:116665/140‑1 (MQ=255)

                       tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                    >  2:353430/1‑139 (MQ=255)

                          agacGCCTATCTCTCAGGCGACAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                 <  1:179096/137‑5 (MQ=255)

                          agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCATTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                                  >  2:371126/1‑140 (MQ=255)

                          agTCGCCTACCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                                  <  1:173036/140‑1 (MQ=255)

                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                >  1:631400/1‑135 (MQ=255)

                           gTCGCCTATATCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                >  2:251746/1‑135 (MQ=255)

                            tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa                                                                                                              >  1:284975/1‑140 (MQ=255)

                            tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑TCTCAATTGCCGCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                               <  1:484475/140‑7 (MQ=255)

                            tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                               >  2:540239/1‑134 (MQ=255)

                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACTCAATAC‑‑aa‑‑aa                                                                                                          >  1:467819/1‑140 (MQ=255)

                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa‑‑aa                                                                                                          <  2:177297/140‑1 (MQ=255)

                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATAATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                             >  1:186500/1‑132 (MQ=255)

                               cctctctttacggCGTAAATTGAGACCCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGTT‑GGTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGCAATAC‑‑AA‑‑AAt                                                                                                         <  2:179563/129‑1 (MQ=255)

                               ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCATATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AAt                                                                                                         >  1:264990/1‑140 (MQ=255)

                               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                            >  1:30862/1‑131 (MQ=255)

                               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                            >  2:324772/1‑131 (MQ=255)

                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           >  2:100951/1‑130 (MQ=255)

                                cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCCGTCGATGTGCTGTTGCGTATTATACGTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           <  2:304254/140‑11 (MQ=255)

                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTCC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTAGCTTAGTAAACACaacaagtaattta                                                                                                          >  1:649359/1‑129 (MQ=255)

                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGATGTTGCGTATAATACTTATTGGCTTGGTAAACGCaacaaataattaa                                                                                                          >  1:451416/1‑129 (MQ=255)

                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                          >  1:307337/1‑129 (MQ=255)

                                  aTCTCTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                         <  2:214678/140‑13 (MQ=255)

                                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattact                                                                                                         >  2:624595/1‑128 (MQ=255)

                                    ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                       >  2:607649/1‑126 (MQ=255)

                                    ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                       <  2:660952/140‑15 (MQ=255)

                                         aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ata                                                                                              >  2:53056/1‑140 (MQ=255)

                                              tccaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGcc                                                                                        <  2:654316/137‑1 (MQ=255)

                                                aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                     >  2:298503/1‑140 (MQ=255)

                                                 aTTTAGGCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                    <  2:117695/140‑1 (MQ=255)

                                                   ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAcc                                                                                  >  2:355254/1‑140 (MQ=255)

                                                         gCAGATACGAGTTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTCGCCGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                            <  1:657444/140‑1 (MQ=255)

                                                           agagagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaaa                                                                          <  1:546191/140‑1 (MQ=255)

                                                            gagagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GAACCAATAAAAAt                                                                         >  1:373733/1‑140 (MQ=255)

                                                             agagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                        >  1:442154/1‑140 (MQ=255)

                                                             agagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                        <  1:581282/140‑1 (MQ=255)

                                                               gggaggGGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCCAAAAAAATcac                                                                      >  2:541456/2‑138 (MQ=255)

                                                               aggaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAAACATAAAAATccc                                                                      >  2:228734/1‑140 (MQ=255)

                                                                ggaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATcccc                                                                     <  1:420496/140‑1 (MQ=255)

                                                                 gaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                    <  2:375683/140‑1 (MQ=255)

                                                                   ggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGAGGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCACTAAAAATCCCCCCa                                                                  >  2:293383/1‑140 (MQ=255)

                                                                    gtgtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                 >  1:556422/1‑140 (MQ=255)

                                                                     tgtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACt                                                                <  1:208517/140‑1 (MQ=255)

                                                                      gtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                               >  1:634105/1‑140 (MQ=255)

                                                                      gtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGAAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCACCTCACTa                                                               >  1:433186/1‑140 (MQ=255)

                                                                      gtA‑TAA‑G‑G‑TGAT‑GCTTATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGAAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                               >  1:354392/1‑140 (MQ=255)

                                                                            a‑G‑G‑TGAT‑GCTTATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑CCCCGGTCAC‑AGAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                          <  2:607591/140‑1 (MQ=255)

                                                                                g‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                        >  2:255216/1‑140 (MQ=255)

                                                                                  tgtt‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcggt                                                       <  1:601680/137‑9 (MQ=255)

                                                                                   taT‑GCTTATTTTTGTTCACATCATAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATa                                                      <  1:301929/139‑1 (MQ=255)

                                                                                        cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa                                                  >  1:60442/1‑140 (MQ=255)

                                                                                         ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGCC‑ATAT‑AGCC‑GCACCAATACACAACCCCTCACTACCGCAAGTaacgcg                                                 >  2:552250/1‑135 (MQ=255)

                                                                                         ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATTAAAATCCCCTCACTACCACAAagagagag                                                 >  2:424809/1‑140 (MQ=255)

                                                                                          tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                <  2:58478/140‑1 (MQ=255)

                                                                                          tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGGTTTGTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                >  1:457144/1‑140 (MQ=255)

                                                                                            ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                              <  1:130121/140‑1 (MQ=255)

                                                                                                  ttCACATCATAGCACCAATCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                        >  2:25765/1‑140 (MQ=255)

                                                                                                        tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCCAAAAGTGAGGGGATTGGTGTATaa                                  <  2:643649/140‑1 (MQ=255)

                                                                                                        tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGAATaa                                  >  2:665917/1‑140 (MQ=255)

                                                                                                         catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                 <  1:146383/140‑1 (MQ=255)

                                                                                                           tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGgtgcat                                    <  1:102347/140‑4 (MQ=255)

                                                                                                               aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                           >  1:185855/1‑140 (MQ=255)

                                                                                                                     catcaGTGGCTGCGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAT‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTAGCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                     <  1:53056/140‑1 (MQ=255)

                                                                                                                      atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                    <  1:240119/140‑1 (MQ=255)

                                                                                                                      atcaGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTATCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                    <  1:343325/140‑1 (MQ=255)

                                                                                                                       tcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTATGTAAACGTAATAC‑‑TA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAC‑TAttttt                  <  2:467819/140‑2 (MQ=255)

                                                                                                                        caGTGGCTGTGCCATTGCGTTTTTTTTTTGTTGGGTGGGAAGACGCAATAC‑‑AA‑‑AAGAGGAGAC‑GTAT‑AGCC‑GCACCAATGAAAATCCCCTCACTACCGCAAATAGTAGGGGGATTGGTGTATAAGGAAATACTTTTTTTc                  >  2:341920/1‑140 (MQ=255)

                                                                                                                        caGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑CA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTc                  >  1:221683/1‑140 (MQ=255)

                                                                                                                               tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATACGGAGGGGGTAGGTGTATAAGTAAATACTTTTTTTTGTTTTCt           >  2:370638/1‑140 (MQ=255)

                                                                                                                                 tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAACACTTATTTTCGTTGTCTTa         <  1:371126/140‑1 (MQ=255)

                                                                                                                                  gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTaa        <  1:637807/140‑1 (MQ=255)

                                                                                                                                    cATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAAAAAAAACCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAAACTTATTTTCGTTGTCTTAAtt      >  1:596654/1‑140 (MQ=255)

                                                                                                                                     aTTGCGTTTTTTTCGTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTTCCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATta     <  1:180334/140‑1 (MQ=255)

                                                                                                                                        gCGTTTTTTTCTCTTTGGCTAAGTAGACTCAATAC‑‑AA‑‑AATAGGGCAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTTCCGTAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTTTCTTAATTATAc  <  2:43113/140‑1 (MQ=255)

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATAC  >  CP000730/1929790‑1930060

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